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ABSTRACT 


This  report  presents  the  results  of  a research  program  that  had  two 
major  objectives.  The  first  objective  was  the  development  of  a prototype 
rotating  stall  control  system  which  was  tested  both  on  a low  speed  rig  and  a 
J-S5-3  engine.  The  second  objective  was  to  perform  fundamental  studies  of  the 
flow  mechanisms  that  produce  rotating  stall,  surge  and  noise  in  axial  flow 
compressors  and  thereby  obtain  an  understanding  of  these  phenomena  that  would 
aid  attaining  the  first  objective.  The  work  is  reported  in  three  separate 
volumes.  Volume  I covers  the  fundamental  theoretical  and  experimental  studies 
of  rotating  stall;  Volume  II  covers  the  theoretical  and  experimental  studies  of 
discrete- tone  aerodynamic  noise  generation  mechanisms  in  axial  flow  compressors; 

and,  Volume  III  covers  the  development  and  testing  of  a prototype  rotating  stall 

\ 

control  system  on  both  the  low  speed  test  rig  and  the  J-85-5  engine. 

Volume  I describes  the  theoretical  and  experimental  investigation  of 
the  influence  of  distortion  on  the  inception  and  properties  of  rotating  stall 
for  an  isolated  rotor  row,  and  the  effects  of  close  coupling  of  a rotor  and 
stator  row  on  rotating  stall  inception.  The  experiments  were  conducted  in  the 
Calspan/Air  Force  Annular  Cascade  Facility,  which  is  a low  speed  compressor 
research  rig.  In  addition,  the  previously  developed  two  dimensional  stability 
theory  for  prediction  of  inception  conditions  was  extended  to  include  the 
effect  of  compressibility  and  the  development  of  a three  dimensional  theory 
was  initiated.  These  studies  led  to  the  following  key  results.  The  experi- 
mental studies  of  distortion  show  that  for  a single  blade  row  the  response  of 
the  blade  row  to  the  distortion  and  rotating  stall  are  uncoupled  phenomena  and 
may  be  explained  on  the  basis  of  a linearized  analysis.  The  experimental 
studies  of  a closely  coupled  rotor-stator  pair  show  that  the  addition  of  a 
closely  spaced  stator  row  downstream  of  a rotor  row  delays  the  onset  of  rotating 
stall.  Moreover,  the  corresponding  theoretical  analysis  predicts  this  trend 
although  quantitative  agreement  is  hampered  by  the  lack  of  appropriate  steady- 
state  loss  and  turning  performance  for  each  blade  row  at  the  required  operating 
conditions.  The  theoretical  investigation  of  the  effects  of  compressibility 
for  wholly  subsonic  flows  outside  the  blade  rows  indicates  that  the  effects  of 
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compressibility  do  not  alter  the  mechanisms  of  rotating  stall  as  deduced  from 
the  incompressible  theory  in  that  the  rate  of  change  of  the  steady  state  loss 
curve  with  inlet  swirl  is  the  dominant  blade  row  characteristic  affecting  its 
stability.  Therefore,  if  the  steady  state  losses  are  known  for  the  compressible 
flow  condition,  the  linearized  stability  analysis  is  expected  to  apply. 

Volume  II  describes  a theoretical  and  experimental  study  of  discrete- 
tone  noise  generation  by  the  interaction  of  a rotor  and  a stator,  and  the  de- 
velopment of  a direct  lifting  surface  theory  for  an  isolated  rotor.  An  approxi- 
mate model  has  been  developed  to  predict  the  sound  pressure  level  and  total 
power  radiated  at  harmonics  of  the  blade  passage  frequency  for  a rotor-stator 
stage.  The  analysis  matches  the  duct  acoustic  modes  for  an  annular  duct  with 
an  approximate  representation  of  the  unsteady  blade  forces  which  includes  com- 
pressibility effects.  Measurements  were  made  of  the  sound  pressure  levels 
produced  on  the  duct  wall  of  the  annular  cascade  facility  by  a rotor-stator 
pair.  Predictions  which  indicated  that  only  the  fourth  and  higher  harmonics 
could  be  excited  at  conditions  achievable  in  the  facility,  were  borne  out  by 
the  experiments.  The  calculations  of  the  sound  pressure  levels  for  the  propa- 
gating modes  were  significantly  below  the  measured  values.  This  discrepancy 
is  believed  to  result  from  inaccuracies  in  existing  models  of  rotor  wake 
velocity  profiles,  which  are  shown  to  have  a strong  influence  on  predictions 
of  the  sound  pressure  levels  of  the  higher  harmonics.  Volume  II  also  contains 
the  formulation  of  a direct  lifting  surface  theory  for  the  compressible, 
three-dimensional  flow  through  a rotor  row  in  an  infinitely  long  annular  duct. 

A detailed  derivation  is  given  for  the  linearized  equations  and  the  corres- 
ponding solutions  for  the  blade  thickness  and  loading  contributions  to  the 
rotor  flow  field.  The  governing  integral  equation  for  the  blade  loading  in  a 
lifting  surface  theory  is  obtained  for  subsonic  flow  and  progress  on  its 
solution  is  reported. 

Volume  III  describes  the  development  and  testing  of  a prototype  ro- 
tating stall  control  system.  The  control  system  was  tested  on  the  low  speed 
compressor  research  rig  and  on  a J -85-5  turbojet  engine.  On  the  low  speed 
research  compressor,  the  control  was  tested  in  the  presence  of  circumferential 
inlet  distortion.  These  tests  were  performed  to  demonstrate  the  ability  of 
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the  control  to  operate  satisfactorily  in  the  presence  of  inlet  distortion  and 
to  aid  in  the  selection  of  stall  sensor  configurations  for  the  subsequent 
engine  tests.  The  control  system  was  then  installed  on  a J-85-5  jet  engine  and 
its  performance  was  tested  under  sea  level  static  conditions,  both  with  and 
without  inlet  distortion.  On  the  engine,  the  stall  control  was  installed  to 
override  the  normal  operating  schedule  of  the  compressor  bleed  doors  and  inlet 
guide  vanes.  The  J-85-5  was  stalled  in  two  ways,  first  by  closing  the  bleed 
doors  at  constant  engine  speed,  and  second  by  decelerating  the  engine  with  the 
bleed  doors  partially  closed  at  the  beginning  of  the  deceleration.  A total  of 
41  compressor  stalls  were  recorded  at  corrected  engine  speeds  between  48  and 
72  percent  of  the  rated  speed.  In  all  cases,  the  control  took  successful 
remedial  action  which  limited  the  duration  of  the  stall  to  325  milliseconds 
or  less.  "v 
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P (1-MXV/2 

dfx)  Dirac  delta  function 

if  Fourier  transform  variable  in  axial  coordinate 
6 azimuthal  coordinate 

/ rT 

fundamental  blade  loading  frequency^  Equation  (12) 
po  undisturbed  density 

O'  dimensionless  radius,  r/>r 

14,  <A  . 

Tnrn  propagation  angle  of  (rt  nr)  mode,  relative  to  dipole  axis.  Figure  14 
jfl  rotor  angular  velocity,  radians/sec. 
uj  acoustic  frequency,  N^BD. 

< > denotes  time  average  over  one  blade-passing  period 
[ j magnitude  of  a vector  or  complex  number 
* complex  conjugate 


XV 


SYMBOLS,  SECTION  III  (Cont'd.) 


Subscripts 

1,1  refer  to  upstream  or  downstream  blade  row,  respectively 
C evaluated  on  viscous  wake  centerline 
2 referes  to  2 th  loading  harmonic 
m refers  to  m th  radial  mode 
n refers  to  n th  azimuthal  mode 

R,5  refers  to  rotor  or  stator  blade  rows,  respectively 

Superscripts 

u,cL  refer  to  upstream  or  downstream  propagating  waves,  respectively 

SYMBOLS  FOR  SECTION  IV,  DIRECT  LIFTING  SURFACE  THEORY  FOR  A COMPRESSOR  ROTOR 
undisturbed  sound  speed 
R vector  defined  in  Equation  (63) 

FI i coefficient  in  polynomial  expansion  for  axial  variation  of  blade  loading 

U) 

Ryr>e-ft  real  Part  of  coefficient  defined  in  Equation  (173),  see  Table  3 
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6 

(V 

C 

C-o. 

ra) 

d,,dz 

D 

r>eQ’gi 

f 

F 

7 

? 

G 


number  of  blades 

imaginary  part  of  coefficient  defined  in  Equation  (173),  see  Table  3 
blade  chord 

axial  projection  of  blade  chord 

coefficient  defined  in  Equation  (173) , see  Table  3 
integrals  defined  in  Equation  (C-SO) , Appendix  C. 
dipole  strength 

quantity  defined  in  Equation  (167) 

unit  vectors  in  cylindrical  coordinate  system 

function  of  radius  appearing  in  blade  thickness  distribution. 
Equation  (110) 

force  per  unit  volume 

net  force  on  control  volume 

axial  dependence  of  blade  thickness  distribution.  Equation  (110) 
Green's  function 

quantity  defined  in  Equation  (178) 
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A ratio  c.V  hub  radius  to  tip  radius 
H ( ) Heaviside  step  function 
I( , lz  surface  integrals  in  Equation  (66) 

1(2)  integral  of  pressure  perturbation  along  streamlines 
integrals  defined  in  Appendix  C 
■k  radial  mode  number 

radial  eigenvalues 

L operator  defined  by  Equation  (58) 

t-mBk  quantity  defined  in  Equation  (173) 

►■vt  multiple  of  blade  number  for  azimuthal  mode,  n = m 8 

oo  mass  flux  to  first  order  in  perturbation  quantities 
M Mach  number  of  undisturbed  axial  stream,  D/cl# 

Mr  Mach  number  based  on  undisturbed  relative  velocity  in  blade  fixed 

coordinates,  Ua  / dt 

OO 

Mt  Mach  number  based  on  rotor  speed  at  blade  tip  radius,  oorr/aa, 

a-  azimuthal  mode  number  or  coordinate  mutually  normal  to  undisturbed 

streamline  and  radial  directions 
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unit  vector  normal  to  undisturbed  streamline  and  radial  directions 
blade  camber  line 

distance  to  upper  and  lower  blade  surfaces 
pressure  perturbation 
undisturbed  pressure 

terms  in  expression  for  pressure  perturbation,  Equation  (148) 

V 

local  blade  loading 
blade  loading  averaged  over  radius" 

terms  in  polynomial  representation  of  blade  loading.  Equation  (190) 
source  strength 

integrals  defined  by  Equation  (C-39) , Appendix  C 
radial  coordinate 
normalized  radial  eigenfunctions 
distance  along  undisturbed  streamlines 
unit  vector  along  undisturbed  streamlines 
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Acjrf'C  ) +1  for  arg  > 0,  -1  for  arg  < 0 

S surface  area 

S(r,Z)  function  defined  in  Equation  (121) 

t blade  thickness 

TmS-k  quantity  defined  in  Equation  (131) 

x7z  axial  component  of  torque  exerted  on  control  volume 

U undisturbed  axial  velocity 

UR  undisturbed  relative  velocity  in  blade-fixed  coordinates 

V perturbation  velocity,  (vr  , vQ  , V£)  or  C vr  , vn  , vs  ) 

l terms  in  normal  component  of  perturbation  velocity,  Equations 
and  (168) 

V volume 

quantity  defined  in  Equation  (117) 

-X 

W total  velocity  in  blade  fixed  coordinates 

X normalized  axial  coordinate,  Equation  (179) 

2 axial  coordinate 


(160) 


xx 


SYMBOLS,  SECTION  IV  (Cont'd.) 


p vW 

2f  local  vortex  strength 

r blade  circulation 

<$(  ) Dirac  delta  function 

6nr^  Kronecker  delta 

£ vanishly  small  quantity 

€ helical  coordinate,  if  = 6 - ujZ/U 

generalized  function  defined  in  Equation  (126) 

Ca./2rT 

& azimuthal  coordinate 

quantity  defined  in  Equation  (71) 
quantity  defined  in  Equation  (114) 
jj.  index  of  Lommel  function 

i)  normal  to  surface 

% Fourier  transform  variable,  or  dummy  variable  for  argument  of  Bessel 
functions 
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poles  of  integrand  in  Equation  (78) 


p density  perturbation 


p m undisturbed  density 


(T  normalized  radial  coordinate,  r l rr 


d)  velocity  potential 


dj n£  coefficient  in  expansion  of  $ in  tei'ms  of  duct  eigenfunctions. 
Equation  (72) 


$ scalor  function  defined  in  Equation  (62) 


$r  cor-T  / U 


y ' (ojr/  U ) 


lu  angular  velocity  of  rotor 


( ) Fourier  transform 


C ) dimensionless  variable,  see  Equation  (186) 


Subscripts 


D property  of  dipole  solution 


H evaluated  at  hub  radius 


SYMBOLS,  StCTLON  IV  (Cont'd.) 


j.  th  blade 
~k  th  radial  mode 
>i  th  azimuthal  mode 
source  coordinates 
property  of  source  solution 
evaluated  at  tip  radius 
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SECTION  I 


INTRODUCTION 


The  useful  operating  range  of  a turbine  engine  compressor  is  greatly 
influenced  by  its  stalling  characteristics.  The  optimum  performance  of  a 
turbo-propulsion  system  is  usually  achieved  when  the  compressor  is  operating 
near  its  maximum  pressure  ratio.  However,  this  optimum  is  generally  not 
attainable  because  it  occurs  close  to  compressor  stall  aid  unstable  flow  con- 
ditions. Because  of  the  serious  mechanical  damage  that  may  result  during  com- 
pressor stall  cycles,  a factor  of  safety  (stall  margin)  must  be  provided 
between  the  compressor  operating  line  and  the  stall  boundary.  This  is  usually 
done  by  prescheduling  the  primary  engine  controls.  However,  the  prescheduling 
approach  can  lead  to  the  requirement  for  a large  stall  margin  in  order  to  keep 
the  engine  from  stalling  under  all  possible  transient  and  steady  state  flight 
conditions.  It  is  clear,  then,  that  an  engine  control  system  that  can  sense 
incipient  destructive  unsteady  flow  in  a compressor  and  take  corrective  action 
would  allow  for  reduced  stall  margins  in  the  design  and  thus  lead  t.o  large 
engine  performance  and/or  efficiency  gains.  Recognition  of  this  fact  has  been 
the  motivation  for  a continuing  program  of  research  that  the  AFAPL  has  sponsored 
at  Cal  span  dating  back  to  1962. 

The  work  at  Calspan  has  been  both  theoretical  and  experimental  in 
nature  and  has  been  aimed  at  obtaining  a sufficient,  understanding  of  the  rotat- 
ing stall  phenomena  such  that  its  onset  and  its  properties  can  be  predicted 
and  controlled.  The  capability  of  predicting  the  onset  of  rotating  stall  on 
isolated  blade  rows  of  high  hub  to  tip  ratios  in  low  speed  flows  was  demon- 
strated in  Reference  1,  In  addition,  the  basic  feasibility  of  developing  a 
rotating  stall  control  system  was  demonstrated  in  the  Calspan/Air  Force  Annular 
Cascade  Facility.  This  present  report  summarizes  the  latest  three  year 
research  program  at  Calspan.  The  specific  goals  of  the  present  program  were 
to  extend  the  fundamental  studies  of  rotating  stall  to  consider  the  effects  of 
compressibility,  blade  row  interaction  and  inlet  distortion;  and  to  extend  the 


1 


fundamental  aerodynamic  and  acoustic  analysis  of  flow  through  a compressor. 

In  addition,  the  rotating  stall  control  system  was  validated  by  successful 
ground  tests  on  a J-85-5  turbojet  engine. 

The  work  is  reported  in  three  separate  volumes.  Volume  I entitled, 
"Basic  Studies  of  Rotating  Stall",  covers  the  theoretical  and  experimental 
work  on  the  effects  of  distortion  and  close  coupling  of  blade  rows  on  rotating 
stall  inception  and  properties.  In  addition,  the  theoretical  analysis  of  com- 
pressibility is  treated  in  the  two-dimensional  approximation  and  the  initial 
development  of  a three-dimensional  theory  is  given.  Volume  II  entitled,  "In- 
vestigation of  Rotor-Stator  Interaction  Noise  and  Lifting  Surface  Theory  for 
a Rotor",  describes  the  development  of  a linearized  lifting  surface  theory  for 
the  subsonic  compressible  flow  through  an  isolated  rotor  row.  In  addition,  a 
theoretical  and  experimental  study  of  the  noise  generated  by  the  interaction 
of  a rotor  and  stator  is  described.  Volume  III  entitled,  "Development  of  a 
Rotating  Stall  Control  System",  describes  the  development  and  testing  of  the 
control  system  installed  on  a low  speed  research  compressor  and  on  a J-85-5 
turbojet  engine. 

In  the  current  three-year  segment  of  the  ongoing  research  program  on 
rotating  stall,  the  scope  of  the  investigation  was  expanded  to  include  studies 
of  the  aerodynamics  and  acoustics  of  axial  flow  compressors.  Volume  II  con- 
tains the  results  of  this  aspect  of  the  program,  which  consisted  of  two  basic 
parts:  the  theoretical  and  experimental  investigation  of  discrete-tone  noise 

generation  by  the  interaction  of  a rotor  and  stator,  and  the  development  of 
the  three-dimensional  direct  lifting  surface  theory  for  a compressor  rotor. 

An  approximate  model  has  been  developed  to  predict  the  sound  pressure  level 
and  total  power  radiated  at  harmonics  of  the  blade  passage  frequency  for  a 
rotor-stator  stage.  Also,  measurements  were  made  of  the  sound  pressure  levels 
produced  on  the  outer  duct  wall  of  the  annular  cascade  facility  by  a rotor- 
stator  pair.  In  the  development  of  a direct  lifting  surface  theory,  the 
governing  integral  equation  relating  the  rotor  blade  loading  to  prescribed 
incidence  and  camber  lines  has  been  formulated  and  progress  made  toward  its 
numerical  solution.  The  experimental  and  theoretical  studies  of  rotor-stator 
interaction  noise  are  reported  in  Sections  II  and  III,  respectively.  The 
direct  lifting  surface  analysis  is  presented  in  Section  IV. 
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SECTION  II 

EXPERIMENTAL  ACOUSTICS  RESEARCH 


As  a part  of  the  work  under  a previous  program,  Contract  AF  33(615) -3357, 
an  annular  cascade  facility  was  designed  and  fabricated.  Its  principal  purpose 
is  to  provide  fundamental  experimental  data  during  and  prior  to  the  occurrence 
of  rotating  stall  in  order  to  improve  our  understanding  of  the  phenomena  and 
for  use  as  a guide  in  improving  the  theoretical  analysis.  In  addition  to  the 
study  of  rotating  stall,  the  facility  has  also  been  used  to  provide  acoustic 
data  for  comparison  with  theory  and  to  evaluate  the  operation  of  a prototype 
rotating  stall  control  system.  The  fundamental  experiments  on  rotating  stall 
are  described  in  Volume  I of  this  report  and  the  control  system  tests  are 
described  in  Volume  III.  This  section  of  Volume  II  presents  the  results  of 
the  acoustic  experiments.  x 

* 

Two  sets  of  experiments  were  performed  in  support  of  the  theoretical 
developments  described  in  Section  III  and  IV  of  this  volume.  The  first  set  of 
experiments  were  designed  to  provide  data  for  correlation  with  the  approximate 
theory  of  Section  III  for  the  prediction  of  rotor-stator  interaction  noise. 

In  these  experiments,  the  far-field  sound  pressure  levels  generated  by  roto- 
stator  interaction  were  measured  in  the  constant  area  annulus  upstream  of  a 
rotor-stator  stage.  These  measurements  were  performed  over  a range  of  rotor 
speeds  for  two  different  stagger  angle  settings  of  the  stator  blades.  The 
second  set  of  experiments  were  intended  to  provide  a measure  of  the  tip  pressure 
loading  on  an  isolated  rotor  for  comparison  with  the  direct  lifting  surface 
theory  of  Section  IV.  In  these  experiments,  time -varying  records  of  the  outer 
wall  static  pressure  fluctuations  caused  by  blade  tip  passage  were  obtained 
for  a number  of  different  chordwise  locations  on  the  rotor  blades.  The  scope 
of  these  tip  loading  experiments  was  very  limited. 

The  order  of  presentation  of  the  experimental  studies  in  the  remainder 
of  this  section  is  as  follows.  A brief  description  of  the  annular  cascade 
facility  is  presented  in  Section  II-A.  Section  I I - B presents  the  results  of 


the  measurements  of  noise  caused,  by  rotor-stator  interaction.  This  is  followed 
by  the  presentation,  in  Section  II-C,  of  a typical  set  of  results  obtained  for 
the  blade  passage  fluctuations  on  an  isolated  rotor.  Finally,  some  concluding 
remarks  on  the  experimental  studies  are  presented  in  Section  II-D. 

A.  DESCRIPTION  OF  ANNULAR  CASCADE  FACILITY 


A detailed  description  of  the  annular  cascade  facility  has  been  pre- 
sented in  Reference  1,  and  further  details  are  given  in  Volume  I of  this 
report.  Only  a brief  description  of  the  facility  is  given  here. 

The  annular  cascade  facility  consists  of  a test  section  built  around 
the  outer  front  casing  of  a J-79  jet-engine  compressor  with  a Calspan  fabri- 
cated  hub.  The  facility  includes  a bell-mouth  inlet  on  the  outer  casing  and  a 
bullet  nose  on  the  hub  to  provide  a smooth  flow  of  air  to  the  test  section. 
Outlet  ducting  is  connected  to  an  independently  variable  source  of  suction  to 
provide  the  required  flow  through  the  annulus.  An  electrically  powered  two- 
speed  axial  flow  fan  is  used  as  the  source  of  suction.  Continuous  control  of 
the  mass  flow  is  achieved  through  the  use  of  variable  inlet  guide  vanes  to  the 
fan  and  a variable  damper  in  the  fan  exit  flow. 

The  test  section  of  the  annular  cascade  forms  a circular  annulus  with 
an  outer  diameter  of  29.35  inches  and  an  inner  diameter  of  23.35  inches  which 
provides  a hub-to-tip  ratio  of  0.80.  The  outer  casing  will  accept  up  to  six 
variable  stagger  angle  stator  rows.  The  hub  has  provision  for  two  rotor  rows 
at  the  third-  and  fifth-stage  rotor  locations  of  the  J-79  compressor.  At  the 
time  the  acoustic  studies  were  performed,  either  rotor  row  could  be  rotated 
while  the  other  was  held  fixed  or  both  rows  could  be  rotated  together.  The 
rotor  assembly  was  powered  by  a 24  horsepower  hydraulic  motor.  Rotational 
speed  was  continuously  variable  in  either  direction  between  zero  and  approxi- 
mately 1500  rpm.  An  external  hydraulic  pump  system  powered  by  a 30  horsepower 
electric  motor  was  used  to  provide  power  for  the  hydraulic  motor. 
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Although  the  speed  range  of  the  rotor  was  restricted  at  the  time  of 
these  tests,  it  was  possible  to  perform  acoustic  tests  for  comparison  with  the 
theory  of  Section  III.  The  flexibility  of  control  and  the  relative  absence  of 
background  noise  in  the  facility  are  features  which  were  highly  useful  in  the 
performance  of  the  tests. 

B.  NOISE  GENERATED  BY  ROTOR-STATOR  INTERACTION 

The  configuration  of  the  annular  cascade  used  for  the  rotor-stator 
interaction  studies  is  shown  in  Figure  1 along  with  some  details  of  the  micro- 
phone-probe assembly  used  in  the  acoustic  work.  The  rotor-stator  stage 
studied  is  designated  Rotor  Set  No.  1 (46  blades)  and  Stator  Set  No.  1 (54 
vanes)  in  Reference  1.  These  are  modified  blade  rows  from  the  fifth  stage  of 
a J-79  compressor.  Their  characteristics  and  performance  when  used  as  iso- 
lated blade  rows  in  the  annular  cascade  have  been  presented  in  Reference  1 . 

The  mean  stagger  angles  (at  mid-annulus)  were  = 40°,  as  = 37.2°  or  28.2°, 

and  the  semi-chords  CR  = .0604  ft.  and  Cs  = .054  ft.  The  axial  separation 
of  the  lid-chord  planes  was  b = .125  ft. 

The  sound  pressure  levels  produced  by  rotor-stator  interaction  were 
measured  on  the  outer  casing  upstream  of  the  rotor.  In  taking  these  data,  the 
fan  system  downstream  of  the  annular  cascade  was  turned  off  and  the  fan  was 
allowed  to  rotate  freely  under  the  influence  of  the  flow  generated  by  the 
rotor  in  the  annular  cascade.  In  addition,  all  dampers  in  the  downstream 
drive  system  were  opened  wide.  It  was  convenient  to  use  this  wide  open  con- 
figuration because  it  was  found  to  provide  nearly  constant  inlet  swirl  angles 
relative  to  the  rotor.  That  is,  the. mean  axial  velocity  in  the  annular 
cascade  increased  in  direct  proportion  to  the  rotor  angular  velocity.  This 
feature  made  application  of  the  rotor-stator  interaction  theory  much  simpler 
because  the  dimensionless  steady  state  loss  and  turning  performance  data  for 
the  rotor  also  remained  approximately  constant  with  changing  rotor  speeds. 


A Bruel  and  Kjaer  (B  6 K)  acoustic-probe  with  its  tip  flush  with  the 
inner  surface  of  the  compressor  casing  was  used  in  combination  with  a B § K 
condenser  microphone  and  sound  level  meter  to  measure  the  interaction  noise. 
The  sound  level  meter  output  was  recorded  on  an  instrumentation  quality  FM 
tape  recorder.  The  recorded  data  were  processed  later  on  a modified  General 
Radio  Type  1921  real  time  analyzer  to  obtain  one-third  octave  spectra  of  the 
signals . 


The  tape  recording  of  the  noise  was  necessary  because  the  rotor-stator 
interaction  harmonic  frequencies  were  above  the  frequency  range  available  on 
the  real  time  analyzer.  The  maximum  frequency  accommodated  by  the  spectrum 
analyzer  is  the  third  octave  centered  on  2500  Hz  while  the  noise  harmonics 
extended  up  to  approximately  5000  Hz.  The  effective  frequency  range  of  the 
noise  spectrum  was  reduced  by  a factor  of  four  by  recording  at  30  ips  and  playing 
back  for  spectrum  analysis  at  7-1/2  ips.  Absolute  calibration  of  the  complete 
system  including  the  tape  recorder  was  obtained  by  recording  a 114  db  signal 
from  a General  Radio  microphone  calibrator  at  the  beginning  and  end  of  each 
tape.  These  calibration  signals  were  used  to  set  the  gain  of  the  overall 
system  when  played  back  through  the  spectrum  analyzer. 

Within  the  frequency  range  of  interest,  the  response  of  the  microphone- 
sound  level  meter  system  was  essentially  flat.  However,  the  acoustic  probe 
assembly  introduced  frequency  dependent  distortion  to  this  flat  response. 

Thus,  it  was  necessary  to  obtain  a calibration  for  the  complete  system  in- 
cluding the  acoustic  probe.  A sketch  of  the  acoustic  probe  calibration 
apparatus  is  shown  in  Figure  2.  All  of  the  apparatus  shown  in  the  sketch  is 
manufactured  by  B § K.  The  calibrator  consists  of  a small  cavity  into  which 
sound  is  introduced  by  an  earphone  speaker.  For  a given  excitation  of  the 
speaker,  the  sound  pressure  level  in  the  cavity  is  measured  by  the  microphone 
in  its  reference  level  location  (position  #1) . The  acoustic  probe  is  in 
place  during  these  measurements  with  a dummy  microphone  inserted  in  its  end. 
Following  the  reference  level  measurements,  the  microphone  is  inserted  in  the 
acoustic  probe  (position  #2)  and  the  dummy  microphone  inserted  in  the  reference 
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level  location  (position  #1).  The  sound  pressure  levels  detected  by  the  micro- 
phone in  the  probe  are  then  determined  for  the  same  speaker  excitation  as  used 
in  the  reference  level  measurements.  The  difference  between  the  two  sets  of 
measurements  provides  a frequency  dependent  calibration  for  the  probe. 

In  the  current  work,  the  probe  was  calibrated  first  without  any  steel 
wool  damping  material  in  the  tube.  This  calibration  showed  a rather  large 
(12  db)  probe  resonant  peak  at  approximately  1900  He.  The  probe  was  then 
calibrated  with  three  different  degrees  of  steel  wool  damping  and  a final 
damping  configurat ion  was  selected  for  use  in  the  noise  measurements.  The 
sound  pressure  levels  measured  with  pure  tone  excitation  of  the  earphone 
speaker  are  shown  in  Figure  5 for  the  frequency  range  of  interest  and  for  both 
microphone  positions  in  the  calibrator.  Both  curves  drop  off  at  high  frequency 
because  of  frequency  response  limitations  on  the  earphone  speaker.  The  dif- 
ference between  the  two  curves  provides  the  calibration  of  the  acoustic  probe. 

The  probe  calibration  is  shown  in  Figure  4.  The  circular  points  were 
obtained  with  pure  tone  excitation  of  the  speaker  at  the  indicated  frequencies. 
These  values  must  be  subtracted  from  the  measured  noise  levels.  A second 
calibration  result  is  also  shown  on  this  figure  as  solid  straight  line  segments. 
This  wa.i  obtained  by  exciting  the  speaker  with  pink  noise  and  analyzing  the 
result  on  the  1/3  octavo  real  tin.c  analyzer.  The  effective  1/3  octave  band- 
widths  of  the  results  are  represented  by  the  length  of  the  horizontal  portions 
of  the  segmented  curve.  If  allowance  is  made  for  the  larger  bandwidth  of  the 
1/3  octave  analyzer,  the  agreement  between  the  two  calibration  curves  is  very 
good.  Since  the  rotor-stator  interference  noise  was  very  close  to  a pure  tone, 
the  pure  tone  calibration  curve  was  used  in  correcting  the  measured  results. 

The  pure  tone  nature  of  the  interference  noise  is  illustrated  in 
Figure  5.  This  figure  shows  photographic  oscilloscope  records  of  the  unfiltered 
output  from  the  sound  level  meter.  In  the  upper  photograph,  the  top  trace  is 
generated  by  the  noise  signal  and  the  lower  trace  is  the  output  from  an  audio 
oscillator  tuned  to  the  same  frequency  as  the  predominant  noise  harmonic. 
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The  jitter  in  the  upper  trace  signals  is  caused  by  the  presence  of  frequencies 
other  than  the  predominant  harmonic.  Tuning  of  the  audio  oscillator  was 
accomplished  by  setting  the  oscilloscope  trigger  mechanism  to  fire  on  the 
audio  oscillator  signal  and  then  adjusting  the  oscillator  frequency  until  a 
stationary  noise  signal  was  attained.  This  procedure  led  to  a number  of 
possible  oscillator  frequencies  for  the  noise  harmonic.  The  correct  frequency 
was  difficult  to  determine  because  of  the  jitter  in  the  noise  signal.  The 
oscillator  frequency  was  then  determined  by  forming  a Lissajous  figure  from 
the  noise  and  oscillator  signals.  The  lower  photograph  in  Figure  5 shows  the 
Lissajous  figure  corresponding  to  the  upper  photograph,  The  multiple  ellipses 
in  the  lower  photograph  show  that  the  noise  and  oscillator  signals  are  cor- 
rectly matched  in  fundamental  frequency.  The  multiplicity  of  ellipses  arises 
from  the  same  reason  as  the  jitter  in  the  upper  photograph. 

The  procedure  described  in  the  previous  paragraph  was  used  to  determine 
the  frequency  of  the  predominant  pure  tone  noise  for  all  tests  where  a strong 
blade  row  interaction  signal  was  obtained.  Thus,  it  was  possible  to  identify 
accurately  which  harmonic  of  rotor  blade  passage  frequency  was  being  excited 
by  the  interaction  mechanism. 

Samples  of  the  one-third  octave  spectra  of  the  interaction  noise  are 
shown  in  Figures  6a  through  6d.  These  figures  are  direct  reproductions  of  the 
output  from  the  spectrum  analyzer.  Absolute  sound  pressure  level  for  any 
one-third  octave  band  is  obtained  by  correcting  the  relative  db  levels  for 
system  gain  and  for  acoustic  probe  response.  The  frequency  dependent  probe 
response  correction  has  been  discussed  previously.  The  system  gain  correction 
is  independent  of  frequency  and  is  noted  on  the  lower  right  hand  side  of  each 
figure.  For  the  particular  examples  of  Figure  6,  system  gain  is  corrected  for 
by  adding  80  db  to  all  relative  levels. 

Figures  6a  through  6d  have  been  chosen  to  illustrate  the  extremely 
sharp  rise  in  interaction  noise  as  rotor  speed  is  increased  over  a very  small 
range  of  rpm  from  below  cutoff  to  above  cutoff  of  the  fourth  rotor  blade 
passage  harmonic.  Figures  6a  and  6b  correspond  to  conditions  slightly  below 
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setoff.  Here  the  noise  spectra  are  generated  primarily  by  turbulent  pressure 
fluctuations  within  the  boundary  layer  on  the  outer  wall  and  by  the  free- 
wheeling fan  downstream  of  the  annular  cascade  test  section,  figures  6c  and 

6d  show  the  sharp  rise  in  fourth  harmonic  interaction  noise  as  the  rotor  rpm 
increased  to  above  cutoff  conditions.  This  rise  can  be  seen  in  the  third 
octave  bands  centered  at  Hz  x 1/4  = 800  and  1000.  Both  of  these  third  octave 
bands  respond  because  the  pure  tone  frequency'  is  nearly  mid-way  between  these 
two  bands.  The  measured  pure  tone  frequency  is  shown  on  the  lower  right  hand 
side  of  Figures  6c  and  6d . It  can  be  seen  that  the  pure  tone  frequency  is 
approximately  3600  Hz  for  both  figures  corresponding  to  Hz  x 1/4  = 900. 

The  correct  sound  pressure  level  for  the  pure  tone  signal  was  obtained 
from  spectra  such  as  shown  in  Figure  6 by  combining  the  two  third  octave  band 
responses  closest  to  the  pure  tone  frequency  and  then  correcting  the  indicated 
SPL  for  acoustic  probe  response  at  the  measured  pure  tone  frequency.  An 
example  of  the  procedure  for  the  data  of  Figure  6d  is  as  follows. 

Relative  Level  for  third  octave  centered  at  Hz  x 1/4  = 800  is  34.75  db 

Relative  Level  for  third  octave  centered  at  Hz  x 1/4  = 1000  is  39,25  db 

Combined  Relative  Level  - 20  log  £log  1 + log’1  = 43.3  db 

Combined  Level  Corrected  for  system  gain  = 43.3  + 80  = 123.3  db 
Probe  response  at  pure  tone  frequency  (3610  Hz)  * -4.7  db  (from  Figure  4) 
Corrected  pure  tone  Sound  Pressure  Level  = 123.3  - (-4.7)  = 128.0  db. 

The  pure  tone  frequency  in  all  cases  where  it  was  measurable  corres- 
ponded very  closely  to  either  the  calculated  fourth  or  fifth  harmonic  of  rotor 
blade  passage  frequency.  For  the  example  given  above,  the  calculated  fourth 
harmonic  of  rotor  blade  passage  frequency  is  3588  Hz  and  the  measured  pure  tone 
frequency  is  3610  Hz.  The  difference  between  the  two  is  within  the  accuracy 
of  the  audio  oscillator  used  in  the  measurements.  For  conditions  below  cutoff 
such  as  shown  in  Figures  &a  and  0b,  the  sound  pressure  levels  for  the  fourth 
blade  passage  harmonic  were  estimated  from  the  data  using  the  calculated  fourth 
harmonic  frequency  and  the  measured  levels  in  the  two  third  octave  bands  closest 
to  this  calculated  frequency. 
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The  results  of  the  interaction  noise  measurements  are  summarized  in 
Figures  7 and  8 for  rotor  speeds  ranging  from  below  cutoff  to  the  maximum 
speed  available.  Figure  7 shows  the  results  obtained  for  a stator  stagger 
angle  of  37,2  degrees  and  Figure  8 show.,  the  results  obtaind  for  28.2  degrees. 
In  both  figures,  an  estimate  of  the  background  level  in  the  third  octave  band 
closest  to  the  measured  pure  tone  frequency  is  also  shown.  This  was  obtained 
by  inspection  of  the  third  octave  bands  on  each  side  of  the  bands  used  to 
calculate  the  pure  tone  levels  and  must  be  considered  highly  approximate  in 
nature.  Any  data  point  for  harmonic  level  of  the  SPL  within  about  10  db  of 
the  approximate  background  level  probably  contains  a certain  degree  of  back- 
ground noise  from  the  mechanical  systems  or  the  turbulent  wall  pressure 
fluctuations.  The  closer  the  data  point  is  to  the  approximate  background 
level,  the  larger  the  contribution  will  be, 


Most  of  the  data  shown  in  Figures  7 and  8 correspond  to  the  fourth 
harmonic  of  rotor  blade  passage  frequency.  The  fifth  harmonic  was  generally 
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a stator  stagger  angle  of  28,2  degrees,  the  fifth  harmonic  was  detectable  at 
moderately  low  rotor  speeds,  where  the  background  noise  level  was  low  but  the 
rotor  speed  was  high  enough  to  allow  propagation  of  this  mode  in  the  annular 
duct  system,  The  results  are  shown  on  the  left  side  of  Figure  8, 


It  was  not  possible  to  obtain  interaction  noise  signals  at  harmonics 
of  blade  passage  frequencies  below  the  fourth  harmonic.  The  generation  of 
lower  harmonics  would  require  either  a greater  rotor  speed  capability  or  else 
blade  configurations  that  are  not  available  to  the  annular  cascade  in  its 
present  configuration.  It  is  unfortunate  that  the  generation  of  lower  har- 
monics was  not  possible  because  the  prediction  of  harmonics  as  high  as  the 
fourth  presents  a test  for  the  interaction  theory  which  may  be  more  severe 
than  is  warranted  in  the  present  state  of  the  theoretical  development.  Never- 
theless, these  data  have  been  used  for  comparison  with  the  theoretical  pre- 
dictions in  Section  III,  and  are  discussed  more  fully  there. 
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ISOLATED  ROTOR  MEASUREMENTS 
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In  addition  to  the  rotor-stator  interaction  noise  experiments  described 
above,  measurements  were  made  of  the  fluctuating  wall  pressures  produced  by  an 
isolated  rotor.  These  measurements,  which  were  intended  to  provide  information 
on  the  blade  tip  loading,  were  taken  at  six  ehordwise  locations  along  the 
outer  casing  wall , 

The  configuration  of  the  annular  cascade  was  similar  to  that  shown  in 
Figure  1 but  with  the  stator  row  removed.  Rotor  Set  No.  1 was  used  with  the 
blade  stagger  angle  set  at  40  degrees  at  mid-annulus  (48  degrees  at  the  tip). 

As  noted  earlier,  this  rotor  set  had  been  studied  previously  to  determine  its 
steady  state  turning  and  loss  performance  (Reference  1) . 

Pressure  fluctuations  produced  by  the  passage  of  the  blades  weze 
obtained  at  rotor  speeds  of  600,  700,  800,  900  and  10QO  rpm,  while  the  mean 
axial  velocity  through  the  test  section  was  held  at  60  ft/sec.  At  the  lowest 
rotor  speed,  the  rotor  blades  were  very  lightly  loaded,  while  at  the  highest 
speed  the  loading  was  large  enough  to  induce  rotating  stall. 

The  pressure  signals  were  measured  by  a Setra  Systems  Model  242  TC 
pressure  transducer  closely  coupled  to  a hole  in  the  outer  casing  of  the  annular 
cascade.  The  time-varying  pressure  signals  were  recorded  photographically  from 
a dual  trace  oscilloscope.  A second  signal  was  used  to  indicate  rotor  blade 
position.  This  signal  was  generated  by  a magnetic  pickup  which  sensed  blade 
tip  passage  past  a point  or.  the  outer  casing  of  the  annular  cascade.  Instan- 
taneous blade  positions  with  respect  to  the  pressure  signals  were  calculated 
from  the  magnetic  pickup  signals  and  the  geometric  relationship  between  the 
magnetic  pickup  and  the  pressure  tap. 

The  six  avial  stations  used  for  the  wall  pressure  measurements  were 
situated  at  15,  18,  32,  49,  64  and  79  percent  of  the  rotor  chord.  A set  of 
results  i.s  shown  in  Figure  9 for  the  station  at  18  percent  of  the  chord. 
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These  results  are  typical  of  the  measurements  at  all  six  chordwise  measuring 
stations.  Parts  (a)  through  (d)  of  Figure  9 show  the  records  obtained  at 
rotor  speeds  between  000  and  900  rpm.  The  relative  inlet  swirl  angle  at  the 
rotor  tip  is  indicated  for  each  rotor  speed.  The  instantaneous  locations  of 
the  blade  pressure  and  suction  surfaces  with  respect  to  the  pressure  signals 
are  shown  also. 

As  expected,  the  maximum  amplitude  of  the  pressure  fluctuations  in- 
creases with  increasing  rotor  speed.  (Note  that  the  scale  for  pressure 
signals  varies  fGr  different  rotor  speeds.  The  scale  in  each  case  is  indi- 
cated to  the  left  of  the  pressure  signals.)  However,  the  shapes  of  the 
fluctuating  pressure  signals  are  unexpected.  The  data  records  show  that  the 
maximum  wall  pressure  is  reached  well  ahead  of  the  blade  pressure  surface  at 
all  rotor  speeds.  Moreover,  a double  pressure  peak  occurs  at  the  lower  rotor 
speeds,  one  ahead  of  the  rotor  pressure  surface  and  one  approximately  at  the 
location  uf  the  pressure  surface.  The  peak  in  pressure  at  the  rotor  pressure 
surface  begins  to  d isnpDC2.r  3. s rotor  speed  is  increased.  The  ressor.  for  the 
unexpected  shape  of  the  pressure  signals  is  not  definitely  known,  However, 
it  is  speculated  that  the  results  may  be  caused  by  three-dimensional  effects 
associated  with  the  relatively  large  clearances  between  the  rotor  blades  and 
the  outer  casing.  As  explained  in  Volume  I,  a blade  tip  clearance  of  approxi- 
mately 0.05  inches  was  required  to  compensate  for  the  slightly  oval  shape  of 
the  production  J-79  casing  which  forms  the  outer  wall  of  the  annular  cascade. 

The  fluctuating  wall  pressures  recorded  during  rotating  stall  are 
shown  in  Figure  10  for  ai.1  six  chordwise  locations.  Although  the  rotating 
stall  phenomenon  is  beyond  the  scope  of  the  direct:  lifting  surface  theory, 
these  data  are  presented  as  being  of  general  interest  to  the  overall  program; 
pressure  signals  such  as  the  ones  shown  in  Figure  10  are  used  as  stall  de- 
tectors in  the  rotating  stall  control  system  which  is  described  in  Volume  III 
of  this  report. 
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The  rotor  speed  for  the  records  in  Figure  10  was  1000  rpm  and  the  mean 
axial  velocity  was  60  ft/scc.  The  oscilloscope  sweep  speed  in  Figure  10  is 
much  slower  than  in  Figure  9 so  that  the  blade  passage  pressure  peaks  are 
highly  compressed.  The  interest  here  is  centered  on  the  pressure  fluctuations 
during  passage  of  a rotating  stall  cell.  Each  photograph  in  Figure  10  shows 
the  passage  of  two  stall  cells.  The  phase  relationships  between  the  separate 
photographs  have  no  meaning  since  the  records  were  not  obtained  simultaneously. 

Figure  10  illustrates  that  the  character  of  the  stall  cell  passage 
signals  change  considerably  with  ehordwise  location  on  the  rotor.  Near  the 
leading  edge,  the  stall  passage  signals  coincide  with  an  amplitude  reduction 
in  the  blade  passage  signals.  In  this  region  the  combined  maximum  amplitude 
of  the  pressure  signals  during  stall  passage  is  almost  the  same  as  the  blade 
passage  pressure  amplitude  between  stall  cells.  In  contrast,  near  the  trailing 
edge  the  passage  of  a stall  cell  coincides  with  an  increase  in  amplitude  of 
the  blade  passage  signals  and  the  combined  maximum  amplitude  is  larger  than 
the  blade  passage  amplitude  between  stall  cells, 

The  rotating  stall  control  system  detects  unusually  large  peak 
amplitudes  in  pressure  signals  such  as  shown  in  Figure  10.  Control  action  is 
taken  when  these  peak  amplitudes  reach  a predetermined  reference  value.  It 
is  required  that  only  those  fluctuating  pressure  levels  associated  with  rota- 
ting stall  should  initial  control  action.  Tests  of  the  control  have  shown 
that  the  best  performance  is  obtained  when  pressure  signals,  due  to  blade 
passage,  are  removed  by  low-pass  filters.  It  is  clear  from  Figure  10  that 
elimination  of  the  blade  passage  signals  will  greatly  improve  the  signal  to 
noise  ratio  for  stall  detection  purposes.  This  is  true  even  for  the  pressure 
tap  locations  near  the  blade  trailing  edge  since  the  blade  passage  pressure 
fluctuations  still  have  an  appreciable  amplitude  between  the  stall  cells, 

D.  CONCLUDING  REMARKS 

The  results  of  two  separate  experimental  investigations  have  been 
presented  in  Section  II.  In  the  first  investigation.  Section  11-3,  the 
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far -l' i.o Id  noise  generated  by  rotor -stator  interaction  was  measured  to  provide 
data  for  correlation  with  the  approximate  theory  presented  in  Section  111. 

In  the  second  study,  Section  II-C,  the  fluctuating  outer-wall  static  pressures 
generated  by  passage  of  the  blade  tips  of  an  isolated  rotor  were  measured  to 
provide  data  for  comparison  ..  Lth  the  direct  lifting  surface  theory  of  Section  IV. 

The  rotor-stator  interference  noise  studies  presented  in  Section  1 1 - B 
consisted  of  measuring  the  sound  pressure  levels  in  the  constant  area  annulus 
far  upstream  of  a rotor-stator  stage.  The  measurements  were  made  at  a series 
of  rotor  speeds  for  two  stagger  angle  settings  of  the  stators,  The  results 
contained  detectable  pure  tone  components  only  at  the  fourth  and  fifth  har- 
monies of  blade  passage  frequency.  The  generation  of  lower  harmonies  was  not 
possible  because  the  rotor  speed  capability  was  limited  at  the  time  of  the 
tests.  Since  the  performance  of  these  tests,  the  rotor  speed  capability  of 
the  test  rig  has  been  increases.  It  is  planned  to  use  this  capability  to 
perform  additional  measurements  for  conditions  in  which  lower  harmonics  of 
blade  passage  frequency  are  propagating.  The  results  of  the  current  tests 
are  compared  with  theoretical  predictions  in  Section  III. 

The  isolated  rotor  studies  of  Section  II-C  were  intended  to  provide 
a measure  of  the  blade  tip  pressure  loading  for  comparison  with  the  predictions 
of  the  direct  lifting  surface  theory  of  Section  IV.  However,  the  theoretical 
development  had  not  reached  the  stage  where  quantitative  predictions  could 
be  made.  Thus,  samples  of  some  of  the  data  are  presented  as  of  interest  in 
themselves.  Inspection  of  the  results  suggest  the  measured  pressures  may 
have  been  influenced  by  rotor  tip  clearance  effects,  As  such,  the  results 
may  not  be  predictable  by  the  direct  lifting  surface  theory.  Nevertheless,  a 
comparison  between  theory  and  experiment  would  still  be  of  interest. 


SECTION  III 

APPROXIMATE  MODEL  OP  ROTOR-STATOR  INTERACTION  NOISE 


The  interactions  between  rot.or  and  stator  blade  rows  have  long  been 
recognized  as  a major  source  of  noise  in  subsonic  axial  flow  fans  and  com- 
pressors. The  interactions  consist  of  fluctuating  forces  which  arise  because 
of  the  motion  of  the  blade  rows  relative  to  one  another,  and  in  turn,  act  as 
acoustic  sources.  Out  of  practical  necessity,  little  detailed  attention  can  be 
paid  to  the  interaction  noise  at  the  design  stage,  where  each,  row  is  usually 
modeled  as  3n  isolated  two-dimensional  cascade  in  a steady  undisturbed  flow. 
Interactions  can  be  kept  to  a minimum  by  spacing  the  rows  several  chord 
lengths  apart,  but  the  designer  is  generally  working  under  size  and  weight 
constraints  as  well.  This  portion  of  the  report  presents  the  results  to  date 
of  a combined  theoretical  and  experimental  program  aimed  at  a better  under- 
standing of  rotor-stator  noise  generation,  and  methods  for  its  alleviation 

Many  investigations  have  appeared  in  the  literature  which  treat,  both 
the  aerodynamic  and/or  acoustic  aspects  of  the  problem,  e.g.,  References  2-26, 
which  is  by  no  means  an  all-inclusive  list.  Because  of  the  complexity  of  the 
problem,  the  various  theoretical  models  represent  several  combinations  of 
simplifying  assumptions  needed  to  make  them  amenable  to  analysis.  These  include 
the  use  of  frec-field  vs.  ducted  boundary  conditions,  two-dimensional  vs.  axi- 
symmetric  blade  rows,  and  various  degrees  of  approximation  to  the  unsteady  aero- 
dynamic processes.  Probably  the  most  universal  approximation,  and  the  most 
restrictive,  has  been  the  use  of  incompressible  flow  theory  to  estimate  the 
fluctuating  blade  loads.  Strictly  speaking,  this  restricts  the  range  of 
validity  of  these  models  to  very  low  flow  speeds . 

As  the  speeds  of  modern  turbomachinery  are  definitely  subsonic,  and 
often  transonic,  the  need  to  include  compressibility  effects  is  obvious.  The 
goaL  of  the  present  theoretical  work  is  to  incorporate  a compressible  flow 
aerodynamic  model  into  the  prediction  of  rotor/stator  interaction  noise. 
Published  investigations,  which  allow  for  compressibility  effects  in  the  pre- 
diction of  unsteady  aerodynamic  loads,  include  those  by  Kaj i and  Okazaki, 
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Kaj  i anil  Okazaki  treat  the  near- fie  Id  aerodynamic  and  far- fi  eld 
acoustic  regions  of  the  flow  in  a unified  linearized  tieatinent.  Their  analysis 
is  the  most  complete  because  it  includes  the  upwash  velocities  on  each  blade 
generated  not  only  by  the  forces  on  that  same  blade,  but  also  by  the  forces 
(both  steady  and  unsteady)  on  all  the  other  blades  in  that  same  row,  as  well 
as  those  from  blades  in  the  neighboring  row.  Unfortunately,  this  requires  one 
to  solve  simultaneously  for  the  loading  distribution  on  both  rows  using  a pair 
of  coupled  singular  integral  equations.  Man i ^ simplified  the  problem  some- 
what by  neglecting  the  influence  of  the  unsteady  loading  on  the  neigh. boring 
row.  It  is  still  necessary  to  solve  a pair  of  integral  equations  to  obtain  the 
unsteady  loading  on  both  rows,  but  the  two  arc  no  longer  coupled.  Similarly, 
the  aerodynamic  analyses  in  References  2S  and  26  require  the  numerical  solution 
of  an  integral  equation  for  the  loading  on  each  row.  These  analyses  are  all 
based  on  two-dimensional  cascade  models. 

It  is  ultimately  hoped  that  our  work  on  the  three-dimensional  lifting 
surface  theory  for  annular  blade  rows  can  be  applied  to  the  rotor-stator 
interaction  problem.  This  would  amount  to  the  extension  of  the  Kaj i and 
Okazaki  analysis  to  include  three-dimensional  effects,  and  is  expected  to  lead 
again  to  a pair  of  coupled  integral  equations,  probably  even  more  complex  than 
theirs,  lienee,  it  is  likely  that  more  experience  with  the  problems  of  steady 
arid  unsteady  flow  through  an  isolated  rotor  (reported  on  in  Section  IV,  below) 
would  be  needed  before  applying  Such  an  analysis  to  rotor-stator  interaction. 

In  the  interim,  however,  it  was  felt  that  our  understanding  of  the  basic 
mechanisms  could  benefit  greatly  from  the  application  of  a simpler  model,  and 
the  comparison  of  its  theoretical  predictions  against  the  acoustic  data  taken 
in  the  annular  cascade  (see  Section  II). 

The  theory  developed  under  the  current  program  avoids  tiie  necessity 
of  solving  integral  equations,  and  thus  considerably  reduces  the  complexity  of 

I 

the  analysis,  as  well  a-,  computing  times.  In  this  simplified  model,  the  aero- 
dynamic and  acoustic  problems  are  treated  individually  and  then  matched  together. 

24  “’l 

The  aerodynamic  analysis  is  that  published  by  Osborne,  *"  which  essentially 
represents  the  compressible  extension  of  the  now  classic  analysis  by  Kemp  and 
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Sears.'’’'1  Osborne's  results  are  closed  form  expressions  for  the  unsteady 
blade  loads  on  both  rotor  and  stator  which  are  easily  accommodated  in  the 
acoustics  calculation,  since  no  numerical  techniques  are  required.  The 
principal  features  of  this  model  are  described  in  Section  III-B  below. 

The  Osborne  aerodynamic  analysis  models  each  row  as  a two-dimensional 
cascade,  which  is  strictly  justifiable  only  at  large  hub/tip  ratios;  even  then, 
McCune"  has  shown  that  in  certain  types  of  transonic  flow,  no  logical  two- 
dimensional  cascade  limit  exists.  The  acoustic  analysis  described  in  Section  1 1 1 - A 
below  employs  axisymmetric  annular  blade  rows  housed  in  an  infinite  hard-walled 
cylindrical  duct  with  a uniform  axial  flow.  Accordingly,  Osborne's  expressions 
for  the  fluctuating  loads  are  applied  on  a strip  theory  basis  at  each  radial 
station.  This  procedure  was  followed  so  that  in  the  event  a truly  three- 
dimensional  aerodynamic  analysis  becomes  available,  it  can  more  easily  be  in- 
corporated in  the  model. 

The  combined  aero-acoustic  analysis  in  Section  III-C  results  in  easily 
evaluated  expressions  for  the  sound  field  in  terms  of  the  same  duct  modes 
studied  previously  by  Tyler  and  Sofrin  and  Morfey.'’’^  A computer  program  has 
been  written  to  evaluate  the  modal  amplitudes  in  terms  of  blade  row  steady 
loadings,  stagger  angles  and  drag  coefficients.  From  these  amplitudes,  the 
total  radiated  sound  power  (either  upstream  or  downstream) , as  well  as  the 
mean  square  pressure  at  any  given  point  in  the  duct  may  be  obtained. 

Section  III-D  contains  comparisons  made  between  the  theoretical  pre- 
dictions and  acoustic  data  taken  in  the  annular  cascade  facility.  Additional 
numerical  results  are  presented  in  Section  III-E  to  better  illustrate  the 
effects  of  compressibility.  Section  III-F  summarizes  the  findings  and  makes 
suggestions  as  to  how  the  theoretical  model  can  be  improved. 


A. 


ACOUSTIC  ANALYSIS 


In  the  acoustic  model  studied,  the  blade  rows  are  assumed  to  be 
housed  in  an  infinitely  long  hard-walled  cylindrical  duct,  as  shown  in 
Figure  11.  The  hub  and  tip  radii  are  denoted  by  rH  and  rr  respectively,  and 
a uniform  axial  flow  at  Mach  number  Mx  is  assumed.  In  describing  the  flow 
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field,  we  will  use  Lighthill's  acoustic  analogy  procedure  as  later  extended 
by  Curle-5^  to  allow  for  the  presence  of  solid  boundaries.  That  is,  rather  than 
having  the  blades  impose  boundary  conditions  on  the  equations  of  motion,  we 
represent  them  as  distributed  sources  of  mass  and  momentum  which  arise  from 
the  blades'  thickness  and  loading,  and  imagine  the  fluid  to  be  otherwise  un- 
bounded in  the  annulus  formed  by  the  duct  walls.  The  linearised,  inviscid 
conservation  equations  of  mass,  momentum  and  energy  in  duct-fixed  coordinates 
are  then: 
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and  p , -p  and  y are  the  perturbation  density,  pressure  and  velocity  fields, 
respectively.  The  quantities  ag  and  Ux  represent  the  undisturbed  sound  speed 
and  the  axial  flow  velocity,  respectively.  Q is  the  rate  of  introduction  of 
mass  per  unit  volume  per  unit  time,  and  F is  a force  per  unit  volume;  both 
are  functions  of  x and  t . Equations  (l)-(3)  are  easily  manipulated  into 
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5 ( x , t ) 


The  only  boundary  condition  imposed  on  this  equation  is  that  the  radial 
velocity  vanish  at  the  duct  walls.  Using  Equation  (2),  this  is  expressed  in 
terms  of  -o  as 
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In  order  for  Equations  (4)  and  (5)  to  be  of  any  use,  one  must  have  a 
priori  knowledge  of  the  acoustic  source  distribution,  S(  x ,t)  . This  is  dis- 
cussed in  the  next  section;  for  now  we  simply  observe  that  solutions  of  (4) 
and  (5)  can  be  obtained  using  a number  of  methods.  The  present  analysis  uses 
a Fourier- Bessel  transform  decomposition.  That  is,  we  define  transform  variables 
pn„(Z,u>)  and  Sn^(%,u>)  , corresponding  to  and  S (.x  , r)  , which 

are  expressed  in  the  cylindrical  polar  coordinates  of  Figure  11  as: 
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Here  r ■ r/rT  is  a convenient  dimensionless  radial  variable  and  -A  * r.  / r_ 
is  the  hub/tip  ratio.  The  functions  C vnm<r)  arc  the  set  of  orthonormal 

radial  eigenfunctions  imposed  on  the  system  by  Equation  (5);  they  are  a linear 
cf'^h. nation  of  the  Bessel  and  Neumann  functions,  and  are  described  more  fully 
in  Appendix  A-  The  corresponding  eigenvalues,  , are  determined  by  the 

condi t ion 
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By  substituting  for  -p  and  S in  liquation  (4j  in  terms  of  the  in- 
verse transforms  appropriate  to  (6)  and  ( 7 J , one  can  show  that  liquation  (4)  is 
then  automatically  satisfied  when  the  tran  form  variables  have  the  following 
relationship  to  each  other, 
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Thus,  once  >5  is  known,  .5„^,,  and  follow  from  Equations  (7)  and  (9), 

respectively.  The  solution  for  the  pressure  field  is  then  obtained  as  the 
inverse  transform  of  Equation  (6).  Combining  these  steps  we  obtain, 
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where  everything  on  the  right  is  presumed  known.  Tor  the  present  study  of 
rotor-stator  interaction  noise,  we  are  concerned  with  the  fluctuating  forces, 
and  so  drop  the  source  term  Q appearing  in  Equation  (4). 


B.  AERODYNAMIC  MODEL 


It  remains  then  to  specify  the  form  of  the  unsteady  forces  F exerted 
by  the  blades  on  the  air.  Lor  the  reasons  discussed  above,  the  linearized 
aerodynamic  analysis  of  Osborne^  was  chosen  for  this  purpose.  Briefly  stated, 
the  principal  assumptions  in  this  analysis  are: 

(lj  Each  blade  is  represented  by  a two-dimensional  cascade  in 
evaluating  its  induced  velocity  field. 
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(j)  The  steady  part  of  the  circulation  about  each  blade,  and  hence 
its  influence  on  a neighboring  blade  row,  is  much  greater  than 
the  unsteady  circulation.  The  latter  is  neglected,  along  with 
any  associated  shed  vorticity. 

(3)  Unsteady  velocity  gusts  parallel  to  the  blade  chord  are 

neglected.  This  is  safe  provided  the  steady  state  angle  of 
31 

attack  is  small,  which  is  implicit  in  the  linearization. 

(4 ) Isolated  two-dimensional  airfoil  theory  is  used  to  estimate 
blade  response  to  gusts  perpendicular  to  the  chord. 

As  in  the  original  analysis  bv  Kemp  and  Scars,  Osborne  considers  a 
single  stage  consisting  of  a stator  followed  by  a rotor  (see  Figure  12a). 

Since  we  have  assumed  a uniform  axial  velocity  in  the  duct  as  the  undisturbed 
state  m the  acoustic  model,  the  linearized  analysis  is  strictly  applicable 
only  to  stages  with  small  stator  stagger  angles.  Three  interaction  mechanisms 
are  considered: 

(a)  Rotor  unsteady  lift  fluctuations  due  to  its  passage  by  the 
steady  upstream  stator  loads. 

(b)  Rotor  unsteady  lift  due  to  its  passage  through  the  viscous 
wakes  shed  by  the  upstream  stator. 

(c)  Stator  unsteady  lift  due  to  the  passage  of  the  steady  rotor 
loads . 

For  brevity,  items  (a)  and  (c)  are  usually  referred  to  as  potential  inter- 
actions since  they  would  be  present  oven  if  the  fluid  were  ideal.  Item  (b) , 
referred  to  as  the  viscous  interaction,  requires  the  presence  of  fluid  friction. 
The  potential  interaction  analysis  involves  the  blade  aerodynamic  response  to 
a generalized  Kemp-type  gust  upwash,  whereas  the  viscous  interaction  is  modeled 
as  the  superposition  of  responses  to  a frozen  sinusoidal  gust. 
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The  sectional  lift  per  unit  span  on  the  zeroth  blade  given  by 
Osborne  for  each  of  the  above  mechanisms  may  be  symbolically  expressed  in 
blade-referenced  coordinates  as 
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where  X = and  similarly  for  and  . 


The  subscripts  R and  s refer  to  a blade  member  of  the  rotor  and 
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equation  numbers  refers  to  the.ir  corresponding  interactions  (a,  b,  c)  above. 
The  Xjj  , Yj.  , and  2t  are  the  amplitudes  of  the  fluctuation  at  each  harmonic 
of  the  fundamental  frequency,  whether  Va  Or  . (Mote:  since  we  are  not 

explicitly  concerned  here  with  the  noise  field  generated  by  the  steady  rotor 
load,  the  1 = 0 term  in  [11a,  b]  is  ignored.)  These  coefficients  are  given 
explicit/  by  Osborne  as  functions  of  blade  row  stagger,  spacing  and  operating 
conditions.  The  full  functional  relationships  would  needlessly  confuse  the 
analysis  that  follows,  and  so  they  are  given  in  Appendix  B.  Since  Osborne's 
model  is  a two-dimensional  one,  x^  , Vj  and  are  independent  of  spanwise 
location  in  his  analysis;  in  the  strip-wise  application  of  his  results  to  our 
annular  duct  model,  they  are  assumed  to  be  implicit  functions  of  the  radial 
coordinate  r shown  in  Figure  11. 


The  fundamental  frequencies  and  correspond  to  the  radian 
frequency  with  which  a rotor  blade  (stator  vane)  encounters  a stator  vane 
(rotor  blade).  In  the  cascade  model,  this  is  related  to  the  rotor  tangential 
velocity  divided  by  the  stator  vane  (rotor  blade)  spacing.  In  the  present 
circular  duct  model,  this  translates  into  the  product  of  the  stator  vane  (rotor 
blade)  number  times  the  angular  velocity  of  rotation,  i.e.. 
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where  U and  B are  the  numbers  of  vanes  and  blades,  respectively,  ^ and 

J a their  spacing,  and  XI  is  the  shaft  frequency  in  radians/sec. 

C.  COMBINED  AERO-ACOUSTIC  ANALYSIS 


5 a 


(12) 


Equation  (4)  shows  that  the  contribution  to  the  source  distribution 
S (x  ,r)  from  these  fluctuating  forces  per  unit  span  is  in  the  form  of 
minus  the  divergence  of  the  fluctuating  forces  per  unit  volume,  expressed  in 
duct-fixed  coordinates.  That  is,  the  blade  rows  arc  imagined  replaced  by  a 
fictitious  distribution  of  acoustic  dipole  singularities,  stationary  in  duct 
coordinates,  whose  magnitude  and  phase  are  in  accordance  with  Equations  (lla-c). 
To  determine  this  distribution,  some  assumption  must  be  made  regarding  the 
spatial  distribution  of  blade  forces  in  the  axial  (x)  and  azimuthal  (6)  co- 
ordinates. For  simplicity,  we  have  initially  assumed  each  blade  (vane)  to  be 
represented  by  a line  dipole  located  at  its  mid-chord.  (More  sophisticated 
models  employing  a chordwise  distribution  of  loading  can  be  treated,  as  this 
assumption  is  not  crucial  to  the  acoustic  analysis.)  Noting  that  Osborne  takes 
the  rotor  as  having  velocity  U (*xir)  in  what  we  define  as  the  negative  6 
direction  (Figure  12a),  the  fluctuating  force  field  corresponding  to  Equations 
(lia-c)  may  then  be  expressed  in  duct-fixcd  coordinates  as 
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where  X,  , x„  are  the  axial  locations  of  the  stator  and  rotor  dipoles.  The 
dipole  axes  are  perpendicular  to  the  rotor  blade  (stator  vane)  chord  lines. 
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which,  in  turn,  are  inclined  at  an  angle  as'-<xs)  from  the  axial  direction, 
as  in  Figure  Ida.  Hence,  the  divergence  of  the  forces  occurring  at  tne  rotor 
can  be  shown  to  be  (neglecting  any  contribution  from  the  radial  direction): 
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Substitution  of  Equations  (13a-c)  into  the  above  expressions  gives  the 
following  contributions  to  S from  each  mechanism: 
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Substituting  Equations  (16a-c)  in  (3)  then  gives  the  expansion  coefficients 
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where  the  sum  over  the  integer  N has  replaced  that  over  .i  . The  notation  £ 
indicates  that  Equation  (1?)  should  be  viewed  as  n being  a fixed  integer  and 
allowing  only  those  teams  in  the  sum  over  N which  give  integral  values  of  X . 
This  restriction  arises  because  of  the  following  identity  which  is  used  in  the 
derivation  of  Equations  (,17a)  and  (17b). 
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Similarly,  in  the  case  of  Equation  (17c),  use  was  made  of 
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Physically,  these  identities  reflect  the  phase  cancellations  in  the  duct  which 
allow  only  a restricted  set  of  Loading  harmonics  to  generate  a given  azimuthal 
mode. 
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1. 


Press u r c_ Mod c Amplitudes 

Substitution  of  liquation  (II7)  into  (lOj  gives  the  desired  solution 
for  the  acoustic  pressure  field.  The  integral  over  uj  is  rendered  trivial  by 
the  presence  of  the  delta  functions  in  Equation  (17).  The  integral  over 
can  be  handled  in  a straightforward  manner  using  contour  integration  in  the 
complex  £,  plane.  The  integrand  exhibits  two  simple  poles  on  the  real  axis. 
The  physical  requirement  that  all  waves  generated  by  a row  should  propagate 
away  from  it  determines  which  half-plane  the  poles  should  bo  considered  to 
lie  in,  ai.d  the  value  of  the  integral  can  then  be  found  from  residue  theory. 


For  the  sake  of  generality,  pressure  waves  generated  in  the  upstream 
direction  by  cither  row  are  referred  to  with  a superscript  "u”,  and  those  in 
the  downstream  direction  with  a "d".  A subscript  of  either  1 or  2 will  be 
used  to  denote  whether  they  originated  at  the  upstream  or  downstream  blade  row 
(see  Figure  13).  (In  tnis  simplified  model,  no  allowance  is  made  for  the  re- 
flection or  transmission  of  waves  by  a neighboring  row  once  they  arc  generated.) 
This  formalism  allows  us  to  present  results  for  either  the  stator-rotor  case 
discussed  thus  far,  or  the  rotor-stator  case  to  be  discussed  shortly,  within 
the  same  context.  Then,  each  harmonic  of  the  pressure  field  is  found  to  con- 
sist of  four  sets  of  waves  as  given  below; 
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i?  tiic  particular  harmonic  of  blade  passage  frequency  being  considered. 

The  values  included  in  the  sum  over  the  azimuthal  and  radial  mode  numbers, 
n and  no  , will  be  discussed  shortly,  liquations  (ISa-d)  hold  regardless  of 
the  relative  positions  of  the  rotor  and  stator.  This  need  not  be  specified 
until  the  mode  amplitudes  ~p  , which  are  in  general  complex,  arc  calculated. 

For  the  stator- rotor  case  discussed  above,  these  are: 
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The  upper  signs  apply  upstream,  and  the  lower  signs  downstream. 

Iri  many  applications,  particularly  high,  bypass  ratio  fans,  the  per- 
tinent geometry  is  that  of  a rotor  upstream  of  a stator  (outlet  guide  vanes; . 
The  same  three  interactions  (a,  b,  c)  listed  at  the  beginning  of  Section  IJ1-B 
are  stili  involved,  provided  one  interchanges  the  words  rotor  and  stator. 
Rather  than  having  to  redo  Osborne's  entire  analysis,  the  appropriate  forces 
can  be  obtained  from  his  present  expressions  using  a simple  transformation 
discussed  in  Appendix  B,  and  the  results  written  down  directly: 
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Equations  (18a-d)  express  the  pressure  field  as  the  superposition  of 
the  same  spiral  duct  modes  studied  by  Tyler  and  Sofrin,4  and  Morfev.S’f>  Note 
that  the  rotor-generated  field  is  made  up  of  modes  with  azimuthal  order 
tif  B - X / , where  N(6>fl  is  the  acoustic  frequency  and  J.VCI  the  blade  loading 
frequency.  The  stator-generated  field  on  the  other  hand  is  composed  of  modes 
with  azimuthal  order  J 3 - NV  , where  both  the  acoustic  and  vane  loading  fre- 
quencies - JBn  • NfBn.  These  rules  for  determining  the  azimuthal  modes 
over  which  one  must  sum  in  Equations  (18a-d)  are  a result  of  the  delta  functions 
involving  uj  in  Equation  (17j  in  conjunction  with  the  special  relationships 
between  -a  , N and  X given  there.  Interestingly,  one  can  see  that  the  rotor 
and  stator  always  will  each,  excite  the  same  set  of  modes.  These  obser- 

vations hold  irrespective  of  which  is  the  upstream  row. 
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Here  we  are  concerned  only  with  the  radiated  field,  and  so  only  those 
modes  above  cut-off,  ie.,  for  which  is  real,  need  be  considered.  This 

allows  one  to  pat  an  upper  limit  on  the  value  of  the  radial  index  m , for  a 
given  frequency  u>  and  azimuthal  mode  number  n . (The  form  of  the  result  for 
modes  which  arc  cut  off  is  t.hc  same  as  (IS),  except  that  -k  is  then 
imaginary;  its  sign  is  determined  from  the  physical  requirement  that  the  mode 
decay,  rather  than  grow,  exponentially.) 

The  mode  amplitudes  given  by  equations  (19a-d)  are  singular  right  at 
the  cutoff  condition  for  any  one  mode,  -k  « 6 y.nrn  , as  first  reported  by 

•7Jj 

McCune”  in  his  study  of  the  disturbances  generated  by  the  steady  loading  or. 
an  isolated  rotor.  In  References  22,  25  and  25,  which  treat  the  aerodynamic 
and  acoustic  fields  simultaneously  in  a two-dimensional  cascade  model,  no  such 
resonance  is  predicted  since  the  unsteady  loads  (which  we  have  assumed  as  given 
a priori;  apparently  vanish  ar  this  condition,  due  to  a strong  interaction 
between  the  two  fields.  The  experimental  results  discussed  below  and  in 
Section  II,  however,  indicate  that  the  sound  levels  are  significantly  higher 
near  such  a condition. 


Sound  Pressure  Level  and  Total  Radiated  Power 


Once  tile  pressure  mode  amplitudes,  -p  , have  been  determined  from 
Lquation  (19),  expressions  for  both  the  mean  square  pressure  at  any  point  in 
the  duct  and  the  total  radiated  power  are  easily  derived  for  any  given  frequency. 
We  denote  time  averages  over  one  blade  passage  interval  by  < > . Then  since 
everything  is  harmonic  in  time, 
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and  i <1 kja  i is  the  magnitude  of  the  sum  of  the  .signals  from  both  rows,  i.e.,  the 
coot  fioient  of  c "A't  in  liquation  [21),  after  substituting  from  [18)  and  (!9J 
as  appropriate.  The  result  xs  easily  shown  to  he: 
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The  last  equality  follows  from  the  fact  that  the  two  blade  rows  always  excite 
the  same  set  of  {n.rn}  modes,  (but  with  different  ampl itudes,  of  course). 

Note  that  <-p1>  U"’L  is  a function  of  x,  *"  - r,  <r  , and  8 , but  is  independent  of 
time;  the  dependence  on  x and  B will  also  disappear  in  situations  where  only 
one  mode  is  excited.  The  sound  pressure  level,  SPL,  then  follows  from 
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To  obtain  the  total  sound  power  radiated  along  the  duct,  P1*^  , we 

mist  evaluate  (Rc-ieronce  0) 
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where  u is  the  axial  component  of  the  acoustic  perturbation  velocity.  The 
last  tern:  is  proportional  to  the  integral  of  Liquation  (22)  over  the  duct  cross- 
section;  it  simplifies  considerably  due  to  the  orthogonality  properties  of  the 
trigonometric  and  Bessel  functions  and  the  fact  that,  for  two  complex 
quantities , 
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The  integral  of  then  reduces  to 
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The  other  terms  in  Equation  (24)  involve  the  axial  velocity  in  some 
form.  For  every  pressure  wave  with  amplitude  there  is  a corresponding 

velocity  wave  with  the  same  spatial  and  time  dependence  and  an  amplitude  of, 
say,  o.  ^ * . The  two  amplitudes  can  be  related  through  the  axial  component  of 
the  momentum  equation: 
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The  quantity  may  be  thought  of  as  the  dimensionless  acoustic  admittance 

of  an  mode.  For  the  plane  wave  mode  with  no  axial  flow,  n r rv>  «■  M K - o , 

it  reduces  to  +1,  as  expected.  Since  is  a real,  constant  multiplicative 

factor  independent  of  x,  r and  8 , the  contributions  to  P from  the  first 
two  terms  in  Equation  (24)  can  be  written  down  directly  from  (25) . 

The  complete  result  can  then  be  reduced  to: 
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which  is  written  in  such  a way  that  P1*  and  P*1  are  both  positive,  i.e.,  power  is 
always  radiated  away  from  the  stage. 
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The  terms  involving  the  squares  of  the  mode  amplitudes  represent  the  power  put 
out  independently  by  each  blade  row,  assuming  only  its  radiation  was  pic-sent. 
The  last  term  represents  an  interference  between  waves  emanating  from  the 
separate  rows;  this  would  be  zero  if  the  two  blade  rows  radiated  independently 
of  one  another,  i.e,  if  (arg  -p  t,™,  ~ a.rg ) were  random. 

If  the  acoustic  interference  effect  is  sufficiently  strong,  it  may 
be  possible  to  ’'tune”  the  rotor-stator  stage  to  minimise  radiation  at  a given 
frequency  by  adjusting  the  inter-row  separation  ix„  - xsl  (sec  liquation  (19)). 

To  achieve  maximum  benefit,  Equation  (28)  indicates  that  i should  ap- 

proximately equal  J -p  j * However,  most  modem  turbofans  are  designed  with 
maximum  rotor-stator  separation,  subject  to  size  and  weight  constraints.  This 
usually  results  in  negligible  potential  interactions  between  rows  (they  decay 
exponentially);  practically  all  the  noise  then  results  from  the  viscous  wake 
impingement  on  the  downstream  row,  i.e.,  in  our  notation,  I I j mj’ 
These  arguments  suggest  that  moving  the  rows  closer  together  until  they  play 
roughl v equal  roles  may  lead  to  noise  reduction  at  s given  frequency,  if  their 
signals  can  be  made  to  destructively  interfere.  This  would  probably  be  useful 
only  in  situations  where  the  acoustic  field  was  dominated  by  a single 
mode.  Whether  the  benefits  of  such  a scheme  could  be  realized  ir.  practice 
would  depend  on  several  factors  not  included  in  the  present  model.  One  would 
also  have  to  consider  if  such  an  optimum  separation  would  have  an  adverse  af- 
fect on  the  other  harmonics  .in  the  spectrum. 

Useful  qualitative  conclusions  regarding  the  variations  in  upstream 
and  downstream  radiated  power  which  result  from  changes  in  axial  Mach  number 
and  frequency  can  also  be  drawn  from  Equation  (28).  The  arguments  below  apply 
to  radiation  from  a single  rotor  as  well  as  from  a rotor-stator  pair,  due  to 
the  absence  of  transmission  and  reflection  effects  of  the  other  row.  Hence 
for  the  sake  of  simplicity  imagine  the  mode  amplitudes  of  one  row,  say  the 
rotor,  to  be  dominant  over  those  of  the  other  row.  First  we  consider  the 
situation  with  no  mean  flow.  In  this  case.  Equation  (27)  shows 
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and  the  power  contained  in  any  one  mode,  call  it  P. 

to  (from  Equation  (28)): 
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Now,  from  Equation  (19)  we  see  that  the  only  thing  differentiating  \-p] 


from 


is  the  sign  of  that  term  in  the  integrand  beginning  with  sin  at 


If  the  integrand  is  not  very  sensitive  to  variations  in  radius,  Equation  (30) 
can  then  be  written  approximately  as 
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This  result  is  in  agreement  with  that  of  Morfey  in  Equation  (09)  of  Reference  6, 
if  there  one  sets  <=  o , interchanges  n and  ^ , and  notes  that  in  his  nota- 
tion, € / tcC.'  equals  cH  &R/(Arr<r)  in  ours.  Equation  (31)  indicates  that  even 
when  * n , equal  amounts  of  power  do  not  radiate  to  either  side  of  the 
blade  row.  This  results  from  the  fact  that  the  relative  angle  between  the 

^ j 

dipole  axis  and  the  propagation  vectors  of  the  acoustic  mode,  say  , is 
different  on  the  two  sides  of  the  blade  row.  The  situation  is  shown 
schematically  in  figure  14. 


The  geometry  of  Figure  14  suggests  that  this  asymmetry  in  radiated 
power  should  disappear  for  modes  propagating  very  near  the  duct  axis,  which 
happens  well  above  cut-or£.  In  this  case 
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and  the  ratio  in  Equation  (31)  approaches  unity. 
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In  fact,  this  can  be  shown  to  be  true  in  the  more  general  case  with 
axial  flow;  Equation  (27)  then  yields, 


0 - r\  U.  ^ 

R = 7 1 (33) 


The  first  factor  in  Equation  (28)  then  becomes 
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and  it  can  be  shown,  again  using  Equations  (19)  and  (32),  that 
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Hence,  while  P^,  ai'id  Pn"t,  will  still  depend  individually  on  Mr  , their  ratio 
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will  always  approach  unity  for  modes  well  above  cut-off,  again  in  agreement 
with  Reference  0.  Care  should  be  taken  not  to  apply  this  result  very  near 
the  condition  of  transonic  axial  flow,  Mx  *=  1 . Here  the  upstream  radiated 
power  will  vanish,  whereas  that  in  the  downstream  direction  remains  finite, 
in  violation  of  Equation  (55).  This  is  because  the  linearized  analysis  on 
which  these  results  are  based  is  invalid  right  at  the  transonic  condition. 


D.  COMPARISONS  OF  THEORY  WITH  EXPERIMENT  AND  DISCUSSION 


This  section  presents  comparisons  between  numerical  predictions  from 
a computer  program  based  on  the  preceding  analysis  and  the  acoustic  data  de- 
scribed in  Section  II.  Unfortunately,  as  pointed  out  in  the  discussion  of 
Section  II,  acoustic  data  could  be  obtained  only  for  the  fourth  and  fifth 
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harmonics  of  blade  passing  frequency.  As  seen  below,  this  seriously  hampers 
one's  ability  to  draw  meaningful  conclusions  concerning  the  model's  validity. 


The  basic  inputs  required  by  the  computer  program  are  the  geometries 
of  the  two  blade  rows,  the  steady  lift  on  each  of  the  blades,  and  the  drag  co- 
efficient of  the  upstream  row,  which  in  the  C3se  of  the  present  experiments, 
is  the  rotor.  Since  at  the  speeds  encountered  in  the  experiment  compressibility 
effects  should  be  negligible,  and  the  hub/tip  ratio  is  large,  a two-dimensional 
incompressible  theory  can  be  used  to  relate  the  lift  and  drag  on  each  blade  to 

the  turning  performance  and  viscous  losses  across  each  row.  For  example, 
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Horloch  has  shown  that  the  drag  coefficient,  C0  , of  each  blade  is  related 
to  the  total  pressure  loss  coefficient,  ACV  , by  the  following 

" AC*r  ^ ^ (36) 


where 
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is  the  mean  flow  angle,  d the  blade  spacing,  c the  semi-chord,  and 
fluid  density.  It  is  important  to  note  that  and  the  total  pressure 
must  be  measured  in  coordinates  fixed  to  the  blade  row.  Since  o^/^. 

being  the  vector  mean  of  the  velocities  upstream  and  downstream  of  the 


v. 


blade  row,  and  o L = 


znr 


Equation  (36)  becomes: 
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where  8 is  the  number  of  blades  and  everything  on  the  right  hand  side  is 
either  known  or  measured. 


Also,  the  lift  coefficient  on  each  blade,  CL  , can  be  obtained  from 
the  turning  performance  through: 
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where  j3t  and  fi  are  the  upstream  and  downstream  flow  angles,  respectively. 
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where  z/  and  v are  the  tangential  flow  velocities,  again  measured  relative 
to  the  row,  and  substituting  for  d , and  C0  from  above,  we  obtain 
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Now  the  lift  force  is  actually  put  into  the  rotor-stater  interaction  noise 
program  in  the  form  of  the  Glauert  coefficients  /)„  , required  by  Osborne's 
theory  (see  Appendix  B) . The  first  two  of  these,  Ro  and  A,  , are  respectively 
the  blade  angle  of  attack  ar.d  the  ratio  of  maximum  camber  to  blade  semi-chord, 
assuming  a circular -arc  profile.  They  are  related  to  JL  by 

Cu  = 2 n (Ra  t fl} ) (43) 

Since  the  angle  of  attack  R # is  simply  the  mean  flow  angle  minus  the  known 
geometric  blade  stagger  angle.  Equations  (42)  and  (43)  are  sufficient  to  de- 
termine flt  ; for  n greater  than  one,  the  Rr  are  assumed  to  be  zero.  In  this 
way  the  Rn  are  computed  in  such  a manner  as  to  use  the  actually  measured 
steady  loads  in  the  Osborne  theory. 


For  the  present  acoustic  experiments  all  quantities  on  the  right 
hand  side  of  Equations  (39) -(43)  are  either  known  or  measured  in  the  case  of 
the  rotor,  for  which  data  were  taken  were  over  a wide  range  of  conditions  in 
the  previous  program.1  Unfortunately,  since  stator  performance  data  were 
obtained  near  rotating  stall  inception,  they  do  not  include  the  range  of  inlet 
conditions  covered  by  the  present  experiments.  Consequently,  stator  turning 
performance  was  inferred  from  design  data  made  available  by  the  manufacturer 
for  the  corresponding  stator  in  the  complete  J-79  engine.  The  outlet  flow 
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angle  at  design  was  assumed  to  apply  for  all  incidence  angles  away  from  the 
design  value  (corresponding  to  a high  solidity  approximation).  The  outlet 
angles  at  the  off-design  stator  stagger  angles  used  here  were  assumed  to 
differ  from  the  design  value  by  the  same  amount  as  the  stagger  angle  differed 
from  the  design  stagger  angle.  Viscous  losses  were  assumed  negligible  for 
the  stator;  this  is  a safe  assumption  since  A has  a relatively  weak  influ- 
ence in  Equation  (42)  and  the  CD  of  Equation  (39)  is  not  needed  for  the  down- 
stream row.  The  rotor  mid-chord  location  was  chosen  as  - 0 . Hence,  the 
coordinates  of  the  piobe  in  Figure  l were  x - -.979  ft.,  r = fr  = 1 2 f ft., 

and  d = 0°. 

Due  to  the  large  hub/tip  ratio  of  about  0.8,  the  dimensionless 
steady  lift  and  drag  coefficients  used  by  the  program  were  input  as  constants 
independent  of  radius.  These  were  computed  from  the  radially  averaged  turn- 
ing and  loss  data  of  Reterence  1.  The  lift  and  drag  forces  do,  of  course, 
vary  with  radius,  because  all  the  velocities  scale  with  the  rotational  speed, 

JT t~  . To  check  the  validity  of  using  the  radially  averaged  coefficients,  one 
case  was  also  run  with  a four-point  radial  distribution  of  these  coefficients 
(also  measured  in  Reference  1).  This  produced  only  minor  differences  in 
the  results. 

Figure  15  compares  the  experimental  measurements  of  r.m.s.  wall  pres- 
sure described  in  section  11  fox  the  fourtli  harmonic  cf  blado  passage  frequency 
with  the  computer  predictions  for  the  mean  stator  staggei  angle  of  37.2°. 

The  axial  flow  velocity  and  rotor  rpm  were  varied  in  such  a way  that  the 
relative  rotor  inlet  angle,  and  hence  also  turning  performance,  of  both*  rows 
was  held  constant  for  all  conditions  shown.  Figure  16  shows  the  same  compari- 
son for  a mean  stator  stagger  angle  of  28.2°. 

The  comparison  of  abso  levels  is  seen  to  be  rather  poor,  although 
the  shapes  cf  both  curves  are  in  general  agreement.  For  example,  in  Figure  15 
both  experimental  and  theoretical  curves  are  characterized  by  two  rath' r well 
defined  peaks,  at  approximately  1170  and  1445  rpm  in  the  experiment,  and  1150 
and  1415  rpm  in  the  theory.  Similar  behavior  is  exhibited  in  Figure  16.  As 
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explained  under  the  discussion  of  Section  II,  di  screte-tone  noise  signals 
could  be  obtained  only  at  the  fourth  and  higher  harmonics  of  blade  passage 
frequency  in  the  accessible  rpm  range  of  the  annular  cascade  experiments. 

The  left  hand  peak  in  the  theoretical  curve  corresponds  to  the  operating 
point  where  the  n - 22  azimuthal  and  m - 1 radial  order  mode  for  the  fourth 
harmonic  of  blade  passage  frequency  exceeds  the  theoretical  cut-off  condition 
and  begins  to  propagate.  The  peak  at  the  higher  rpm  corresponds  to  where  the 
n = 22,  m - 2 mode  begins  to  propagate.  These  peaks  are  a result  of  the 
acoustic  resonance  that  occurs  at  the  cut-off  frequency  of'any  mode,  due  to 
the  vanishing  of  in  the  denominator  of  the  modal  amplitudes  (see 

Equation  (19)) . The  height  of  the  theoretical  curve  is  limited  only  by  the 
precision  of  the  computer.  The  level  and  shape  of  the  experimental  peaks 
in  these  regions  are  heavily  influenced  by  factors  not  included  in  the  ideal- 
ized model,  e.g.,  non-linearities,  viscous  effects,  and  radial  variations  in 
axial  velocity.  In  view  of  this  a small  discrepancy  between  the  locations 
of  the  experiment”)  and  theoretical  peaks,  here  about  20  rpm,  is  to  be  expected. 

Note  that  the  pressure  levels  rose  slightly  from  Figure  15  to 
Figure  16,  probably  a result  of  tilting  the  stator  load  vector  closer  to  the 
angle  of  propagation  of  the  acoustic  mode,  as  well  as  a modest  increase  in 
flow  velocity.  Again,  it  is  encouraging  to  note  in  this  regard  that  theory 
and  experiment  rose  by  roughly  the  same  amount 

Figure  17  contains  limited  data  for  the  fifth  harmonic  noise  of  the 
same-  rotor-stator  pair  as  Figure  lo.  The  fifth  harmonic  excites  propagating 
modes  at  a much  lower  rpm.  In  this  case  only  the  n =14,  - 1 mode  is 

present,  the  rpm  at  which  it  theoretically  exceeds  cut-off  being  just  to  the 
left  of  that  shown  in  the  figure.  Here,  the  shapes  of  the  two  curves  are 
only  marginally  similar,  and  the  discrepancy  in  absolute  levels  is  even  wider 
than  in  Figures  IS  and  16. 

This  discrepancy  in  absolute  level  between  theory  and  experiment 
bears  further  discussion.  For  the  experimental  rotor-stator  Stage  having 
46  blades  in  the  rotor  and  54  blades  in  the  stator,  the  n - 22,  m - 1, 


33 


2 modes  in  Figures  i5  and  16  are  generated  by  the  third  loading  harmonic,  on 
the  rotor  in  concert  with  the  fourth  loading  harmonic  on  the  stator  (.see  Equa- 
tions (19c,  dj).  Moreover,  computer  results  indicate  that  at  these  conditions 
the  potential  interactions  between  rows  are  negligibly  small  in  comparison 
with  the  viscous  interaction.  Hence,  the  question  is  reduced  to  how  well  tne 
model  can  predict  the  noise  resulting  from  the  fourth  loading  harmonic  on  the 
stator  induced  by  its  passage  through  the  rotor  viscous  wakes.  (Figure  17 
involves  the  fifth  harmonic  of  stator  loading,  and  hence  is  an  even  more 
severe  test  of  the  model.)  Since  the  velocities  involved  are  all  low  enough 

to  be  considered  incompressible,  we  are  essentially  dealing  with  the  original 
2 3 

Kemp-Sears  model  ’ to  which  Osborne's  analysis  reduces  in  this  limit.  It  is 
felt  that  the  following  are  the  most  likely  sources  of  the  discrepancy: 

• Poor  modeling  of  the  viscous  wake  defect  at  the  higher  harmonics 

• Inaccuracies  in  turning  performance  and  loss  data  used  as  inputs 

• Inadequacy  of  strip-wise  application  of  isolated  airfoil  theory 
in  predicting  unsteady  loads 

Each  of  these  is  discussed  individually  below,  beginning  with  the  last. 

The  ad  hoc  application  of  a strip  theory  aerodynamic  analysis  for 
predicting  unsteady  loads  in  a three-dimensional  annular  flow  is  open  to 
question,  although  the  large  hub/tip  ratio  and  large  number  of  blades  in  this 
case  would  seem  to  justify  it.  Also,  since  the  solidity  ratio,  2c/d  is 
near  unity,  the  use  of  isolated  airfoil  theory  in  estimating  these  loads 
introduces  some  error.  These  questions  are  particularly  crucial  with  respect 
to  phase  variations  in  loading  along  the  span,  as  these  critically  affect 
the  radial  integration  involved  in  the  modal  pressure  amplitudes,  -p  in 
Equation  (19).  It  is  difficult  to  assess  the  importance  of  these  approxima- 
tions in  the  absence  of  any  experimental  data  on  the  load  fluctuations,  or 
a truly  three-dimensional  annular  cascade  theory  for  unsteady  flow  through  a 
blade  row. 
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The  inaccuracy  involved  in  estimating  the  steady  loading  and  viscous 
losses  on  the  rows  is  also  difficult  to  assess.  The  principle  reason  is  that 
such  data  as  were  obtained  were  measured  with  the  row  operating  in  isolation, 
since  at  the  time  their  use  in  the  present  investigation  was  not  envisioned. 
Hence,  the  performance  of  both  rows  may  have  been  somewhat  different  from 
that  assumed.  This  could  significantly  affect  the  potential  interactions. 
Since  these  were  found  to  be  negligible  as  noise  sources  at  these  conditions, 
a strong  variation  in  loading  would  probably  be  needed  to  appreciably  change 
the  acoustic  predictions.  The  prediction  of  viscous  wake  interaction,  how- 
ever, might  have  been  considerably  improved  if  more  detailed  knowledge  of  the 
blade  wake  structure  had  been  available.  This  is  discussed  at  more  length 
below. 


With  regard  to  the  wake  modeling,  both  Osborne's  analysis  and  the 
previous  work  of  Kemp  and  Sears  made  use  of  empirical  laws  for  the  wake  width 
and  velocity  defect  suggested  by  Silverstein  et  al,'--5  based  on  isolated  air- 
foil data.  In  particular  the  decay  of  the  wake  centerline  dynamic  pressure, 
q^  . is  assumed  to  follow 


±c_  , _ 

°°  - 0.4 


(44) 


where  ^ is  the  free-stream  dynamic  pressure,  x the  distance  downstream  of  the 
blade  mid-chord  position,  and  c the  semi-chord,  from  this,  the  centerline  ve- 
locity defect,  l'„c  , was  inferred,  assuming  it  to  be  a small  perturbation  of 
the  free-stream  flow.  In  the  present  experimental  configuration,  there  is 
roughly  a one  semi- chord  separation  between  _ . rotor  trailing  edge  and  the 
leading  edge  of  the  stator,  liquation  (39),  in  conjunction  with  the  loss  data, 
indicates  a rotor  CD  of  0.16.  For  such  a high  drag  coefficient.  Equation  (44) 
predicts  a negative  value  of  dynamic  pressure,  and  is  clearly  inapplicable. 


Partially  for  this  reason,  and  also  because  it  was  feit  that  the 
radially  averaged  loss  data  were  perhaps  unduly  influenced  by  viscous  losses 
at  the  inner  and  outer  duc.t  walls,  the  value  of  CD  used  in  the  calculations 
presented  here  was  lowered  to  0.10.  To  indicate  the  uncertainty  involved  in 
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estimating  uc  , an  alternate  correlation  proposed  by  Dittmar  (Reference  17 
liquation  (".,),  ar.d  based  in  part  : . the  cascade  data  of  Lieblcin  and  Roudebush , 
was  evaluated  for  the  same  conditions.  It  indicated  an  increase  in  o-c  by  a 
factor  of  about  1.6.  Since  the  unsteady  loading  is  proportional  to  u-c  , and 
only  one  mode  is  propagating  over  most  of  Figures  IS  and  16,  one  would  expect 
this  to  raise  the  theoretical  SPL  by  about  4 dB. 


In  addition  to  this  uncertainty,  a surprising  variation  was  noted  in 
the  calculations.  That  is,  when  C0  was  lowered  from  .16  to  .10,  there  was  a 
rise  in  the  predicted  viscous  interaction  force,  and  hence  duet  wall  pressure, 
in  sharp  contrast  to  the  conclusion  of  Kemp  and  Sears'5  that  such  loading  ap- 
pears to  be  linearly  proportional  to  cD  . The  observed  behavior  can  be  traced 
to  an  exponential  factor  in  the  load  prediction  (sec*  Appendix  B (c))  whose 
argument  is  propoitional  to  minus  the  product  of  Cp  times  the  square  of  the 
load  harmonic  index.  Clearly,  for  the  high  drag  coefficient  and  load  harmonic 
index  of  interest  here,  this  exponential  decay  dominates  any  linear  variation. 


The  above  anomalous  behavior  strongly  suggests  that  the  viscous  wakes 
are  not  modeled  accurately  enough  to  predict  the  higher  loading  harmonics  on 
the  downstream  row  when  the  rotor  and  stator  arc  in  close  proximity  to  one 
another,  and  the  wake  structure  is  accordingly  changing  rapidiy.  Table  I 
further  supports  this  contention.  The  first  two  entries  in  the  first  column 
give  the  velocity  defect  profile  assumed  by  Kemp-Sears-Osborne,5’ 24  and  that 
originally  suggested  by  Silverstein  et  al,  as  a fit  to  their  isolated  airfoil 
data.  Here  u is  the  velocity  defect  at  any  distance  y from  the  wake  centerline, 
and  Y is  the  wake  half-width.  The  last  entry  is  suggested  by  the  present 
authors  as  fitting  the  Silverstein  data  nearly  as  well,  at  least  within  the 
experimental  scatter  of  about  i.10  in  u/ue  shown  in  Figure  42  of  Reference  .33. 
All  three  are  normalized  in  such  a way  as  to  produce  the  same  total  momentum 
defect  in  the  wake  flow.  The  second  column  gives  the  relative  magnitude  of 
the  ith  harmonic  in  the  Fourier  expansion  of  an  infinite  train  of  such  pulses 
spaced  a distance  o-5  apart.  Since  this  is  the  only  change  that  would  be  in- 
volved if  a different  profile  were  assumed,  the  influence  cn  the  theoretical 
predictions  can  be  easily  assessed.  The  first  entry  is  the  exponential  discussed 
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above,  tut  written  m tents  ot'  Y rather  than  Ca  for  the  sake  of  simplicity. 

The  next  two  entries  show  how  this  factor  would  change  if  one  or  the  other  of 
tiie  substitute  profiles  had  been  assumed.  Th“  resulting  changes  in  the  pre- 
dictions are  shown  in  the  last  column  for  V/a,  wa,  (actually  V/a. „ w»  a, 
for  the  present  case  of  a rotor-stator ) = .265.3  appropriate  to  the  present 
conditions,  and  the  first,  fourth,  and  tenth  loading  harmonics.  Note  that  the 
first  harmonic  is  relatively  insensitive  to  the  choice  of  profile,  but  at 
higher  values  of  2 such  as  needed  here,  the  choice  is  critical.  Had  the  last- 
profile  been  used  in  the  calculations,  in  conjunction  with  the  alternate 
estimate  ot  u4  discussed  above,  it  is  seen  that  theoretical  curve  in 
Figure  15  would  be  almost  12  dli  higher. 

It  would  be  presumptuous  to  suggest  that  this  constitutes  a better 
model  for  the  wake  properties  in  the  present  case,  or  any  other,  simply  because 
it  correlates  better  with  our  acoustic  data.  But  it  is  cleai  'hat  the  prediction 
of  iiigher  harmonics  is  too  sensitive  to  the  details  of  wake  structure  to  be  re- 
liable, and  for  these  we  have  had  to  rely  on  isolated  airfoil  data.  Moreover, 
recent  experimental  work2'-1  ’ ^ ^ has  shown  that  the  wake  properties  for  cascad- 
ed airfoils  and  rotor  blades  can  differ  significantly  from  those  for  ar:  isolated 
ait  foil.  Hence  it  is  felt  likely  that  the  theory  could  benefit  from  a refined 
viscous  wake  model  which  better  reflects  the  true  situation.  Also,  if  they 
should  prove  possible,  comparisons  of  the  theory  against  acoustic  data  for  the 
lower  harmonics  in  the  spectrum  would  be  vex>  valuable.  These  are  loss  sen- 
sitive to  detailed  changes  in  wake  structure,  and  would  allow  a better  assess- 
ment of  the  validity  of  other  portions  of  the  theoretical  model. 

E.  ADDITIONAL  NUMERICAL  RESULTS 

In  addition  to  the  above  comparisons  with  experiment,  which  all  in- 
volve incompressible  flow,  some  calculations  have  been  performed  to  show  the 
influence  of  compressibility  on  the  radiated  noise.  One  of  the  test  cases 
presented  by  Clark,  et  al.  (Reference  1 9 , Vol.  Ill,  p.  59)  involved  tne  cal- 
culation of  the  upstream  pressure  amplitudes  generated  by  the  potential  inter- 
action on  the  upstream  s?  tor  of  an  inlet  stator-rotor  pair.  The  present 
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program  was  run  using  the  same  input  conditions  as  given  there  (see  test 
case  #2,  subroutine  package  AAR/Jt) . Briefly,  the  stator  and  rotor  were  housed 
in  an  annular  duct  with  hub/tip  ratio  of  .35,  and  axial  Mach  number  of  0.5. 

The  rotor  (15  blades)  had  a tip  Mach  number  of  .875  and  was  loca.ted  a distance 
rT  /$  downstream  of  the  stator  (10  vanes).  In  this  case  the  potential  inter- 
action on  the  upstream  stator  was  dominant.  The  amplitudes 
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as  used  by  Clark,  et  al.*9  are  compared  wit.,  the  present  predictions  in 
Table  2. 


Table  2 


COMPARISON  OF  MODAL  AMPLITUDES  COMPUTED  BY  CLARK  ET  AL.  19 

WITH  PRESENT  THEORY 


Clark  et  al . 19 

Present  Theory 
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-2 

i 

3 

7.46  X 10 
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1 

' , 

1 .34  x ID"2 

-3  i 

5.36  x 10 

1 

- 1 i 

: +5 

1 

2 

S.07  x 10  ‘ 

1.91  x 10  ' 

-2 

3 

b . 94  x 10 

3.84  x 10  " 

This  particular  calculation  in  Reference  19  also  assumed  a line 
dipole  model  for  the  blades,  and  so  the  differences  in  the  predicted  values 
are  largely  a result  of  using  the  Osborne  rather  than  Kemp- Scars  aerodynamic 
analysis.  As  one  might  expect,  the  inclusion  of  compressibility  (the  relative 
Mach  number  at  the  stator  was  about  0.5)  had  a noticeable  effect,  generally 
decreasing  the  modal  amplitudes. 
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In  the  acoustic  experiments  performed  under  the  current  program,  the 
fundamental  blade  passing  tone  could  not  propagate  owing  to  power  constraints 
in  the  annular  cascade  (see  Section  II).  However,  calculations  indicate  that 
it  should  begin  to  propagate  above  roughly  1700  rpm.  To  study  its  behavior, 
and  give  some  idea  of  how  the  theoretical  predictions  vary  over  a wider  range 
of  conditions  than  those  above,  the  calculations  in  Figure  18  are  presented. 

The  inputs  are  the  same  as  in  Figure  IS,  except  that  the  rpm  and  Mx  were  in- 
creased, again  proportionately,  into  the  subsonic  regime.  Thus,  when  this 
tone  first  begins  to  propagate  at  the  left  »f  the  figure,  the  relative  Mach 
numbers  of  the  rotor  and  stator  are  approximately  .IS  ar.d  .11,  respectively. 

When  8700  rpm  is  reached,  these  are  .77  and  .56,  and  the  first  four  radial 
modes  of  the  n = -8  circumferential  node  are  propagating.  The  rpm  at  which 
each  of  these  starts  to  propagate  is  shown  at  the  bottom  of  the  figure. 

The  upper  curve  gives  the  variation  in  SPL  at  the  same  location  as 
before.  As  e-xjectod,  there  is  a large  increase  in  signal  over  the  fourth 
harmonic  of  Figure  15.  The  lower  curve  in  Figure  18  is  the  total  power  radiated 
in  the  upstream  direction.  Straight  line  segments  have  been  used  to  connect 
the  computed  points  for  ease  of  visualisation  only,  ar.d  should  not  be  assumed 
to  accurately  depict  the  true  curve.  The  curves  exhibit  a rise  just  above  each 
cut-off  frequency,  due  to  the  acoustic  resonance  for  tiiat  mode,  as  discussed 
previously;  again  it  should  be  emphasised  that  the  levels  computed  very  near 
these  frequencies  arc  oper.  to  question  because  they  represent  singular  points 
in  the  solution.  In  fact,  the  rise  becomes  harder  to  resolve  as  rpm  is  in- 
creased. This  results  because  as  increasingly  higher  order  radial  msecs  get 
turned  on,  the  integrand  in  the  radial  integration  in  liquation  (19)  oscillates 
more.  The  value  of  the  integral  itself  thus  approaches  zero,  and  so  one  must 
be  extremely  close  to  the  cut-off  frequency  before  the  pressure  amplitude  begins 
to  diverge. 

The  behavior  of  the  SPL  curve  in  Figure  id  is  much  more  unpredictable 
than  that  in  Figure  15.  This  i.s  because  over  most  ol  the  rpm  range  more  than 
one  mode  is  prepagat ing , and  so  the  interference  between  modes  described  by 
liquation  (25)  comes  into  play.  No  doubt  at  some  other  point  in  the  duct  this 
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curve  would  look  quite  different,  due  to  its  complicated  dependence  on  observer 
position  and  frequency.  At  the  higher  rpm  values  the  curve  begins  to  snoot, 1 
out  somewhat,  probably  because  as  the  signal  is  distributed  amongst  more  modes, 
the  constructive  and  destructive  interferences  tend  to  cancel  one  another. 

The  curve  representing  total  radiated  power  is  much  more  well  behaved, 
since  it  represents  an  integral  over  the  duct  cross  section,  thus  removing  any 
dept .'.dence  on  position.  With  the  exception  of  the  resonance  peaks,  it  exhibits 
a smooth  rise  as  relative  velocity  is  increased,  as  expected.  The  interference 
effects  between  waves  from  the  rotor  and  stator  described  by  Equation  (28)  do 
not  arise,  since  the  computer'  results  again  indicate  that  -p  “ >> 

The  leveling  off  and  subsequent  decay  of  the  curve  in  the  higher  rpm  region 
is  not  surprising.  Upstream  radiated  power  must  go  to  zero  as  a sonic  axial 
flow  velocity  is  approached;  hence  the  presence  of  a maximum  in  the  subsonic 
regime  is  to  be  expected.  However,  as  Osborne's  theory  becomes  more  approxi- 
mate as  the  relative  Mach  number  is  increased,  whether  the  location  of  the 
peak  is  accurately  predicted  is  open  to  question.  Subsequent  to  the  completion 
of  the  present  work,  Ajniet'-'7  presented  an  approximate  analysis  of  compressible 
unsteady  airfoil  theory  which  also  makes  use  of  asymptotic  techniques.  His 
analysis  complements  Osborne's  in  that,  it  becomes  more  exact  as  the  Mach  number 
and/or  reduced  frequency  increases,  and  so  should  prove-  useful  in  the  higher 
subsonic  regime. 

The  computing  time  needed  to  generate  the  numerical  results  in 
Sections  ili-D,  E above  on  an  IBM  370/103  varied  between  approximately  1 and 

10  seconds  per  point,  depending  primarily  on  the  numb  .-r  of  propagating  modes. 

* 

I'.  CONCLUDING  REMARKS 

A compressible  two-dimensional  analysis  of  the  aerodynamic  forces 
generated  by  rotuT -stator  interaction  has  been  matched,  on  a strip-theory  basis, 
to  the  well  known  three-dimensional  annular  duct  acoustic  modes.  Expressions 
have  been  derived  for  both  ho  mean  square  pressure  at  any  position  in  the  duct 
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and  the  total  power  radiated  at  harmonics  of  blade  passage  frequency.  The 
mean  square  pressure  signal  is  seen  to  be  subject  to  interferences  between 
the  various  propagating  modes.  The  total  radiated  power  is  the  sum  of  inde- 
pendent contributions  from  each  mode,  but  is  still  subject  to  interference  ef- 
fects between  waves  generated  by  each  row  within  the  same  mode.  These  are 
evaluated  easily  on  a computer  in  terms  of  stage  operating  conditions  and 
geometry . 


Experimental  acoustic  data  were  taken  in  the  annular  cascade  facility, 
as  described  under  Section  II,  to  provide  a means  of  assessing  the  validity 
of  the  theoretical  model.  Data  on  blade  row  steady- state  performance  and 
losses  had  been  taken  in  the  annular  cascade  facility  under  a previous  program. 
These  were  used  as  inputs  to  the  theory  to  calculate  SPL  at  the  outer  wall, 
and  the  results  were  compared  with  the  acoustic  data  for  two  configurations  of 
a rotor-stator  pair.  Unfortunately,  power  constraints  in  the  rig  allowed  only 
tiie  fourth  and  higher  harmonics  to  be  studied,  which  is  felt  to  be  a rather 
severe  test  of  the  theory.  Although  it  was  able  to  predict  cutoff  frequencies 
and  follow  trends  in  the  data  adequately,  the  theory  significantly  under- 
estimated absolute  levels. 

The  most  likely  source  for  the  discrepancy  is  felt  to  be  in  the  viscous 
wake  modeling,  which  at  present  is  based  on  isolated  airfoil  data.  Recent 
experimental  work  by  others  has  suggested  that  the  structure  of  cascade  and 
rotor  wakes  can  be  significantly  different,  and  this  is  expected  to  have  a 
strong  influence  on  the  higher  harmonics.  It  is  felt  that  the  theoretical 
model  could  benefit  from  a closer  examiriat  on  of  such  cascade  or  rotoi  wake 
data  as  are  available,  as  well  as  an  improved  aerodynamic  load  model  at  high 
reduced  frequencies  and  Mach  numbers. ^ It  is  also  hoped  to  obtain  more 
acoustic  data  in  the  annular  cascade  facility,  which  has  been  repowered  and 
will  now  allow  the  propagation  of  lower  harmonics  in  the  spectrum. 
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SECTION  IV 

D1RCCT  LUTING  SURFACE  THEORY  FOR  A COMPRESSOR  ROTOR 


Increased  emphasis  on  reduction  of  the  size,  weight,  and  noise  out- 
put of  axial  flow  turbomachinery  demands  improved  understanding  of  the  flow 
through  high-speed  fan  and  compressor-  blade  rows.  As  more  detailed  questions 
are  asked  about  modern  blade  row  performance,  the  essentially  three-dimen- 
sional character  of  the  flow  takes  on  increased  importance.  The  task  of 
calculating  the  fully  nonlinear,  three-dimensional , viscous  flow  through  inter- 
acting blade  rows  is  a formidable  one  indeed.  Consequently,  some  approxima- 
tions are  required  in  order  to  obtain  a tractable  model,  the  most  familiar 
being  the  idealization  of  inviscid  flow  through  a two-dimensional  cascade. 

A linearized  analysis  of  the  steady,  inviscid,  three-dimensional  flow  through 
an  isolated  rotor  contains  important  features  not  present  in  the  ct rresponding 
two-dimensional  cascade  approximation.  For  example,  though  restricted  to 
lightly  loaded,  thin  blades,  it  does  include  disturbances  induced  by  the 
trailing  vortex  wakes  which  result  from  spanwise  variations  in  the  blade 
circulation.  Also,  modem  fans  and  compressors  often  operate  with  a transonic 
inflow  that  is  subsonic  over  the  inner  portion  of  the  blade  span  and  super- 
sonic over  the  outer  portion.  At  Mach  numbers  where  shock  losses  arc  small, 
the  linearized,  three-dimensional  analysis  remains  valid,  but  it  does  not 
have  a well-defined  two-dimensional  limit." 

The  small-perturbation  approach  to  three-dimensional  compressor 
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flows  was  pioneered  by  McCune"  who,  in  the  spirit  of  linearized  wing  theory, 
separated  the  thickness  and  loading  contributions  to  the  rotor  disturbance 
field.  In  his  original  papers,  McCur.e  treated  the  thickness  problem  for  the 
subsonic,  transonic,  and  supersonic  flow  regimes.  Later,  Okurounmu  and 

McCune  employed  a vortex  representation  of  the  blade  row  to  solve  the 
indirect  lifting  problem,  i.e.,  that  in  which  t he  disturbance  field,  along 
with  the  blade  geometry  needed  to  produce  it,  are  determined  from  a prescribed 
distribution  of  blade  loading.  Lordi  ' has  investigated  the  acoustic  im- 
plications of  both  of  these  analyses;  the  sound  radiation  from  an  isolated 
rotor  operating  at  supersonic  tip  speeds  was  evaluated  for  various  thickness 
and  loading  distributions. 
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While  the  thickness  and  loading  contributions  to  the  rotor  disturbance 

field  can  be  treated  separately  in  the  linearized  analysis,  they  are  not 

entirely  independent,  in  order  that  there  be  no  loading  contribution  to  the 

flow  field  associated  with  a given  distribution  of  blade  thickness,  this 

thickness  must  be  distributed  about  an  unknown  camber  line.  The  camber  lines 

required  for  the  rotor  to  be  unloaded  must  be  determined  as  part  of  the 

solution  to  the  thickness  problem  , in  much  the  same  way  as  the  camber  lines 

corresponding  to  a specified  loading  distribution  are  computed.  Erickson, 
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Lordi,  and  Rae  presented  results  for  thickness-induced  camber  lines  at  high 
subsonic  tip  speeds.  They  also  have  computed  the  camber  lines  reo.uired  to 
produce  given  loading  distributions,  as  have  Okurcunmu  and  McCune  ^ 

A major  objective  of  the  preset.-  program  has  been  to  develop  a 
linearized  analysis  of  the  direct  lifting  surface  problem  for  the  compressible, 
three-dimensional  flow  through  a rotor  row.  In  contrast  to  the  indirect  load- 
ing problem,  the  direct  problem  refers  to  the  situation  where  the  blade  inci- 
dence and  camber  lines  are  given,  and  the  resulting  blade  loading  must  be 
determined.  Once  the  blade  loading  is  known,  the  evaluation  of  the  rotor 
flow  field  is  accomplished  in  the  same  way  as  in  the  indirect  case.  The 
solution  of  the  direct  problem,  together  with  the  existing  solution  of  the 
thickness  problem,  would  permit  computation  of  the  aerodynamic  and  acoustic 
performance  for  a rotor  row  of  given  geometry  ac  specified  operating  conditions. 
In  addition  to  providing  a predictive  capability  for  steady  flow  at  off-design 
conditions,  the  successful  analysis  of  the  direct  lifting  surface  problem  pro- 
vides the  basis  for  examining  three-dimensional  flows  which  are  unsteady  in 
rotor  coordinates. 

In  the  thickness  and  indirect- loading  analyses,  the  flow  field 
quantities  are  found  by  evaluation  cf  integrals  which  express  the  superposi- 
tion of  the  appropriate  singularity  solutions  e.g.,  sources  in  the  thickness 
cast  and  vortices  in  rhe  loading  case.  However,  direct  lifting  surface 
calculations  entail  t.he  solution  of  an  integral  equation,  which  for  a ducted 
rotor  has  a very  complicated  kernel  function.  Some  progress  on  the  direct 
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loading  problem  has  been  reported  by  other  investigators.  Dharwadkar  and 
MeCunc  ' have  obtained  a solution  using  a lifting- line  approximation.  Namba^ 
has  reported  a direct  lifting  surface  analysis  including  some  numerical 
examples.  Nmnba’s  approach  closely  parallels  that  presented  here  and  frequent 
reference  to  it  is  made  throughout  the  report. 

Our  initial  approach  to  the  direct  lifting  surface  problem  was  to 

generalise  the  formal  relationships  between  the  blade  camber  lines  and  the 

blade  loadings  which  were  derived  i.n  Reference  46.  In  this  way,  an  integral 

equation  for  the  blade  loading  was  obtained  in  the  framework  of  the  vortex 
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theory  of  Okurounmu  and  McCune.  “ This  formulation  of  the  direct  lifting 
surface  theory  produced  an  integral  equation  which  not  only  contained  the 
expected  complexities  in  the  kernel  function,  but  moreover  would  require 
repeated  evaluations  of  both  ordinary  and  modified  Bessel  functions  of  large 
order.  The  appearance  of  both  kinds  of  Bessel  functions  is  a consequence 
of  the  way  in  which  the  trailing  vortex  waxes  are  handled  separately  from 
the  bound  vorticity  at.  the  blade  surfaces. 

At  this  point,  an  alternative  formulation  of  the  lifting  surface 
problem  was  sought  by  using  a pressure  dipole  rather  than  a vortex  representa- 
tion of  the  blade  surfaces.  In  principle,  the  two  approaches  should  result  in 
the  same  governing  integral  equation  for  the  blade  loading.  However,  it  was 
felt  that  an  alternative  derivation,  as  is  often  the  case  in  potential  theory 
for  complicated  geometries,  would  yield  an  alternative  form  for  the  result. 

A simpler  form  did  indeed  result  from  the  pressure  dipole  approach,  the  main 
difference  being  the  absence  of  the  modified  Bessel  functions  of  large  order 
which  appear  in  the  vortex  formulation. 

At  about  the  time  we  began  to  pursue  the  pressure  dipole  representa- 
lion,  N'amba  ’ s initial  report  appeared.  There  arc  important  differences 
between  our  formulation  and  his.  Furthermore,  the  initial  sulution  for  the 
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pressure  field  which  we  obtained  with  the  dipole  representation  failed  to 
yield  correct  limits.  In  particular,  the  pressure  rise  between  points  far 
upstream  and  downstream  of  the  rotor  did  not  agree  with  that  from  the  vortex 
theory,  or  with  Mamba's  version  of  the  dipole  formulation.  These  discrepan- 
cies prompted  a complete  review  of  our  analysis.  A crucial  revision  has  been 
made  in  the  pressure  dipole  solution  which  led  to  a limiting  pressure  rise 
across  a lifting  rotor  that  now  agrees  with  the  vortex  theory  result,  but 
still  differs  from  Namba's. 

In  the  material  which  follows,  a detailed  account  is  given  of  the 
linearized  solution  for  the  three-dimensional,  compressible  flow  through  a 
rotor.  This  formalism  is  presented  not  only  because  of  its  relevance  to 
establishing  the  correct  integral  equation  in  a direct  lifting  surface  theory, 
but  also  because  of  its  general  importance  to  other  aspects  of  the  flow 
through  rotating  blade  rows.  In  Section  IV-A,  the  derivation  of  the  line- 
arized equations  is  reviewed  ar.d  a formal  solution  of  them  is  obtained  basec 
on  Green's  theorem.  In  Section  IV-B  solutions  are  found  for  the  disturbance 
fields  of  a point  source  of  mass  and  a pressure  dipole.  Next  these  singu- 
larity solutions  are  used  as  the  Green's  functions  to  determine  the  flow- 
fields  produced  by  rotor  thickness  and  loading  distributions.  Discussion 
of  the  thickness  contribution  is  included  here  because  cf  its  role  in  deriving 
a correct  solution  for  the  entire  flow  field,  and  also  to  point  out  omissions 
in  earlier  results.  In  the  respective  presentations  of  the  thickness  and 
loading  solutions  in  Sections  IV-C  and  IV-D,  the  flow  fields  are  shown  to 
display  the  correct  discontinuous  behavior  across  the  blade  surfaces  and 
trailing  vertex  wakes.  Demonstrations  that  the  results  for  a mass  source, 
a pressure  dipole,  and  the  thickness  and  loading  problems  satisfy  global  mass 
and  momentum  balances  appropriate  to  the  non  inertial  reference  frame  are 
collected  in  Appendix  C.  In  the  final  part  of  the  main  text.  Section  IV-E, 
the-  governing  integral  equation  for  the  direct  lifting  surface  theory  is 
derived  and  the  progress  made  on  its  solution  is  reported. 
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A. 


DERIVATION  OF  MODEL  EQUATIONS  AND  GREEN'S  FUNCTION  SOLUTION 


In  this  section,  the  linearized  equations  and  a formal  integral 
representation  of  their  solution  are  developed  for  the  flow  through  an  isolated 
rotor  row  in  an  infinitely  long  annular  duct.  The  key  assumptions  in  the 
analysis  are  that  the  undisturbed  axial  velocity  is  uniform  and  subsonic, 
and  that  the  disturbance  field  of  the  rotor  is  a small  perturbation  about 
the  resulting  helical  inflow  seen  by  an  observer  in  blade-fixed  coordinates. 
Both  compressibility  and  three-dimensional  effects  are  included.  While 
not  a fundamental  restriction  in  the  analysis,  attention  is  confined  to 
subsonic  relative  tip  speeds.  The  required  extensions  to  supersonic  tip 
speeds  arc  indicated  at  several  points  in  the  derivations.  As  a consequence 
of  the  linearization,  the  blade  surface  boundary  conditions  can  be  separated 
into  thickness  and  camber  line  contri'r  utions  and  their  associated  flow  field 
solutions  superimposed  to  find  the  overall  disturbance  field.  This  separa- 
tion is  effected  by  requiring  the  rotor  blades  to  be  local i>  unloaded  in 
the  thickness  case  and  by  assuming  that  the  blades  have  vanishing  thickness 
in  the  loading  case. 

The  geometry  of  the  blade-fixed  coordinates  is  illustrated  in 
Figure  19  for  a rotor  rotating  in  the  negative  9 direction  with  angular 
velocity  aj  . The  full  nonlinear  equations  for  the  flow  through  such  a blade 
row  in  a cylindrical  coordinate  system  fixed  to  the  rotor  have  been  given  by 
Wu.  ^ These  equations  can  be  linearized  by  writing  the  velocity  in  blade- 
fixed  coordinates  as 

W = •+  V (45] 

where  Uu  is  the  undisturbed  velocity, 

UA  - U 2 2 + cdree  (46) 

and  v is  the  perturbation  velocity  with  components  vr  , vg  , ira  . The  fluid 
pressure  and  density  arc  expanded  in  a similar  way;  the  undisturbed  quantities 
are  denoted  by  ^ and  p„  , the  perturbation  quantities  by  and p . Sub- 
stituting these  definitions  into  Wu's  equations  and  linearizing  the  result  by 
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neglecting  the  products  of  perturbation  quantities  leads  to  the  following 
vector  forms  of  the  stead)'  flow  continuity  and  momentum  equations,  to  first 
order. 

U*  ■ VP  + Poo  ( V ■ v)  - o (47) 

( t?  • V ) + ( UR  • V)  v + X <x>  ( 7,r  - i/r  CQ)  - - — V -p  t4gj 

p*> 

Paying  careful  attention  to  the  proper  differentiation  of  the  un  ,.t  vectors, 
the  component  form  of  the  momentum  equation  can  be  expressed  as 
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where  UR  = and  5 is  the  coordinate  along  the  helical  undisturbed 

streaml incs . 


From  this  basic,  set  of  conservation  equations,  we  can  develop  the 
governing  equations  for  either  the  perturbation  pressure  or  the  velocity 
potential,  the  solutions  of  which  can  then  be  used  to  obtain  the  remaining 
flow  field  quantities.  First,  the  9 and  2 components  of  the  momentum  equa- 
tions are  cast  in  terms  of  distance  along  the  undisturbed  streamlines  and 
the  direction  normal  to  the  streamline  and  radial  directions.  Unit  vectors 
in  these  directions,  which  are  illustrated  in  Figure  20,  are 

S « e0  ^ iff  < ea  c*,  y (sou) 

" = eB  coo  w - AW  y (50b) 

where  ifr  - tan  ' ujr/  U . The  directional  derivatives  along  5 and  n are 
related  to  the  partial  derivatives  with  respect  to  3 and  & by 
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With  these  results  the  $ and  n components  of  the  linearized  momentum  equations 
become 
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This  form  of  the  linearized  momentum  equations  is  useful  in  relating  results 
for  the  pressure  and  velocity  fields.  For  subsonic  flow  where  disturbances 
decay  far  upstream,  Eq.  (52a)  can  be  integrated  along  the  streamlines  to 
obtain 

= 'P»  vs  (53  J 

The  normal  momentum  equation,  Eq.  (52b),  plays  a central  role  in  the  direct 
lifting  surface  theory.  It  contains  the  upwash  velocity,  , which  is 
related  to  the  blade  camber  line  in  order  to  derive  the  integral  equation 
for  the  blade  loading.  The  above  equations  have  the  same  appearance  as  their 
counterparts  in  isolated  airfoil  theory  because  the  linearized  substantial 
derivative,  Ua  , is  independent  of  radius.  It  s ould  be  noted  that  the 
nonir.ertia!  acceleration  terms  enter  only  through  the  curvature  of  the  un- 
disturbed streamlines. 


The  governing  equation  for  the  perturbation  pressure  is  derived  by 
introducing  the  assumption  that  the  disturbance  flow  is  isentropic  so  that 
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where-  a„  is  the  undisturbed  sound  speed.  Using  this  relationship  to  eliminate 
the  density,  and  then  comb’ning  Up  of  the  linearized  continuity  equation 
with  the  divergence  of  the  momentum  equation  leads  to 

j,  U d p 
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There  are  several  ways  to  demonstrate  that  the  velocity  potential 
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satisfies  the  same  governing  differential  equation.  McCune,  starting  in 
flui  '-fixed  coordinates  where  the  disturbance  potential  satisfies  the  wave 
equation,  obtains  this  result  by  transformation  to  the  blade-fixed  system. 

cn 

IVu  u derives  the  nonlinear  potential  equation  which,  when  linearized,  reduces 
to  the  same  result.  Here  we  note  that  the  linearized  momentum  equations, 
in  the  streamline-normal  coordinates,  indicate  that  the  velocity  components 
are  proportional  to  the  gradient  of  the  integral  of  the  pressure  along  the 
undisturbed  streamlines,  a scalui  function.  Accordingly,  n a scalar  velocity 
potential,  defined  such  that  v * V (f  , is  introduced  into  Eqs.  (47)  and  (52a), 
the  results  may  be  combined  with  Eq.  (54)  to  obtain 

a Jo  dX(P 
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This  equation  is  the  same  as  that  given  for  -p  in  Eq.  (55),  and  when  expressed 
in  cylindrical  coordinates  leads  to  the  same  form  of  the  governing  equation 
for  the  velocity  potential  that  McCune  obtained. 

The  formal  solution  of  the  governing  differential  equation  for  d! 
(or  -p  ) is  derived  from  Green’s  theorem,  which  states  that51  for  two  scalar 
functions  <f  and  G , 
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where  the  surface -CT  encloses  the  volume  V and  V is  the  normal  to  the  surface 
S directed  into  V . The  following  development,  while  carried  out  for  the 
velocity  potential  applies  as  well  to  the  pressure  perturbation.  As  shewn 
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above,  <f  satisfies  the  equation  L dj  * 0 where  the  operator  L is  defined  by 
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and  where  wo  have  introduced  the  undisturbed  relative  Mach  number  seen  in 
blade-fixed  coordinates,  M,-  = 'Jit/a . If  wc  eliminate  F^in  terms  of  the 
operator  /.  in  Green's  theorem,  we  get 
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The  second  volume  integral  on  the  right-hand  side  can  be  rewritten 


as 


[(cm;  - cm;  £!* w . ... 

J K * as*  * js*'  j * 2s 


= f Ma  4-  (GM„  4-  ~ 

a s 


3 u 


ft  - it) 


dV 


CA/V i 

'•*'  W i 


In  this  form,  the  volume  integral  can  be  converted  to  a surface  integral  by 
applying  the  divergence  theorem  to  the  product  of  a scalar  function  $ and 
a vector  A . 
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The  volume  integrals  in  the  divergence  theorem  can  be  identified  with  the 
volume  integral  on  the  right-hand  side  of  Eq.  (60)  if 
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where  the  cylindrical-coordinate  form  of  Eq.  (63)  readily  shows  that  V-A  vanishes 
Using  Lqs.  (61)- (63)  in  Eq.  (60),  Green's  theorem  takes  the  form 
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where  41  and  R are  given  by  Fqs.  (62)  and  (03). 


This  relationship  can  be  used  to  express  the  velocity  potential  in 
terms  of  surface  integrals  which  bound  the  region  of  solution  by  the  usual 
Green's  function  technique.  The  scalar  function  G in  Eq.  (64)  is  chosen 
to  be  the  Green's  function  having  the  property  that  it  is  the  solution  to 
the  governing  differential  equation  for  a point  disturbance.  If  r and  r0 
denote  the  observation  point  and  the  source  point,  G(r,rc  ) satisfies  the 
equat ion 

L G ( r , rB  ) * 6 ( r - r0  ) IbS) 

where  d’  denotes  the  Dirac  delta  function.  If  the  integrations  in  Lq . (64) 
are  taken  to  be  over  the  source  coordinates,  and  within  the  region  of  inteiest 
Le  cf  ( ) *0,  then  the  following  integral  expression  is  obtained  for  the  ve- 

locity potential . 


The  term  tf(r)  is  the  result  of  integrating  cfira ) La  & (r,  ra  j over  the  volume, 
which  requires  that  the  Green's  function  determined  from  Eq.  (65)  also  has 
the  property 

Lo0  (?,?,)*  6 (ra  - 7)  (o7 

An  alternative  approach  to  the  integration  over  the  region  containing  the 
singular  point,  which  may  be  more  rigorous,  is  to  exclude  this  point  from 
the  volume  integral  by  surrounding  it  with  a vanishingly  small  surface. 

Then  the  functions  C?  and  G are  continuous  and  differentiable  throughout  the 
region  of  interest,  and  the  volume  integral  in  Eq.  (64)  vanishes.  However, 
there  would  then  be  a contribution  from  the  integration  over  the  surface 
enclosing  the  point  f ■=  <r9  . It  has  been  demonstrated,  using  the  Green's 

function  determined  in  the  subsequent  section,  that  the  integral  over  such 
a surface  yi  Ids  the  same  result  as  the  volume  integral  over  the  delta 
function.  Throughout  the  derivations  which  follow,  the  generalized  function 
approach  is  used  to  handle  point  singularities  because  the  presentations 
are  thereby  shortened  considerably. 


In  the  present  application,  the  surfaces  over  which  the  integrations 
must  be  done  in  Cq.  f66I  include  the  blade  surfaces,  the  duct  walls,  and  the 
surtaxes  normal  to  the  duct  axis  at  large  distances  upstream  and  downstream 
of  the  blade  row.  The  evaluation  of  the  surface  integrals  is  simplified 
considerably  through  the  use  of  a Green's  function  which  satisfies  the  sane 
boundary  conditions  at  the  duct  walls  as  the  velocity  potential.  Then  the 
surface  integrals  over  the  duct  walls  vanish  identically  and,  in  a direct 
parallel  with  isolated  wing  theory,  the  velocity  potential  for  a rotor  can 
be  found  by  superimposing  the  solutions  for  singularities  which  are  distributed 
only  over  the  blade  surfaces.  The  Green's  function  &(r\ra)  can  be  identified 

•*  r 

with  a mass  source,  and  then  ~~  corresponds  to  a fluid  doublet.  In  the 
linearized  analysis,  where  thickness  and  loading  effects  are  separable,  the 
source  solutions  can  be  used  to  repiesent  the  flow  about  a non!,  it  ting  rotor 
and  the  doublet  solutions  employed  for  a lifting  rotor.  However,  it  is  more 
convenient  to  treat  the  loading  case  in  terms  of  the  perturbation  pressure 
because  the  blade  ocundar}  conditicr.s  arc  expressed  m a simpler  form  and 
further,  integration  over  the  blade  wakes  is  avoided.  Since  the  perturbation 
pressure  satisfies  the  same  equation  as  the  velocity  potential,  the  formal 
solution  given  in  Lq.  (Go)  also  applies  to  the  pressure  field.  Here  the 
Green's  function  is  interpreted  a pressure  monopole  and  i*  , in  turn,  is 
differentiated  to  find  the  disturbance  field  of  a pressure  dipole.  The  pres- 
sure field  of  a lifting  rotor  can  then  be  found  bv  superposition  cf  the 
pressure  dipole  solution. 

In  the  next  section,  the  solutions  for  a point  source  of  mass  and 
. piessure  dipole  are  found  for  the  ducted  geometry. 


B.  MAKS  SOURCE  AND  PRESSUKL  DIPOLE  SOLUTIONS 

The  governing  differential  equation  ioi  the  velocity  potential  given 
in  Eq  (56)  cm  be  expressed  in  cylindrical  coordinates  by  using  Lq.  (Sla). 
Then,  the  velocity  potential  due  :o  a source  located  at  r„  , 6a  , in  a 
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rotating  reference  frame  and  having  a mass  addition  rate  of  satisfies 
the  equation 
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where  M is  the  Mach  number  based  on  the  undisturbed  axial  velocity,  U/a„ 
and  /S a . In  this  torn,  the  homogeneous  equation  is  separable  and, 

4 1 

as  originally  shown  by  McCune,  it  possesses  the  following  eigenfunction 
solutions  when  the  boundary  condition  of  no  flow  through  the  walls  is  enforced. 
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The  quantity  is  a normalised  combination  of  the  Bessel  and  Neuman 
of  order  n as  described  in  Apperjix  A and  cr  * <~/rT.  Kn^  is  the  -k'”  eigen- 
value of  the  equations  which  n suit  from  the  boundary  condition  that  S ^ / i n 
vanish  at  the  duct  walls. 
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where  k * r ,.>r  . The  quant; ry is  defined  by 
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with  Mr  - u>rx  / a.m  . Transformed  to  duct-fixed  coordinates,  these 
homogeneous  solutions  represent  the  duct  acoustic  modes,  Eq.  (71)  contains 
the  so-called  cutoff  condition  for  the  propagation  of  these  modes.  When 
i~lj  > £>  (^,i»  't : , A „ ^ becomes  imaginary  and  the  solutions  in  tq.  (b9)  cor- 
respond to  propagating  waves.  The  cutoff  condition  can  be  stated  approximately 
as  requiring  that  the  relative  Mach  number  at  the  tip  radius  must  be  super- 
sonic for  the  source  to  excite  propagating  modes.  Here  we  restrict  attention 
to  the  subsonic  case  where  the  modes  decay  with  distance  from  the  source. 
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In  order  to  solve  Eq.  (68)  the  form  of  Eq.  (69)  suggests  we  assume 
an  expansion  for  S °f  the  form 

o.-  » m(d-et) 

<?a  - l l (72) 
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where  will  denote  the  solution  for  a mass  source.  Here  = i for  u >o 
> 

but  *,*0  for  in  order  to  include  the  nontrivial  zero  eigenvalue  re- 

quired to  make  the  zeroth-order  Bessel  functions  a complete  set.  In  addition 
we  introduce  the  expansions  of  the  delta  functions  in  terms  of  the  azimuthal 
and  radial  eigenfunctions . 
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Substituting  these  expansions  into  Eq.  (68),  using  the  differential  equation 
satisfied  by  the  radial  eigenfunctions, 
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and  making  use  oi  the  orthogonality  properties  of  the  azimuthal  and  radial 
eigenfunctions  leads  to 
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where  R„i_(<r,)  is  introduced  as  a shortened  notation  for  O • 

The  solution  of  this  equation  for  ^*(2)  can  be  found  using  Fourier  transform 
techniques.  Kith  the  following  definition  of  a transform  pair 
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taking  the  transform  of  Eq.  (75)  yields 
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Solving  for  |ni 

and  taking  the 

inverse  transform. 
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The  integral  in 

Eq.  (To)  can  be 

evaluated  by  residue 
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denominator  being 
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For  2 » the  contour  is  closed  in  the  upper  half-plane  and  encloses  the 
pole  at  t't  , for  £ ' i,  the  contour  is  closed  in  the  lower  half-plane  and 
encloses  the  pole  at  if  - f,,*  . This  procedure  ensures  that  the  solution 
decays  rather  than  diverging  for  1 - t (or  for  supersonic  tip  speeds,  cor- 
respo1  's  to  outward  moving  waves). 


The  case  n-o,-k'C  deserves  special  attention.  For  those  values  of 
n and  A the  integrand  in  Eq.  (78)  has  a second-order  pole  at  f « o . The 
contribution  from  this  pole  is  included  in  the  contour  which  encloses  the 
upper  half-plane  (corresponding  to  2 > i.  ) and  excluded  from  the  contour  en- 
closing the  lower  half-plane  (corresponding  tc  2 <■  2.0  ) . This  choice  is  made 
on  the  grounds  that  there  can  be  no  steady  perturbation  at  upstream  infinity. 
The  integrated  results  for  cPn^  are 
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where  H (£-,%)  is  the  Heaviside  step  function  and  £0)0 
hit h this  result  for  <?„$  and  the  definition  of  in  6q.  (71),  Eq.  (72)  be- 
comes 
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The  first  tern  in  the  soi  rce  solution  has  been  omitced  in  previous 
41  44-40 

treatments.  ’ Except  for  the  presence  of  this  term,  the  above  result  can 

be  integrated  in  the  radial  direction  to  recover  the  line  source  solution 
which  McCune  used  to  solve  the  rotor  thickness  problem.  The  omission  of  this 
term  in  the  mass  source  (or  pressure  monopole)  solution  also  affects  the  fluid 
doublet  (or  pressure  dipole)  solution.  The  implications  that  omitting  this 
term  from  the  source  anc  dipole  solutions  would  have  to  the  results  for  the 
rotor  thickness  and  loading  problems  are  elaborated  upon  below  where  those 
solutions  are  presented. 


As  a result  of  discovering  the  omission  of  this  term  in  previous 
analyses,  several  checks  were  made  on  the  revised  source  solution.  The  fiist 
test  made-  on  the  solution  for  <?,  was  to  substitute  it  back  into  Eq.  (08) , 
and  to  verify  that  it  was  indeed  the  correct  solution.  In  addition,  the  pres- 
sure and  velocity  fields  associated  with  the  mass  source  solution  were  obtained 
from  t he  velocity  potential,  and  it  lias  been  verified  that  the  solution  displays 
the  properties  of  a mass  source.  The  flow  field  quantities  required  for  this 
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demonstration  arc  the  perturbation  pressure,  -p}  , and  the  axial  and  tangential 
velocity  components  (u.)s  and  (ir0}t  • The  velocity  components  are  just  the 
appropriate  derivatives  of  the  velocity  potential.  ihe  pressure  pertuibat Ion 
can  be  expressed  in  terns  of  tr,  and  by  using  Eqs.  (51a)  and  (5J)  together 
with  the  definition  of  c?  . 
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The  resulting  expressions  for  the  flow  Held  variables  are 
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where  o^»i  (H-2,1  is  ■•1  for  2 > and  -1  for  £ < £t  . 


(85c) 


Now  consider  a control  surface  which  encloses  the  point  source  and 
is  bounded  by  planes  normal  to  the  duct  axis  at  '£  * Z^t  t and  by  the  duct 
walls.  Next,  the  linearized  perturbation  in  the.  mass  flow  rate  through  the 
surfaces  bounded  by  this  control  volume  is  computer.  Since-  the  source  solution 
satisfies  the  boundary  conditions  of  no  flow  through  the  duct  wall,  the  only 
mass  flux  is  through  the  surfaces  normal  to  the  duct  axis.  Letting  ■'w 
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represent  the  first  order  terms  in  the  axial  component  of  the  mass  flux, 
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The  density  perturbation  can  be  expressed  in  terms,  of  and  lr£  using  F.qs.  (34) 
and  (82).  As  a result  becomes 
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If  7n  Q is  the  rate  of  introduction  of  mass  by  the  source  then  we  must  have 
for  this  control  volume 
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Taking  the  limit  of  this  expression  as  6 — o and  l — • ► , and  noting  that  v6 
is  continuous  at  2„  , then 
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Using  the  series  expansions  l'or  the  delta  functions  given  in  Eq.  (73),  this 
last,  equation  can  be  written 
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and  hence  the  mass  balance  in  Eq.  (87)  is  satisfied. 


(89) 


In  addition  to  the  demonstration  of  mass  conservation  for  the  above 
control  volume,  a similar  check  has  been  made  for  the  control  volume  obtained 
by  withdrawing  the  surfaces  normal  to  the  duct  axis  to  the  location  7.  = - » . 
It  also  has  been  shown  that  the  results  for  the  flew  field  of  a point  source 
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satisfy  the  axial  component  cf  the  linear  and  angular  momentum  balances  for 
this  latter  contr  >1  volume.  These  tests  on  the  source  solution  are  given  in 
Appendix  C together  with  comparable  tests  for  the  pressure  dipole  and  rotor 
flow  field  solutions. 


The  velocity  potential  of  a fluid  doublet  can  be  obtained  from  that 
for  a source  in  the  conventional  way.  and  the  flow  field  produced  by  rotor 
loading  could  be  found  by  superimposing  these  doublet  solutions.  However,  as 
discussed  in  the  previous  section,  it  is  more  convenient  to  work  with  the  per- 
turbation pressure  in  the  loading  problem.  Such  an  approach  is  analogous  to 
the  acceleration  potential  method  in  isolated  airfoil  theory  where  the  dis- 
turbance field  of  an  infinitesimally  thin,  lifting  surface  is  constructed  by 
the  superposition  of  pressure  dipoles. 


Since  the  rotor  pressure  field  satisfies  the  same  equation  as  the 

velocity  potential,  the  solution  for  a pressure  monopole  is  mathematically 

the  same  as  the  solution  for  a mass  source.  The  corresponding  dipole  solution 

is  found  by  differentiating  the  monopole  solution.  If  the  pressure  monopole 

is  to  be  used  as  the  Green's  function  in  the  formal  solution  for  the  pressure, 

then  the  required  orientation  of  the  dipoles  is  normal  to  the  blade  surfaces, 

or  in  the  linearized  analysis,  to  the  undisturbed  stream  direction.  For  a 

dipole  of  strength  D located  at  the  point  rt  and  oriented  in  the  positive 

direction  shown  in  Figure  20,  the  pressure  field,  -p  , is  given  bv 
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where  G-ir,  T\)  is  the  source  (or  monopole)  solution  given  in  F.q.  (81)  with 
unit  strength.  Performing  the  indicated  operations  in  Eq.  (90), 
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As  with  the  result  for  . this  solution  for  -p  has  been  substituted  into  the 
governing  differential  equation  to  verify  that  it  is  the  correct  solution. 
Forming  the  quantity  L-p  , where  L is  the  operator  defined  by  Eq.  (58),  yields 
the  result 
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when  the  series  expansions  in  Eq.  (73)  are  used  for  the  delta  functions. 


The  velocity  components  associated  with  the  dipole  field  are  obtained 
by  integrating  the  momentum  equations,  Eq.  (49),  along  the  undisturbed  stream- 
lines. Consistent  with  our  generalized  function  approach  in  treating  singular 
points,  a delta-function  body  force  is  included  in  the  momentum  equations. 

Then,  the  resulting  expressions  for  the  velocities  are  valid  everywhere  in  the 
duct,  including  points  which  lie  on  the  streamline  that  passes  through  the 
dipole  location.  Otherwise,  the  expressions  would  not  be  valid  in  a small 
region  enclosing  this  streamline. 


It  is  convenient  to  work  with  the  streamwise  and  normal  velocity 
components,  both  because  fewer  integrals  have  to  be  evaluated  and  because  the 
body  force  term  appears  only  in  the  normal  momentum  equation.  The  dipole 
exerts  a force  per  unit  volume  on  the  fluid,  F0  , which  is  in  the  negative  r\ 
direction  and  expressed  by 
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The  streamwise  velocity  component  is  simply  proportional  to  the  pressure  by 
Eq.  (53),  which  was  derived  by  integrating  the  streamwise  momentum  equation. 
Integrating  the  radial  and  normal  momentum  equations  along  the  undisturbed 
streamlines  yields 
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The  integrations  along  the  undisturbed  streamlines  can  be  done  by  expressing 
as  in  terms  of  d.  £ . Along  these  streamlines,  r and  the  helical  variable,  % 
which  is  defined  by 
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remain  constant,  and  the  arc  length  along  the  stx-eamlines  is 
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Making  this  change  in  the  integrals,  Eqs.  (94)  and  (95)  become 
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The  results  arc 
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These  expressions  can  be  substituted  into  those  fort  Wr)0  andCir^;^  to  complete 
the  derivation  of  the  velocity  components.  Also,  the  velocity  components 
(^rn)D  and  are  related  to  (ir$)D  and  ( Un ) D by 

Vq  ~ V5  A*,-*  y + Vn  ^ c 1 o 1 a) 

V - vs  * y/-  - tr,  ys  (101b) 


These  results  for  the  flow  field  of  a pressure  dipole  have  been 
shown  to  satisfy  the  mass  and  momentum  balances  given  in  Appendix  C.  The 
dipole  solution  possesses  the  appropriate  properties  that  it  does  not  intro- 
duce any  mass  into  the  flow,  and  exerts  a force  D on  the  fluid.  In  the  next 
two  sections,  the  source  and  dipole  solutions  in  Eqs . (81)  and  (91)  are  used 
as  the  Green's  functions  to  construct  the  flow  field  produced  by  a rotor  with 
distributed  thickness  and  loading. 


C.  FLOW  FIELD  OF  A NONLIFTING  ROTOR  (THICKNESS  PROBLEM) 

1 . Solution  for  the  Veloc ity  Potential 

In  the  previous  two  sections,  the  foundation  has  been  laid  to-  develop 
the  solution  for  the  thickness  contribution  to  the  flow  field  of  a rotor  in 
an  annular  duct.  The  source  solution  given  in  Eq.  (81)  can  be  used  as  the 
Green's  function,  &(.?,  r ) in  the  integral  representation  of  the  velocity 
potential  in  Eq.  (66).  The  solution  for  a source-  of  unit  strength  (Q  * t ) 
satisfies  Eq.  (65)  and  also  has  the  property  required  by  Eq.  (67).  The 
integrations  in  Eq.  (6t>)  must  be  done  over  the  following  surfaces:  (i)  the 

duct  walls,  (ii)  the  surfaces  normal  to  the  duct  axis  at  large  distance  from 
the  rotor,  and  (iii)  the  blade  surfaces.  In  describing  the  integrations 
over  each  of  these  subsurfaces,  I,  and  I,  will  refer  to  the  first  and  second 
surface  integrals  ir  Eq.  (66),  respectively. 
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(i)  Along  the  inner  and  outer  duct  walls  the  boundary  condition  on 
G (.■'„)  is  that  the  normal  derivative  vanish  corresponding  to  no 
flow  through  the  walls.  Since  the  Green's  function  we  have  found 
satisfies  the  same  boundary  condition,  I,  vanishes  for  this  surface. 
For  the  outer  wall  i>a  =-  er  while  for  the  inner  wall  va  - er  . 

From  Eq.  (65)  we  see  that  Vc  ■ R - o along  both  the  inner  and  outer 
walls,  and  so  1 ^ also  vanishes  at  the  duct  walls. 


(ii)  For  the  surfaces  normal  to  the  duct  axis  at  z and 

is-<?a  and  + respectively.  Over  these  two  surfaces  the  sum  of 
the  integrals  1 , and  Ix  become 
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Consider  the  integrand  at  ZQ—^t<n  first.  In  this  case  2 < Z„  and  so 

the  first  term  (^-<5,  * O ')  in  G ( P , ) is  zero.  The  remaining 

terms  in  <J(r,r.)  decay  exponential  ly . The  same  behavior  holds  for 

-r-r- and  -r—  ■ Since  d? r ) and  its  derivatives  must  be  bounded,  the 
d Ue  v £ & 

integi-arid  vanishes  as  £<,  — *•*<*>  . Next,  consider  the  integrand  evalu- 
ated at  £„—*•- ® . Here  2 > 2,  and  the  /!•<?,  ii  * O term  contributes 

2 r 2 /*  o r 

to  both  G and  ; the  remaining  terms  in  G , ~~  , and  decav 

it,  d 9. 

exponential  1)'.  lor  subsonic  flow,  dj  and  the  velocity  components 
obtained  from  its  derivatives  are  required  to  vanish  far  upstream  of 
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the  blade  row.  Note  that.  G diverges  linearly  as  x)  -*•  and  so  the 
velocity  field  must  fall  off  faster  than  this  in  order  for  the  inte- 
grand to  vanish.  We  shall  sec  that  the  velocity  tj.old  decays  expo- 
nentially upstream  of  the  rotor.  Hence,  the  integrands  van i. si.  for 
2„ -*-<*>  also,  and  the  integrals  2,  and  lz  make  no  contribution  along 
the  surfaces  normal  to  the  duct  axis.  For  supersonic  tip  speeds, 
acoustic  waves  propagate  away  trom  the  blade  row  and  these  surface 
integrals  should  be  re-examined. 


( i i i ) In  the  linearized  analysis,  the  blade  row  is  assumed  to  make  only  a 
slight  perturbation  of  the  iree-stream  flow.  Consiste  with  this 
assumption,  the  blade  surface  boundary  conditions  are  applied  along 
the  undisturbed  stream  direction.  In  this  approximation,  the  nor- 
mals to  the  upper  and  lower  blade  surfaces  are,  respectively, 


t)  m + r' 
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where  y,  - tan*' (tor,  / u)  . From  Eqs.  (63)  and  (T">3)  it  can  be 
seen  that  ^ vanishes  and  thus,  in  the  linearized  approximation, 
] contains  no  contribution  from  the  blade  surfaces. 


The  separation  of  the  rotor  flow  field  into  the  thickness  and  load- 
ing contributions  is  made  by  prescribing  that  there  be  no  pressure  difference 
across  the  blade  surfaces  in  the  thickness  case,  ox  that  the  blades  are  lo- 
cally unloaded.  If  the  pressure  is  continuous  across  the  blade  surface,  then 

J Q 

and  (f  are  also.  Hence,  because  of  the  opposite  signs  of  on  the 

av,  * 

upper  and  lower  surfaces,  only  the  part  of  the  integrand  in  I,  containing  G j 
contributes  to  the  integration  over  the  blade  surfaces. 

Except  for  this  contribution  from  integrating  f,  over  the  blade 
surfaces,  all  the  surface  integrals  in  F.q.  (66)  vanish.  Thus,  the  expression 
for  (f  has  been  reduced  to 

cP(r ) * [ G^ri  d 50  (104) 
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where  Se  denotes  the  surface  area  of  the  rotor  blades  projected  on  tne  un- 


disturbed stream  surface,  and  — ] represents  the  difference  in  i.ormal 

w n,  / 

velocity  across  each  blade  surface.  This  expression  is  the  same  as  the 
familiar  result  in  isolated  airfoil  theory  that  the  effects  of  wing  thickness 
can  be  represented  by  the  superposition  of  sources  whose  strength  is  equal 
to  the  discontinuity  in  vn  at  each  point. 


The  i in ear iced  form  of  the  blade  boundary  conditions  is 
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= U, 


as 
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where  and  are  the  distances  to  the  upper  and  lower  surfaces,  measured 
noimal  to  the  undisturbed  stream  direction,  S • The  quantities  <fyu_  and  9t\^ 
can  be  expressed  in  terms  of  a blade  rni  ,kness  and  a blade  incidence  plus 
camber  line  in  the  conventional  way  as  illustrated  in  Figure  21.  However, 
for  a rotor  the  blade  thickness  and  camber  arc  no'-  independent  because,  as 
noted  earlier,  in  order  for  the  blades  to  bo  unloaded,  they  must  be  cambered 
to  account  foi  loadings  which  would  otherwise  be  induced  by  blade  intcrferer.ee 
effects.  If  t(i,r)  represents  Vl^O.r)  - then. 


= <Ylc(s,r)  * -j  t (s,r) 
X(s>  r)  " 'ttc'S,  r?  - y t (5  , r) 


(10b) 


where  l!Mt  (S,  r)  is  the  camber  line,  which  consists  of  two  parts,  a thickness 
part,  H4t  and  a loading  port,  . From  bqs.  (105)  and  (10b;,  the  dis- 
continuity in  the  noimal  velocity  across  r.kc  blade  surface  is  then  related 
to  the  thickness  distribution  by 


A K, 


it 
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For  a rotor  with  5 equally  spaced  blades,  the  blades  lie  on  the 
surfaces  % • J.jn/3,  j'o,  3 •/  where  if  is  the  helical  variable  defined 
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by  Lq  (9o).  It  i convenient  to  do  trie  integration  over  these  surfaces  by 
projecting  them  on  the  r0  , plane  so  that 
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With  the  blade  leading  edges  located  at  z„  = o , the  velocity  potential  is 
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While  not  crucial  to  the  ensuing  analysis . it  is  convenient  to  assume  that 
the  axial  projection  of  the  blade  chord  is  a constant,  C , and  that  the 
radial  and  axial  variations  of  At/,,  (and  t ) are  factorable  so  that 
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With  these  assumptions,  the  expression  for  the  velocity  potential  becomes 
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llic  summation  over  the  number  ot*  blades  can  be  done  by  using  the  identity 
«"  fit?  ) fc  r,/ 

B r.  - m Q I (112) 


y c^i-; 

L « 


where  mis  an  integer.  The  result  for  the  velocity  potential  becomes 
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where  we  have  introduced  the  notation 
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2 . Expressions  for  tlic  Velocity  Components  and  Pressure  Perturbation 

Having  this  solution  for  the  velocity  potential,  the  other  flow 
field  variables  can  be  found  by  taking  the  appropriate  derivatives.  The  re- 
sults for  the  velocity  components  are  given  in  terms  of  v*  , v„  , and  vr 
which  are  related  to  c?  by 
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The  pressure  is  simply  proportional  to  vs  by  Tq.  C53).  The  resulting  expres- 
sions for  the  velocity  components  are 
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where  the  prime  on  uenotes  differentiation  wich  respect  to  the  argument 

and 
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also  has  been  introduced. 


These  results  for  the  flow  field  of  a nonlifting  rotor  are  shown 
to  satisfy  the  mass  and  momentum  balances  in  Appendix  C,  For  subsonic  flow, 
there  should  be  no  net  foice  or  work  done  by  the  rotor.  The  derivations 
presented  in  the  appendix  have  verified  that  this  is  so.  Another  check  which 
has  been  made  on  these  expressions  is  that  they  display  the  correct  behavior 
at  the  blade  surfaces 

3 . Behavior  of  the  Velocity  Components  at  the  blade  Surfaces 

The  expressions  given  for  the  disturbance  velocity  field  contain 
doubly  infinite  series  expansions  in  the  duct  eigenfunctions.  In  order  to 
demonstrate  that  the  velocity  components  display  the  correct  behavior  at  the 
blade  surfaces,  the  convergence  properties  of  these  series  need  to  be  con- 
sidered. Terms  in  the  cm  summat ions  which  are  of  order  (mS)  arc  expected 
to  lead  to  divergent  series,  or  series  which  do  not  converge  uniformly. 

The  manipulations  performed  below  are  aimed  at  identifying  those  terms  wh:rh 

- I 

are  0 Cm  5 ) and  hence  are  potentially  divergent.  Furthermore,  within  th-  ■ 
group  of  terms,  we  wish  to  isolate  a series  for  which  the  ^ sununation  can  be 
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done  analytically  and  for  which  the m summation  does  not  converge  uniformly 
foi  all  values  of  the  t,  coordinate.  It  can  be  anticipated  that  such  a series 
produces  the  discontinuities  in  the  surface  quantities  which  occur  as  changes 
in  £ are  made  which  correspond  to  crossing  a blade  surface.  For  the  remainder 
of  this  section  our  attention  is  confined  to  points  within  the  blade  row, 
i . e.  , for  O 2.  iz  Co_  . 

First  we  wish  to  show  that  vs  is  continuous  across  the  blade  sur- 
faceo.  This  demonsti at  ion  also  provides  a check  on  the  assumption  that  the 
blades  are  locally  unloaded  by  virtue  of  the  relationship  between  -p  and  Vg. 

The  first  step  in  examining  the  value  of  irs  at  the  blade  surfaces 
is  to  carry  out  an  integration  by  parts  in  the  integral  ever  . After  the 
first  integration  by  parts  the  expression  for  is 
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Were  it  not  for  the  A -dependence  of  the  factors  other  than  % rr,  el  > 
the  A -summations  in  this  expression  would  be  just  the  Fourier-Dessel  series 
expansion  of  the  form  (see  Appendix  A) 
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Since  -f  (r)  miut  be  independent  of  the  order  of  the  Bessel  functions  in  tins 

expansion,  the  product  need  not  be  considered  in  determir int; 

the  m -dependence  of  each  term  in  the  double  scries.  The  coefficients  in 

the  series  containing  and  ^(c^l  are  inversely  proportional  tc 

which  is  O(rnB) , but  the  exponential  factors  in  these  terms  prevent  the 

series  from  diverging,  except  as  the  leading  or  trailing  edge  is  approached. 

There,  unless  the  slope  of  the  thickness  profile  vanishes,  the  series  diverge. 

The  divergence  of  these  scries  produces  the  singularities  in  the  pressure 

which  typically  occur  at  the  leading  and  trailing  edges  of  subsonic  airfoils. 

The  divergence  of  the  surface  pressure  at  these  points  is  evident  ii  the 

4 1 

original  results  of  McCune,  as  well  as  in  the  surface  pressure  results  ac- 
companying the  thickness-induced  camber  lines  presented  in  Reference  46. 


Each  integration  by  parts  over  introduces  another  factor  of 
(<r>8)  in  each  term  of  the  double  series.  A second  integration  by  parts  in 
the  u5  expression  produces  the  following  results. 
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The  second  integration  by  parts  yields  integrated  terms  which  are 
proportional  to  (j'a,)  , evaluated  at  » o , z , and  Cn  . These  terms 

are  0 (r^a)'2  and  so  the  corresponding  series  arc  uniformly  convergent.  Hence, 
away  from  the  leading  and  trailing  edges,  all  the  series  in  zrs  are  uniformly 
convergent  for  all  values  of  2f  , including  those  at  the  blade  surfaces.  It 
can  therefore  be  concluded  that  lts  and  p are  continuous  across  these  surfaces. 

ipr 

The  same  is  not  true  ot  un  , however,  which  must  be  discontinuous  at  ; - — 
by  an  amount  dictated  by  the  boundary  condition,  Eq.  (107). 


We  have  already  evaluated  the  ir.  contribution  to  the  expression 
given  for  i/n  in  Eq.  (115b),  and  shown  it  to  be  continuous.  Any  discontinuity 
in  vn  must  then  come  from  the  contribution.  After  performing  a single  in- 
tegration bv  parts  on  the  z integral,  the  expression  for  4v  becomes 

* * 5 


2 A , 


V. 


r.  ?. 


-4  re 


V l 


Z Z 

* '*>  * * 1 


( 4 niB) 


Zqd> 


<}  (<0  £ 


-A-rnSllU  ,0)2 


(i/J  [ {37  r A™57 


A „ 


U 


J / £ Z1 


uj  A 


m0 


* Z dl° 


ot?. 


(120) 


O ' O 


where  the  prime  on  the  summation  indicates  the  m * o term  is  excluded.  The 
leading  and  trail ing-edge  terms  behave  as  discussed  above,  while  the  re- 
maining integral  is  0 (nn  &)  as  in  the  V't  expression.  However,  the  first 
of  the  integrated  terms  is  0(rnB)~'  and  does  not  contain  the  exponential 
factors  which  appear  in  the  leading  and  trailing  edge  terms.  This  term, 
which  may  diverge,  or  at  best  not  converge  uniformly,  will  be  found  to  produce 
a discontinuity  in  vn  ■ 
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Next  we  consider  the  behavior  of  the  first  term  in  ) by  definin 
the  double  series  Si V/)  to  which  th:s  term  is  proportional. 
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Substituting  Eq. 
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If  the 
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The  4 summation  in  the  first  line  of  F;q.  (123)  is  just  the  Fourier- 
Eessel  expansion  for  the  bracketed  function  in  terms  of  the  radial  eigenfunctions. 
Performing  this  summation  and  using  the  equation  satisfied  by  the  radial 
eigenfunctions,  Eq.  (74),  to  substitute  in  the  second  integral,  S(r,  %)  can 
be  written 
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Now,  the  eigenvalues  are  all  greater  than  mg,  being  O(r^B)  for  large 

rr.  S , and  the  remaining  integral  can  be  integrated  by  parts. 
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The  integrated  term  vanishes  because  each  radial  eigenfunction  identically 
satisfies  the  boundary  conditions  at  the  duct  walls.  Another  integration  by 
parts  could  be  done  but  it  does  not  appear  worthwhile.  The  important  point 
to  make  is  that  each  m term  of  the  double  series  in  Eq.  (124)  is  at  most 

of  order  (mfl)’1  . Here  again  this  series  is  then  a regular  series,  and  the 
first  series  in  Eq.  (124)  is  the  only  remaining  contribution  to  S(r,n;) 
which  could  contain  a discontinuity.  The  single  summation  on  the  right-hand 
side  cf  Eq.  (124)  is  proportional  to  the  Fourier  series  expansion  of  the 
generalised  function,  , which  is  defined  by 
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This  function,  illustrated  in  Figure  22  was  first  used  by  Ueissncr  ‘ in  his 
representation  of  a propeller  wake.  It  can  be  shown  that  the  Fourier  series 
representation  of  is 
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Kith  those  results,  the  expression  for  Sir, % ) becomes 
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As  can  be  seen  from  Figure  22,  the  generalized  function  £ has  a jump  of  [ 3 

as  ^ crosses  a blade  location  moving  in  the  direction  of  increasing  if  . The 

contribution  to  j~\ , from  the  term  containing  S(r,f)  is  simply  j(e)  5(r,tf) 

Combining  this  result  with  Eq.  (115b)  to  find  the  corresponding  term  in 
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and  trie  definition  of  the  functions  f ?ud  ^ in  terms  of  £ vn  given  in  Eq.  (110) 
is  retrieved.  Since  the  other  contributions  to  i rn  are  continuous  for  all 
values  of  < , and  for  values  cf  I away  from  the  leading  and  trailing  edges, 
we  have  succeeded  in  showing  analytically  that  our  result  for  the  velocity 
field  does  in  fact  contain  the  correct  jump  in  normal  velocity  across  the 
blade  surfaces. 


If  the  results  for  v$  and  are  collected  to  find  the  whole  ex- 


pression for  vn  according  to  Eq.  (115b),  the  answer  is 
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The  first  term  in  ir„  contains  only  the  symmetric  discontinuity  and 
hence  makes  r.o  contribution  to  the  continuous  part  of  vn  at  the  blade  surfaces. 
When  evaluated  at  a blade  surface  the  second  term  in  vn  vanishes  and  the  re- 
maining terms  represent  the  continuous  part  of  vn  , or  the  thickness-induced 
camber  lines.  The  third  term  in  the  above  expression,  and  the  first  term  in  I'.o. 
(119),  come  from  the  n •>  O , -k  - r term  omitted  from  previous  treatments  of 
the  thickness  problem.  The  camber  line  calculations  of  Ref.  46  and  the  surface 
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pressure  calculations  of  Reference  H should  he  cortceted  for  the  presence  of 


these  Terms.  Since  Lh 


dur,  sj 


• t j ; . , i . , i 3ti  r\  2 ) 

- — , which  may  be  written  as  u — ~-£ — > is  pro- 


portional to  ?. ) , these  contributions  to  both  i-r,  and  f>  are  proportional  to 
the  chordwise  variation  in  the  local  thickness  distribution.  Furthermore, 
these  terms  make  no  contribution  to  the  disturbance  field  upstream  or  downstream 
of  the  blade  row,  and  therefore , the  acoustic  calculations  of  References  44  anu 
4;  are  unaffected  bv  their  inclusion. 


The  above  manipulations  of  the  expressions  for  the  perturbation  ve 
loc.ities  apply  for  points  within  the  confines  of  the  rotor  row,  i.e.,  for 
0 £ ?.  £ C . If  the  integrations  by  parts  are  done  for  field  points  upstream 
and  downstream  of  the  blade  row  the  singularities  encountered  at  the  loading 
and  trailing  edges  remain  but  the  velocities  everywhere  else  in  the  flow  are 
found  to  be  continuous. 


The  final  expressions  presented  for  v%  and  va  in  Eqs.  (11R)  and 
(130)  have  been  snown  to  have  the  correct  behavior  at  the  blade  surfaces. 
Moreover,  in  the  course  of  demonstrating  this,  these  results  have  been  put  in 
a form  which  siould  facilitate  the  computation  of  the  surface  pressure  and  the 
tnickness- induced  camber  lines  for  nonlifting  rotor  blades. 


D.  FLOW  1- 1 ELD  OF  A LIFTING  ROTOR  (LOADING  PROBLEM) 

1 . Solution  for  the  Perturbation  Pressure 

The  determination  of  the  pressure  field  produced  by  a lifting  rotor 
closely  parallels  the  solution  procedure  for  the  velocity  potential  in  the 
thickness  problem.  The  Green's  function  is  the  same,  though  now  its  interpreta- 
tion is  in  terms  of  a pressure  inonopolc.  The  formal  solution  for  p is  identical 
to  that  given  in  Eq.  (66)  for  q , except  for  the  replacement  of  if  with  p . 

In  addition,  both  of  the  integrations  over  the  duct  wails  and  the  surfaces  at 
I os  vanish  as  before,  as  does  the  integration  of  i)  ■ ( fi  $ ) over  the  blade 

b" 


s 


surfaces.  The  remaining  integral  over  the  blade  surfaces  distinguishes  the 

a t 

lifting  case  iron  the  thickness  case.  Here  the  magnitude  of  is  the  same 

on  both  sides  of  the  biades  because  vn  must  be  continuous.  Since  the  normals 

to  the  upper  and  lower  surfaces  lie  in  opposite  directions  and  G ‘ r , ) is 

continuous  across  the  blades,  the  integral  of  G K~  over  the  blade  sur- 

0 

faces  vanishes.  Thus  the  integral  representation  of  the  solution  for  the 
pressure  field  reduces  to 


-p(r)  = / Ap(re 


10. 


where  SB  is  the  projection  of  the  blade  surface  on  the  undisturbed  stream  di- 
rection and 

ciuo  = p (r  £ £ ~ — f-  ) - *o  ( r Z £ = - 77 ^ ) 

■ B ) ' \ ' ° ' J'  3 J (133) 

The  sign  of  the  integrand  in  Eq.  (132)  has  been  reversed  from  that  in  Eq.  (6b) 

so  that  .1;.  is  defined  as  a positive  n umbel  . 

If  w?  substitute  the  pressure  dipole  solution  (Eq.  (91))  with  unit 

strength  for  i— - and,  exactly  as  was  done  ford  it.  the  thickness  problem, 
a rm 

carry  out  the  summation  over  the  number  of  blades,  we  get  the  following  result 
for  the  pressure  field  of  the  rotor 


mfilU 
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where  we  have  introduced 


jA  I 
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and  the  quantities  -A-  and  (r„  . Z . ) are  del  ined  in  Eqs . (114) 

and  (ir). 


The  first  property  of  this  solution  to  be  examined  is  the  limtiing 
pressure  rise  between  points  far  upstream  and  far  downstream  of  the  blade  row. 
Since  for  subsonic  flow  the  solution  decays  as  Z -*•  - « , the  limiting  value 
of  the  static  pressure  rise  is  simply  the  limiting  value  of  p for  £ — • * » . 

All  terms  except  the  first  decay  and  so 

■^B  r C‘~ 

f -**)  * : — / Ap(i,)dre  (.156) 

P.rt/3  U Jo 

This  result  agrees  with  that  obtained  from  the  vortex  theory  of  Okurounmu  and 
4 ■* 

McCune.  ' In  the  present  formulation  its  origin  is  in  the  n -o,  l.j  term  in 
the  dipole  solution.  When  we  initially  considered  the  pressure  dipole  formula 
tion  of  the  loading  problem,  this  tern  was  omitted,  and  the  erroneous  result 
of  no  limiting  pressure  rise  was  found.  After  tracing  the  difficulty  to  a 
missing  term  in  the  source  or  monopole  solution,  that  solution  was  revised 
as  discussed  in  Section  IV- B.  Mien  the  coi responding  correction  was  made  in 
the  dipole  solution,  the  result  for  the  pressure  rise  given  in  Eq.  (136)  was 
obtained . 


N3mba  reports  a limiting  static  pressure  rise  even  though  his  mono- 
pole  and  dipole  expressions  do  not  contain  the  n*o  , -rr  « o terms  found  here. 
However,  his  result  for  the  pressure  rise  differs  from  Eq.  (136)  in  two  respects. 
First  of  all,  he  finds  a nonvanishing  pressure  perturbation  for  Z — - ■»  , one 
half  of  the  static  pressure  rise  occurring  between  upstream  infinity  and  the 
rotor  with  the  other  half  occurring  downstream.  Also,  his  result  contains 
a "scale  factor"  which  introduces  an  additional  radial  dependence  into  the 
integrations  over  radius  in  Eqs.  (134)  to  (136).  Based  on  the  tests  we  have 
made  on  the  singularity  solutions,  and  those  described  below  for  the  flow 
field  of  a lifting  rotor,  we  have  concluded  that  this  factor  should  not  be 
present . 
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2 . Expressions  for  the  Velocity  Components 


In  the  pressure  dipole  representation  of  the  blade  row,  the  velocity 
field  must  be  found  by  integrating  the  momentum  equations  along  the  undisturbed 
stream  direction.  This  situation  in  the  dipole  approach  somewhat  compensates 
for  avoiding  the  distribution  of  singularities  over  the  blade  wakes,  as  would 
be  necessary  if  fluid  doublets  were  used  to  find  the  velocity  potential.  As 
for  the  thickness  case  the  velocity  field  will  be  found  in  terms  of  the  com- 
ponents trs  , trn  , and  vr  . 


for  subsonic  flow  the  streajnwi.se  velocity  perturbation,  irs  , is 

found  from  Eq.  (53).  If  y~-  is  expressed  in  terms  of  the  derivatives 

and  -ry)  then  integration  of  the  normal  momentum  equation,  Eq.  (52b), 
o C,  / i , r 

along  the  undisturbed  streamlines  yields 
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where  £ is  a dummy  variable  of  integration.  Note  that  c and  are  held 
constant  ir.  this  integration.  In  keeping  with  the  generalized  function  approach 
to  the  singularities,  a body-force  term  representing  the  blade  forces,  Fp  , 
has  been  included.  Since  this  force  acts  normal  to  the  undisturbed  stream 
direction  in  the  linearized  theory,  it  appears  only  in  the  zr„  expression. 

Before  proceeding  further,  this  contribution  to  wn  is  evaluated.  Each  blade 
contributes  a force  per  unit  volume  of 


- n 


A i t>  <r  S ) cf  ! n - xi 
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where  rt  is  the  unit  vector  normal  to  the  undisturbed  stream  surface  (see 
Figure  20),  <1  is  a coordinate  in  this  direction,  and  n is  the  value  of  n at 
the  -j  th  blade.  The  total  force  on  the  fluid  due  to  the  ^ blade,  , is  then 
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Returning  to  the  term  in  question,  the  coordinate  r is  expressed  in  terms  of 
K and  r by 


A m 


r K 


and  so  we  must  have 
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Substituting  these  results,  the  term  in  vn  containing  the  body-force  repre- 
sentation of  the  blade  forces  is 


P« 


Fe  d 


L V U / J e» 


J 

— CB 


Upr 


T 

J "O 


(f 


where  we  have  accounted  for  the  fact  that  the  blades  are  located  on  the  sur- 
faces < = ^4-  , between  the  axial  stations  c-o  and  I * . As  we  shall 
o 

see,  this  term  will  be  cancelled  by  other  terms  in  . For  now  the  remaining 
terms  in  vn  will  be  designated  vn  . If  the  blade  forces  were  treated  as 
surface  forces,  then  only  ir'n  would  be  present  at  this  point;  the  delta  func- 
tion terms  still  present  would  be  excluded  from  the  value  of  on  the  blade 
surface  on  the  grounds  that  the  singularities  are  within  this  surface. 


The  expression  for  vr  which  results  from  integrating  the  radial 


momentum  equation,  Eq.  (49a),  is 

, r z 
i J 
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o r 


X,  Z')  d ? ' 


(142) 


Eqs.  (137)  and  (142)  show  that  the  calculation  cf  the  velocity  components  re- 
quires the  integration  of  p along  the  stream  direction.  The  required  integral 
is  denoted  by  l (Z)  where 
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Substituting  for  ^ from  Eq.  (134)  and  doing  the  £ integration  yields 
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Now  the  results  for  f>  and  1(e)  can  be  combined  according  to  Eq.  (137) 
in  order  to  obtain  vn  ■ Since  the  explicit  representation  of  the  blade  force 
is  omittea,  the  remaining  terms  are  denoted  Vn  • 

(£r.)  ru 
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The  expression  for  isr  , obtained  by  combining  Eqs.  (142)  and  (144) 
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Particular  note  should  be  paid  to  those  terms  in  vn  and  vr  which  do 
not  decay  downstream  of  the  rotor.  These  terms,  which  are  preset,  within  the 
blade  row  and  downstream  of  it,  represent  the  contribution  of  the  trailing 
vortex  wakes  to  the  velocity  field.  The  flow  field  produced  by  these  wakes 
has  a helical  pattern  and, as  a consequence  of  the  linearization,  the  wakes 
coincide  with  the  undisturbed  stream  surfaces  on  which  the  blades  lie.  There 
are  no  wake  terms  in  ? or  v}  , which  should  be  continuous  across  these  sur- 
faces . 


The  results  given  in  this  section  for  the  lifting-rotor  flow  field 
satisfy  the  global  mass  and  momentum  balances  presented  in  Appendix  B.  In  the 
next  section  it  is  demonstrated  that  the  velocity  perturbations  found  for  the 
lifting  rotor  also  exhibit  the  proper  behavior  at  the  blade  surfaces. 

3 . Behavior  of  Velocity  Components  at  the  Blade  Surfaces  and 
across  the  Trailing  Vortex  Wakes 

The  properties  of  the  above  solution  for  the  velocity  field  of  a 
lifting  rotor  can  be  examined  in  much  the  same  way  as  was  done  in  the  previous 
section  for  a nonlifting  rotor,  although  the  loading  case  is  somewhat  more 
complicated  because  of  the  presence  of  the  trailing  vortex  wakes.  Again  the 
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straamwise,  normal,  and  radial  components,  irs  , ^•/I  , and  vr  , are  considered. 
First  the  behavior  of  vi  is  treated,  followed  by  y,  . The  examination  of  va 
will  be  done  last  because  it  provides  a natural  transition  to  the  discussion 
of  the  direct  lifting  surface  theory  in  the  next  section. 

a.  Streamwise  Velocity  Component 

Since  vs  is  proportional  to  -p  , it  must  contain  a discontinuity 
across  the  bl-*de  surfaces  which  is  in  the  same  proportion  to  ths  blade  loading, 
. Therefore,  it  is  equivalent  to  demonstrate  the  proper  behavior  of  v\ 
at  the  blade  surface  by  showing  that  the  solution  for  the  pressure  field 
contains  the  correct  discontinuity.  Again  this  is  done  by  ordering  the  series 

, r 

expansions  in  (mB)  and,  in  particular,  isolating  a series  fer  which  the  4 
summation  can  be  done  and  for  which  them  summation  yields  the  generalized 
function  £ . 


As  before,  the  first  operation  on  the  expression  for  jo  is  to  per- 
form an  integration  by  parts  on  the  20  integral.  The  result,  after  some 
rearrangement,  is 
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In  order  to  d scuss  the  behavior  of  the  various  terms  in  this  form,  it  is 
convenient  to  write 
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where  -fi  refers  to  each  of  the  series  in  braces  in  Eq.  (147),  numbered  Sequen- 
tially. Each  of  the  m,4  terms  in  -pt  and  *p  are  OCmS)’1  and  are  multiplied 
by  the  loading  at  the  leading  and  trailing  edges,  respectively.  In  addition, 
the  terms  in  -b  contain  the  exponential  factor  £ ,viti  L‘  while  those  in 

<*> x r*  * 4'  r, - 3 S 

■P;  contain  e"  w • . Away  from  the  leading  and  trailing  edges, 

the  presence  of  these  factors  accelerates  the  co:  '"'T,<7ence  of  the  series. 
However,  at  the  leading  or  trailing  edge,  one  or  the  other  pf  the  exponentials 
approaches  unity  and  the  series  c.an  diverge.  At  the  trailing  edge  Ac  should 
vanish  according  to  the  Kutta  condition,  and  the  convergence  properties  of 
the  series  in  -p  depends  on  the  behavior  of  Ap  as  r — ► . At  the  leading 

edge,  the  linearized  analysis  contains  a singularity  in  the  loading,  and  so 
the  series  for  -ft  can  be  expected  to  diverge  there.  The  behavior  at  the  lead- 
ing and  trailing  edges  i;  more  difficult  to  understand  than  it  was  in  the 
thickness  case  where  the  slopes  of  the  thickness  profile  are  well  drained. 
Because  of  their  impact  on  lifting  surface  calcul  tions,  these  terms  arc 
examined  further  in  the  next  section. 


The  series  represented  by  -t>,  contains  the  integral  which  remains 

after  the  integration  by  parts.  Integration  of  the  exponential  factors  over 

• / 

2„  introduces  a factor  which  is  O(mQ)  and  sc  each  term  in  ps  is  at  most 
O (wd)  . It  has  be>  n verified  that  successive  integration  by  parts  produces 
terms  of  higher  order  in  6 t>8)  1 and  so  this  series  is  convergent. 


Co  r>  that 
First  the 


Mext  consider  the  series  denoted  by  . It  is  this  contribution 
will  be  shown  to  contain  the  discontinuity  across  the  blade  surfaces 
expression  for  is  rewritten  in  the  form 
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The  expression  for  is  now  seen  to  be  -S(r,£), 
defined  in  Eq.  (121),  if  we  identify  f (ra ) in  Eq. 
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the  double  series  function 
(121)  with  the  quantity 


When  the  final  expression  obtained  for  S(r,t?,'  in  Eq.  (128)  is  used  here, 
then 
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Another  integration  by  parts  may  be  done  in  the  second  tern  in  p+  . The 
integrated  term  vanishes  as  a consequence  of  the  boundary  condition  that 
vanishes  at  the  hub  and  the  tips.  The  final  result  for  -p>4  is 
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The  term  containing  the  generalized  function  £■  is  the  term  which 
contains  the  discontinuity  in  -p  ; all  of  the  remaining  terms  converge  uniformly 
for  all  t?  away  from  the  leading  or  trailing  edge  and,  hence,  are  continuous 
across  the  blade  surfaces.  If  the  difference  in  p across  the  blade  surfaces 
is  formed,  p ( r,  2,  (r,  H , £ ’ — ) , then  because  of  the  nature 

of  the  function  £ ' , it  can  be  seen  that  the  definition  of  given  in  Eq.  (133) 
is  recovered. 
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Collecting  the  evaluations  of  the  various  terms,  the  expression  for 
the  perturbation  pressure  at  points  within  the  blade  row  is 
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When  the  first  two  terms  are  omitted,  this  expression  for  -p  is  valid 
upstream  (if  < 0)  and  downstream  (2  > cft)  of  the  blade  row.  Thus,  away  from 
the  leading  and  trailing  edges,  -p  (and  hence  vs  ) is  continuous  in  these  re- 
gions, and  in  particular,  across  the  blade  wakes. 

b.  Kadial  Velocity  Component 

Next,  the  properties  of  the  expression  for  the  radial  velocity  are 
examined,  focussing  on  the  terms  due  to  the  trailing  vortex  wakes.  The  radial 
velocity  is  tangential  to  the  wake  surfaces  and  should  be  discontinuous  across 
them.  This  behavior  can  be  demonstrated  simply,  if  we  consider  points  far 
downstream  of  the  rotor  where  all  other  terms  except  the  wake  term  have  de- 
cayed and  Eq.  (146)  becomes 
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Taking  the  derivative  of  the  expression  for  in  Eq.  (151]  and  using  the 
result  in  Eq.  (153) , then 
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this  point  we  introduce  the  blade  circulation,  f , which  is  defined 
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where  j(r,S)  is  the  local  vortex  strength. 
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Substituting  for  dVs  in  terras  of  Ap  using  Eq.  (53)  and  transforming  the  stream- 
wise  integration  to  an  integration  over  £ , the  circulation  is  expressed  in 
terms  of  Ap(r,Z)  as 


93 


f c*- 

/ £\  -p  ( r,  ?. ) a c 


CI5' 


Thus 


Vr(Z 


— » 4 CO  ) 


B_  d£_ B_ 

2TI  J-  dr  n U 


— 

z z 


* i • i u i + / uj'r  v-1 1 

L \ rrB  J \ U ) J 

j ir~c  [r°  ^ R„s*(<r0)  d 


^(rr.Rt f)  f'  r j J'  3 
( mfl;  1 


r 

o 

(158) 


The  term  produces  the  expected  discontinuity  in  trr  across  the  wakes  while 
the  second  term  is  continuous.  Forming  the  difference  in  vr  across  the 
blade  wake  locations  according  to  the  same  convention  adopted  for  A£>, 

' 0 H~4  - X , Z?t  \ 

AVr  - Vr  (£  “ -g*-  ) - vr  (tr  = -g~  > ..  we  get 

dP 

Aur  = (159) 

Notice  that  in  this  convention  r is  negative  when  &p  is  positive  and  work  is 
done  on  the  fluid.  This  result  for  the  jump  in  vr  across  the  wakes  is  the 
same  as  for  an  isolated  wing.  The  difference  here  is  that  the  trailing  vortex 
wakes  are  helical  rather  than  plane  surfaces. 


c,  Normal  Velocity  Component 


The  final  task  of  this  section  is  to  develop  further  the  expression 
for  the  normal  component  of  the  perturbation  velocity.  As  with  the  expression 
for  p , the  terms  in  vn  are  numbered  to  facilitate  the  following  discussion. 
From  tq.  (145) 
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The  term  labeled  ( irn)  contains  the  wake  term,  and  can  be  rewritten  as 
follows  : 
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where  it  is  recalled  that  the  prime  on  the  nn  summation  denotes  that  there  is 
no  m ,0  term,  The  -^-summation  can  be  rearranged  using  the  same  techniques 
employed  to  evaluate  the  series  function  5‘  (r,  t)  in  Section  jv-C-3.  The  resulting 
expression  for  C>rn)l  , is 
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The  A summation  which  includes  the  first  term  in  braces  is  just  the  Pourier- 
Bessel  expansion  for  Afi(r,  2„)  . The  corresponding  ■>!  summation  also  can  be 
evaluated  if  the  appropriate  m«o  term  is  added  and  subtracted  sc  that  the 
following  relation  can  be  used. 
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Then  < 1 , , becomes 


95 


1 * 


( lr„  ). 


/ L>r  Y c / u‘r  \ 

T^fr-  ?/(<■  Hr)[  *«'-*■> t*  70 


r: 


’ rr,  3-k 


.2.)  H<-2 

-2.) 

oif,  t 

i m 8 £ 

i')  e 

h'(£ 

-2.) 

r/* 

)\ 

L'  m0  > 

1 + 

^ u ) . 

B 


A A 


£ £ 

Kn  * - f>  -A  * >' 


(164) 


—■  K 77-  t-  (r, ) | R^e^(0'.')  dr.di* 
ec  L er<. 


The  first  term  precisely  cancels  the  blade  force  term  in  vn  , represented  as 
a volume  distribution  of  dipoles  in  Eq.  (141).  Thus,  we  can  drop  the  prime 
on  vn  and  the  expression  is  valid  anywhere,  including  on  the  blade  surfaces, 

The  next  operation  performed  on  xrn  is  an  integration  by  parts  over 
I.  in  (vn\  • 
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Except  near  the  leading  and  trailing  edges,  the  first  three  terms 
in  usn)2  are  well-behaved.  In  the  integrated  terms  evaluated  at  i K - 2.  , 
the  last  three  terms  in  Eq.  (165),  the  r -dependence  has  been  factored  in  a 
special  way  in  order  to  demonstrate  that  vn  is  continuous  across  the  blade 
surfaces.  The  first  of  these  terms  comes  from  the  ^-contribution  to  ir, 
Eq.  (137);  the  other  two  of  these  terms  derive  from  the  Id)  contribution  to 

vn  . The  first  of  the  last  three  terms  is  just  (u> '‘/(J)  ~~r  times  the  dis- 

Ptm  Ufi 

continuous  part  of  -p  , (see  -p4  in  Eq.  (H9J.  Hence  if  vn  is  to  be  continuous 
across  the  blade  surfaces,  then  the  last  two  terms  must  cancel  this  discon- 
tinuity. in  order  to  simplify  the  demonstration  that  such  is  the  case,(v„) 
is  further  subdivided. 
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where  ( !/„ ) includes  the  terms  containing  L rnBi.  and  can  be  written, 

using  the  definition  of  in  Eq.  (lt>7)  , as 
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Now  the  series  containing  the  first  two  tens  in  braces  are  in  the  s 
as  the  doubly  infinite  series  evaluated  already.  The  third  tern  in 
can  be  manipulated  in  much  the  same  way;  the  factor  containing 
written  using  Eq.  (122)  and  the  quantity  (~ —)  / [^  + ] is 

subtracted  to  it.  When  these  operations  arc  completed,  the  form  of 
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If  the  results  for  in  Eq.  (121)  and  -p,  in  (149)  are  used  to  evaluate 

the  series  corresponding  to  the  first  two  terms  in  braces,  the  respective  con- 
tributions to  (ir„)4  which  contain  the  generalized  function,  , cancel.  The 

differential  equation  satisfied  by  the  radial  eigenfunctions,  Eq.  (74).  can 
be  used  to  substitute  in  the  third  term  and  then  an  integration  by  parts 
carried  out.  The  final  expression  for  <z/„)4  is 
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The  series  in  Cv^.'i  ^ are  uniformly  convergent,  as  are  the  remaining 
terms  in  vn  , and  so  we  conclude  that  vn  is  continuous  across  the  blade  sur- 
faces. As  a prelude  to  evaluating  the  full  expression  for  vn  on  the  blade 
surfaces,  in  order  to  establish  the  integral  equation  for  the  blade  loading, 
we  collect  the  result  for  v„  . 
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where  (vn)a  , giver,  in  Eq.  (171),  is  not  written  out  because  it  vanishes  on 
the  blade  surfaces.  The  quantity  3 *.(r,  ra ) , defined  by  Eq.  (166), 
has  been  rationali7ed  and  written  in  the  form 
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where  is  defined  in  Eq.  (167).  The  coefficients  of  the  different  radial 

factors  m Eq.  (173),  c'^3  * (i , i, ) , arc  complex  quantities,  i.e., 
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The  real  and  imaginary  parts,  Hm0t  and  (e,  i, ) , are  written  out  in 

Table  3. 


The  formulation  of  the  loading  problem  with  the  pressure  dipole 
approach  is  now  complete.  Expressions  given  in  Sections  IV-D1  and  IV-D2  for 
the  pressure  and  velocity  fields  have  been  shown  to  satisfy  global  mass  and 
momentum  conservation,  and  to  display  the  correct  behavior  at  the  blade  and 
wake  locations.  The  result  given  in  Eq.  (172)  fox  the  normal  or  upwash  velocity 
is  used  next  to  develop  a direct  lifting  surface  analysis. 
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TABLE  3 


Coefficients  in  the  Expression  for  the  Normal  Velocity 
Component  of  a Lifting  Rotor,  Eq.  (172) 
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E. 


DIRECT  LIFTING  SURFACE  THEORY 


The  formal  solution  for  the  loading  contribution  to  the  flow  field 
of  a rotor  can  now  be  applied  to  the  direct  lifting  surface  problem.  An 
integral  equation  is  dc^iced  which  relates  the  unknown  blade  loading  to  a 
specified  blade  geometry.  This  derivation  is  accomplished  by  evaluating  the 
normal  component  of  the  perturbation  velocity,  Eq.  (1723,  at  the  blade  surface 
and  using  the  flow  tangency  condition  to  relate  it  to  the  blade  camber  line. 
The  resulting  linear  integral  equation  for is  then  nondimensionalized,  and 
what  appears  to  be  a promising  solution  procedure  is  discussed. 


1 . Formulation  of  the  Integral  Equation 


Due  to  the  periodic  nature  of  Eq.  (172),  and  the  fact  that  all 
blades  are  assumed  identical,  it  makes  no  difference  on  which  of  the  blade 
surface  streamlines,  w - y - 0,  1,  2...,  we  choose  to  specify  cor.diticr 

Two  simplifications  are  immediately  apparent  when  vn  is  evaluated  at  a blade 
surface:  Cv„)«  , as  given  by  Eq.  (171),  vanishes,  and  in  the  other  terms 
of  Eq.  (172),  the  complex  exponentials  involving  £ become  unity. 


At  this  point  some  discussion  of  the  integrated  terms  proportional 
to  o)  and  df  (r0  , Ca ) is  in  order.  Analogous  to  subsonic  isolated 

airfoil  theory,  the  Kutta  condition  is  assumed  to  be  satisfied  at  the  trailing 
edge;  i.e.,  the  loading  d-p  vanishes  so  that  this  term  need  no  longer 

be  carried.  The  leading  edge  term,  on  the  other  hand,  poses  a problem  of 
interpretation.  The  loading  on  a subsonic,  isolated  airfoil  is  known  to  ap- 
proach infinity  as  the  minus  one-half  power  of  the  distance  from  the  leading 
edge.'’3  Again,  we  car.  expect  analogous  behavior  in  the  present  situation,  im- 
plying that  the  quantity  is  singular.  Tnis  also  implies  that  the 

integrand  containing  ih-p/dla  now  contains  a minus  three-halves  power 
singularity,  which  is  not  integrable  in  the  usual  sense.  However,  note  that 
the  original  expression  for  i /„  before  performing  the  integration  by  parts 
with  respect  to  , Eq.  (145) , involved  only  integrations  over  A-p(rc,zo)  . 


102 


Hence  it  contained  only  an  integrable  singularity,  and  would  thus  give  a finite 
result.  This  strongl>  suggests  that  the  infinities  in  Eq.  (172)  are 
self-car.celling,  and  that  only  the  finite  part  of  the  expression  remains. 


Following  Mangier,  the  singularity  in  the  integral. d can  Do  djs 
played  explicitly,  with  the  other  factors  being  expended  in  Taylor  series  in 
the  vicinity  of  the  leading  edge.  Such  a procedure  shows  that  the  infinities 
do  indeed  cancel,  and  that  the  appropriate  way  to  write  the  tent;  in  question  is 
with  the  leading  edge  term  proportional  to  Ap\r0,£),  and  the  last  integral 
taken  with  f as  the  lower  limit  of  integration,  where  t‘  is  some  small  but 
finite  number.  This  will  always  yield  a finite  result,  and  amounts  to  neg- 
lecting the  contribution  from  a strip  of  width  6 near  the  leading  edge.  How- 
ever this  contribution  is  easily  shown  to  be  of  order  £ 1/1  and  hence  can  be 
made  negligible  by  chocsing  £ sufficiently  small. 


It  is  also  useful  to  eliminate  any  complex  quantities  from  the 
expression  for  , since  on  physical  grounds  it  must  be  a rcai  quantity. 

This  can  be  verified  from  Hq.  (172)  by  first  recognising  that  the 
are  real  and  are  odd  functions  of  the  index  m ; thus,  the  product  RrrSttf) 
is  even  with  respect  to  on  . From  here  it  is  easy  to  show  that  the  1 •nn.’i  terms 
are  the  conjugates  of  the  corresponding  (♦m)  terms,  and  so  the  doubly  in- 
finite sum  over  m will  always  yield  a real  result  for  7 rn  . This  suggests  the 
calculation  can  be  shortened  considerably  by  considering  only  nr,  > o , i.e-. 
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where  <5rr,  „ , the  Kronecker  delta, 
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is  used  to  avoid  counting  the  zeroth  tern  'Nice. 
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The  ne^t  stop  is  to  relate  ir,  on  the  left-hand  side  of  liq.  (112) 
to  the  prescribed  nlade  geometry,  which  'n  the  loading  case  means  the  blade 
camber  line.  (As  used  here,  the  camber  line  includes  any  mean  angle  of  attack). 
As  pointed  our  in  Section  IV-C  i,  this  is  not  as  straightforward  as  it  is  in 
isolated  airfoil  theory  because  of  the  so-called  "camber  due  to  thickness". 

A row  of  blades  with  this  camber  distribution  is  by  definition  unloaded,  and 
so  represents  the  zero  reference  in  any  loading  calculation.  Thus,  v„  in  the 
loading  case  must  be  continuous  across  the  blades  and  satisfy  the  flow  tangency 
condition  with  respect  to  the  following  camber  line, 

*nCL<t-',r)  - yc<'s.r/  - ,r)  (174) 

where  4l{  is  the  geometric  camber  of  Figure  21  and  “fy  t is  the  thickness-induced 
camber  given  by 
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Here  (lr„  )eowT  is  the  portion  of  vr  in  Eq.  (130)  which  is  continuous 
across  the  blades,  and  the  integral  is  understood  to  be  along  at  constant 
radius . 


For  a blade  with  thickness,  it  is  assumed  now  that  has  been 

determined  so  that  is  known;  alternatively,  one  can  consider  the  results 

below  as  applying  to  a blade  with  zero  thickness,  in  which  case  4lcr  is  ob- 
viously zero.  In  either  case,  the  left  hand  side  of  F.q.  (172)  is  replaced  by 
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to  injure  that  the  flow  remains  parallel  to  the  loading  cam!  er  line.  Here 

2 . r">  is  understood  to  be  ^?Ci_  (5  •=  yi  + [~j~)  ?-  . r ) when  both  s and  £ are 
assumed  to  originate  at  the  leading  edge  of  the  reference  blade,  and  use  has 
been  made  of  bq.  (Sla). 


Whan  all  the  above  steps  are  taken,  Kq.  (172)  yields  the  following 
integral  equation  for  &■{)  : 
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where  for  convenience  we  have  defined 

Jt>  , . n'-t)  rrtQ  wC  , _ . f'nB  u;  , 

^ mg-k  C°<>  JJ  i2'2.'  - ^rrtOi  < 20  ^ I.  H -H„  ) (.178) 


for  X « 1,2,3;  the  fl(^aJk's  and  6^ 4's  are  as  given  in  Table  3. 
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Before  considering  the  nature  of  this  equation,  it  is  convenient 
to  nondimensionalite  the  variables.  In  the  radial  direction,  we  retain  the 
use  of  <r ■—  as  the  dimensionless  variable.  For  the  axial  direction,  in 

't 

keeping  with  the  use  1 convention  in  isolated  airfoil  theory,  we  nondimen- 
sionalize  by  C^/x  and  shift  the  origin  to  midchord.  Thus,  we  define 


c » 


x 

so  that  now  the  blades  lie  between  A * cr  ± > and  -i  $ x * / . For  convenience 
in  specifying  the  blade  row  geometry,  we  also  define  the  parameters 
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<PT  is  the  inverse  of  the  usual  flow  coefficient  at  the  tip;3*"  f}r  can  be  re- 
lated to  either  the  aspect  ratio  or  the  solidity  at  the  tip  (based  on  Cx 
rather  than  c ) : 


Aspect  ratio 
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Solidity  at  tip  - — 

X rc  r tt 


(182) 


Ke  next  note  that  since  quantities  like  , Lrmg^  , ft**ait  , 8^3* 

and  o'm&A  are  already  dimensionless,  we  need  only  rewrite  them  in  terms  of 
the  above.  Hence,  we  get 
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The  expressions  defining  and  8^g*  can  be  used  as  they  stand  in  Table  3. 


As  for  the  dependent  variables,  such  as  velocity  and  pressure,  it 
is  natural  to  normalize  them  by  the  axial  velocity,  U , and  dynamic  pressure, 
ip«>U  i respectively.  Thus  we  define 
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Kith  these  definitions,  Eq.  (177)  becomes: 

* ’ r1 

i d\  d <r.  Ap  (o„  , x„ ) <ra 


e 4>*  nr  r 


e,<r 


xu  fi  0 - Ax)  (t*$;  <TJ)  2 J 


J, 


3 fi  { h ^ 2 ^ ^ ^ i ^ i r m ~ 

-----  — *j  / £p  (<r, , x„)  i.'cx-x#)  d xo  ~ ~ - 2 I c/x  £ 2. 

[J  t (&9r)  J,  m.i  *•> 


<*  n <r 


R^B  + Ccr)  H (X  - X 


D 


ry  Sp 


6 


* dr.(c 


•»  tt  4>t  p*  (i  + <t>*  <r*  )r/l  J 
I A^D  (oi  ,-1+e)  e 


^ A 
d A*t> 

- ^ 
? A 

i?<ra  / 

*7  «*» 

T<* 

I I 

*o 

A a*  D 


PT  (*  * d r a 

£»,**<*.  - O(<r<r0) 


r + ’ 
/ 

3 Ap  ^ - ^ B <r 

J 

9 

X 

-i  +e 

u> 

/ <r 

T \ -J  <3>  ,1 

~)  - <Pr  Qm9t<x,-’)(r<r.)- 


d> 


\<Pr  G~s*(><> 


U)  /x-  ff~  \ -l  fj>  j I 

+ x#)  ( — + — y + <#>r  X.)(rrm)  J 


107 


(187) 


The  next  step  is  to  consider  the  behavior  of  the  integrand.  In 
isolated  airfoil  theory,  singularities  generally  appear  as  the  source  ( <T0  , >„ 1 
and  field  point,  approach  one  another,  and  we  must  expect  such  behavior 

here.  In  our  case  however,  they  would  appear  as  divergent  series,  so  that 
the  nature  of  the  singularity  is  not  immediate ly  obvious.  Preliminary  con- 
sideration of  the  series  in  Lq.  (187)  when  the  source  and  field  points  coincide 

- t 

suggests  that  for  fixed  A , the  terms  in  the  integrand  decay  only  as  (mfi)  , 
and  for  fixed  mfi  as  (4)''  , thus  indicating  potential  convergence  problems. 

(It  should  be  noted  here  that  this  in  no  way  contradicts  the  earlier  statement 
that  the  rvi  scries  occurring  in  Eq.  (172)  for  v„  is  uniformly  convergent. 

There  we  were  considering  the  integrated  expression  for  -t/„  , and  hence  could 
make  the  use  of  the  fact  that  the  integration  over  z„  of  the  exponential  argu- 
ments would  yield  an  extra  factor  of  (r^B)  Here  we  are  considering  the  be- 

havior of  the  integrand  near  x » X0  , at  which  point  the  exponentials  ail 
reduce  to  unity,  and  hence  cannot  aid  in  the  convergence). 


The  nature  of  these  singularities  has  not  yet  been  determined,  and 
so  a detailed  discussion  of  the  convergence  of  these  series  is  deferred  until 
then.  Once  the  form  of  the  singularity  is  determined,  it  is  anticipated  that 
its  influence  can  be  subtracted  out  and  isolated,  much  as  the  discontinuous 
portions  of  and  Vr  were  isolated  in  Section  IV-D-3.  In  addition  to  facili- 
tating the  evaluation  of  the  improper  integrals  which  occur,  the  isolation  of 
the  singularity  will  serve  another  very  useful  purpose.  In  any  numerical 
scheme,  the  infinite  series  over  the  azimuthal  and  radial  mode  numbers  must  be 
truncated  after  a finite  number  of  terms;  since  all  scries  would  be  convergent 
after  the  above  separation,  a much  better  estimate  of  the  truncation  error 
should  result. 

2 . Progress  Toward  the  Solution 

Despite  the  fact  that  this  singularity  has  not  yet  been  isolated, 
some  progress  has  been  made  toward  the  inversion  of  pq.  (187).  This 
equation  is  obviously  too  complex  to  hold  any  hope  for  an  analytical  solution, 
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and  so  various  numerical  schemes  employed  successfully  in  past  1 ifttng-suriace 
studies^ ' have  been  reviewed.  The  schemes  all  appear  to  fall  into  one  ot 
two  basic  categories.  The  first  involves  the  idealization  of  the  distributed 
loading  as  a lattice  of  discrete  loading  elements  (j.e.,  vortices,  dipoles, 
etc.)  whose  strengths  are  constant,  but  initially  unknown.  The  second  ap- 
proach is  to  represent  the  continuous  loading  as  a double  series  of  suitable 
functions  in  the  chordwise  and  spanwise  variables,  the  coefficients  in  which 
are  initially  unknown.  This  is  also  sometimes  referred  to  as  a kernel  func- 
tion approach. 

In  either  case,  one  then  requires  that  the  velocity  be  parallel  to 

the  blade  surface  at  each  of  a set  of  suit.abl>  chosen  collocation  points. 

This  reduces  the  integral  equation  to  a set  of  simultaneous  linear  algebraic 

equations  which  can  be  expressed  as  a matrix  equation.  The  rank  of  the  matrix 

equals  the  number  of  unknowns,  which  in  turn  equals  the  number  of  collocation 

points  for  the  system  to  be  deteiminant.  One  can  also  choose  the  number  of 

collocation  points  to  exceed  the  number  of  unknowns,  in  which  case  the  boundary 

56  5"1 

condition  at  the  surface  can  only  be  satisfied  in  a least-squares  sense.  ' 

The  solution  of  a matrix  equation,  whether  for  the  unknown  strengths 

of  the  source  lattice  or  for  the  unknown  coefficients  in  the  series  expansion, 

is  thus  common  to  both  approaches.  This  portion  of  the  solution  is  relatively 

straightforward  and  rapid.  It  is  in  determining  the  elements  of  the  coeffi- 

£6 

cient  matrix  that  the  most  time  is  usually  consumed,  and  it  is  here  that  the 
two  techniques  differ  significantly,  with  each  having  advantages  in  certain 
problems . 

The  greatest  advantage  of  the  lattice  methods  is  that  no  numerical 
quadratures  are  needed  to  determine  the  coefficient  matrix.  Moreover,  pro- 
vided the  source  and  collocation  points  are  never  coincident,  the  problems 

S8 

mentioned  above  in  treating  the  singularity  are  evidently  avoided.’  However, 
since  one  only  solves  for  the  loading  at  a set  of  discrete  locations,  there  is 
some  ambiguity  involved  in  how  best,  te  interpolate  for  t.ho  loadings  at  other 
positions.  (A  general  discussion  of  this  point  is  given  in  Chapter  3 of 
Reference  5V.) 
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We  feel  that  the  series  expansion  approach  has  several  overriding 
advantages  in  the  ; resent  problem,  the  principal  one  being  that  it  enables 
one  to  express  the  loading,  as  well  as  its  derivatives,  uniquely  at  every 
point  oil  tile  surface.  Through  judicious  choice  of  the.  loading  functions  in 
the  series  expansion,  one  can  also  include  the  appropriate  behavior  at  the 
leading  and  trailing  edges  automatically.53,5'’  The  principal  disadvantage  is 
that  numerical  quadratures  usually  are  required  lor  the  elements  of  the  co- 
efficient matrix,  which  for  our  kernel  function  would  likely  be  rather  time- 
consuming  . 


fortunately,  there  appears  to  be  a strong  possibility  that,  once 
the  singularity  has  been  treated  separately,  the  remaining  spanwisr  integrations 
can  all  be  done  analytically.  To  see  this,  note  that  the  <r„  integrals  ir. 

Eq.  (187  j .involve  only  products  of  </„'  ' with  A-f  or  its 

derivatives . Furthei , if  we  -xpand  the  radial  dependence  of  Ap  in  a simple 
power  series,  then  the  integrals  vi.ll  still  involve  only  products  of  Rmg^  with 
integral  powers  of  if*  . Then,  since  the  RMtft  arc  linear  combinations  of 
and  (see  Appendix  A),  use  can  be  made  of  the  indefinite  integral:' 
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where  ix  and  n are  t-ach  integers,  the  n are  boriuv.ei's  functions  -nd  % is  a 
dummy  vai  iafcle  for  K,,*  <r,  . Reference  59  gives  an  easily  used  asymptotic  ex- 
pansion for  the  evaluation  of  the  „ valid  for  large  argument.  Since 
r'*t.  » where.  * m3  , can  be  shown  to  be  of  the  same  order  as  the  number  of 
blades,  this  should  be  applicable  in  the  majority  of  cases  encountered  in  the 
present  problem. 
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One  further  point  needs  to  be  made  regarding  the  radial  variation 
of  Ap  , and  that  concern, s its  behavior  as  the  hard  walls  at  cr  * A and  1 are 
approached.  In  isolated  airfoil  theory  the  loading  at  the  tips  must  vanish 
since  no  pressure  differential  can  exist  off  the  airfoil.  However,  in  the 
present  problem,  since  we  have  not  allowed  for  th : effects  of  any  clearance; 
between  the  blades  and  the  walls,  the  hard-wall  boundary  conditions  require 
instead  that  there  be  no  radial  pressure  gradient  there,  so 
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This  requirement  can  be  satisfied  easily  by  a power  series,  and  is  another 
advantage  of  using  such  an  expansion  in  the  radial  coordinate.  In  fact,  if 
one  assumes  an  expansion  with  N unknown  coefficients,  two  of  these  can  be 
expressed  as  linear  combinations  of  the  remaining  (N~2)  using  the  two  condi 
tions  embodied  in  Eq.  (189).  This  allows  one  to  represent  the  radial 
variations  with  a power  series  two  orders  higher  than  the  number  of  unknown 
coefficients  one  is  willing  to  include. 


As  for  the  chordwise  variation,  it  is  very  doubtful  that  a simple 
power  series  expansion  in  would  be  adequate.  This  I:  because  of  the  need 
to  represent  the  behavior  near  the  blade  leading  and  trailing  edges,  where 
Ap  and/or  have  singularities.  Accordingly,  an  expansion  which  represents 

such  behavior  explicitly  will  likely  prove  necessary,  i.e. 

A-M-tT.*.)  ■=  P0(tr.)  '[p— ~ + Vt  -xf  £ fl  (o O90) 

V 1 4 A „ L 

where  each  of  the  is  a power  series  in  as  discussed  above,  and  the  P 
ir  a suitarly  chosen  set  of  finite  polynomi  ns,  each  of  order  i,  . Such  a 
representation  automatically  exhibits  the  appropriate  (i  * xo)*^  behavior 
near  the  leading  and  trailing  edges,  respectively.^  Unfortunately,  it  also 
requires  the  use  of  numerical  integrations  with  respect  to  xe  . The  effi- 
ciency and  accuracy  of  these  integrations  can  be  affected  by  the  form  as.umcd 
for  the  set  of  polynomials  P and  various  alternatives  discussed  in  Refs.  55 
and  60  are  being  considered.  Once  tne  terms  containing  the  singularity  that 
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occurs  at  the  coincidence  of  the  source  and  field  points  has  been  isolated  and 
the  mathematical  o.iar  icter  of  the  nonsingular  portion  of  the  integrand  has 
been  established,  a choice  for  the  form  of  the  P will  be  made  and  efforts  to 
obtain  numerical  results  will  begin. 

F.  CONCLUDING  REMARKS 

A linearized  direct  lifting  surface  theory  has  been  formulated  for 

the  compressible,  three-dimensional  flow  through  a rotor  of  specified  geometry 

in  an  annular  duct.  In  arriving  at  this  formulation,  the  overall  subject  of 

the  linearized  analysis  of  three-dimensional  compressor  flows  was  reviewed  in 

detail.  Tills  review  was  made  necessary  by  difficulties  encountered  in  early 

phases  of  the  work.  Initial  attempts  to  derive  the  loading  contribution  to 

the  flow  field  of  a rotor  by  the  superposition  of  pressure  dipole  solutions 

met  with  two  problems.  The  solution  for  the  rctor  pressure  field  did  not 

yield  a net  pres sue  rise  between  points  far  upstream  and  downstream  of  the 
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rotor.  Furthermore,  it  did  not  agree  with  the  result  obtaind  by  Mamba  in  a 
similar  investigation. 

The  error  in  the  pressure  rise  was  traced  to  the  omission  of  a term 

in  previous  versions  of  the  fundamental  solution  for  a point  source.  The 

missing  term  in  the  source  solution  also  led  to  an  error  in  the  dipole  solution. 

When  these  singularity  solutions  were  corrected  and  the  result  for  the  pressure 

field  revised  accordingly,  then  a nonvanishing  pressure  rise  across  the  rotor 

was  found  which  agreed  with  the  vortex  theory  result  obtained  by  Okurounmu 

4; 

and  McCune. 

Both  the  thickness  and  loading  contributions  to  the  flow  field  of  a 
rotor  were  rederived  using  the  corrected  singularity  solutions.  The  implica- 
tions of  the  additional  term  in  the  source  solution  to  published  reports  on 
the  thickness  problem  have  been  described.  The  solution  for  the  loading  problem 
has  been  applied  to  formulate  a direct  lifting  surface  theory.  There  are  two 
remaining  differences  between  our  formulation  and  that  given  by  Nainba. 

First  of  all,  there  is  the  difference  which  results  because  t he  additional 
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term  we  found  to  be  required  in  the  dipole  solution  is  not  present  in  his. 
Secondly,  our  formulation  does  not  contain  his  scale  factor,  a function  of 
radius  which  modifies  the  dipole  strength.  Wo  ha."  concluded  that  this  factor 
should  not  be  included  on  the  basis  of  several  checks  made  on  our  analysis. 

Both  the  thickness  and  loading  contributions  to  the  rotor  liow  field  have  been 
shown  to  satisfy  the  global  conservation  of  mass  and  of  the  axial  components 
of  linear  and  angular  momentum.  In  addition,  the  velocity  components  have 
been  shown  to  display  the  correct  behavior  at  the  blade  surfaces  and,  in  the 
loading  case,  across  the  trailing  vortex  wakes. 

Progress  has  been  made  toward  obtaining  numerical  solutions  of  the 
integral  equation  for  the  unknown  blade  loading.  A kernel  function  method  lias 
been  chosen  as  the  best  approach;  the  loading  is  expanded  in  a double  series 
of  suitable  functions  in  the  chordwise  and  spanwise  variables,  the  coefficients 
in  which  are  unknown.  The  expected  loading  behavior  near  the  blade  loading 
and  trailing  edges  is  ensured  by  the  inclusion  of  appropriate  factors  outside 
the  double  sum.  The  choice  of  the  form  of  the  expansion  functions  used  to 
represent  the  axial  and  radial  variations  is  deferred  until  the  singularities 
in  the  integral  equation  have  been  isolated.  These  singularities  occur  in  the 
kernel  function  when  the  source  and  field  points  coincide.  It  is  felt  that  the 
same  techniques  tied  to  isol-  _c  the  discontinuities  in  the  pressure  and  velocity 
fields  can  be  ose-.  to  determine  the  nature  of  these  singularities,  and  work 
along  these  lines  has  begun.  The  isolation  of  these  singularities  should 
accelerate  the  rate  of  convergence  of  the  series  in  the  remaining  portions  of 
the  integrand,  and  therefore,  facilitate  numerical  evaluation  of  the  solution. 
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SUCTION  V 

SUMMARY  AM'  CONCLUSIONS 


A theoretical  and  experimental  investigation  of  discrete-tone  rotor- 
stator  interaction  noise,  and  the  development  of  a three-dimensional  li.fting 
surface  theory  for  a rotor  have  been  presented.  In  the  experimental  part  of 
the  program,  the  sound  pressure  levels  produced  on  the  outer  duct  wall  of  the 
annular  cascade  facility  were  measured  for  a rotor-stator  pair.  The  micro- 
phone probe  was  located  in  the  far  field  upstream  of  the  blade  rows  and  data 
were  taken  for  two  stator  stagger  angle  settings.  In  both  sets  of  data  the 
rotor  speed  was  varied  continuously  up  to  about  1450  rpm.  A pure  tone  signal 
was  observed  above  1100  rpm  which  was  identified  as  the  fourth-harmonic  of 
blade  passage  frequency.  For  one  of  the  stator  settings,  some  additional  data 
were  obtained  on  the  fifth  harmonic  at  about  600  rpm.  Lower  harmonics  could 
not  be  excited  at  the  rotor  speeds  attainable  in  the  facility  at  the  time  of 
these  experiments. 


A limited  amount  of  experimental  data  were  taken  on  the  time  varying 
pressure  signal  at  several  chordwi.se  locations  adjacent  to  the  blade  tips  of 
an  isolated  rotor.  These  data  were  taken  in  ar,  attempt  to  obtain  information 
on  the  loading  at  the  blade  tips  with  which  to  compare  lifting  surface  calcu- 
lations. However,  the  data  appear  to  contain  tip  clearance  effects  which  are 
not  included  in  the  theory  and  would  complicate  such  an  eventual  comparison. 


Ail  approximate  model  has  been  developed  to  predict  the  discrete-tone 
noise  produced  by  the  interaction  of  a rotor  and  a stator.  In  this  analysis 
the  aerodynamic  and  acoustic  aspects  of  the  problem  are  treated  separately. 

The  sound  pressure  levels  in  the  duct  and  total  radiated  power  at  the-  harmonics 
of  blade  passage  frequency  are  computed  in  terms  of  the  propagating  duct 
acoustic  modes.  The  amplitudes  of  these  modes  are  related  to  the  blade  row 
parameters  through  an  approximate  two-dimensional  representation  of  the  unsteady 
blade  forces.  The  main  advance  over  previous  approximate  treatments  of  rotor- 
stator  noise  is  the  incorporation  of  an  aerodynamic  model  that  accounts  for 
compressibility  effects. 
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Sample  calculations  have  been  presented  for  sound  pressure  levels  and 
totii  radiated  power  over  a wide  range  of  subsonic  Mach  numbers  and  comparisons 
of  the  inode  amplitudes  were  made  with  those  predicted  by  an  incompressible 
model.  In  addition,  theoretical  predictions  were  made  of  the  discrete-tone 
r.oise  to  be  expected  in  the  acoustic  experiments.  The  calculated  sound 
pressure  levels  were  sigiiif  icantly  lower  than  the  measured  values.  Lvidenco 
has  been  presented  showing  that  the  sound  pressure  levels  at  the  higher 
harmonics  of  blade  passage  frequency  arc  sensitive  to  the  uncertainties  in  the 
data  for  wake  mean  velocity  profiles.  This  should  not  be  true  for  the  funda- 
mental and  the  lower  harmonies  and  so,  if  the  inaccuracies  in  predicting  the 
sound  pressure  levels  of  the  fourth  harmonic  are  due  to  an  inadequate  wake 
model,  then  better  agreement  with  experiment  should  result  for  lower  hannonics. 
The  annular  cascade  facility  has  been  repowered  so  that  higher  rotor  speeds 
are  now  available,  further  acoustic  experiments  are  planned  in  order  to  obtain 
data  on  the  lower  harmonics. 

Finally,  a direct  lifting  surface  theory  has  been  formulated  for  the 
compressible,  three-dimensional  flow  through  an  isolated  rotor.  A pressure 
dipole  representation  of  the  blade  row  was  used  to  derive  the  integral  equation 
relating  the  blade  loading  to  a specified  blade  shape.  During  this  development, 
tlie  omission  of  a term  in  previously  published  versions  of  the  fundamental 
source  and  pressure  dipole  solutions  was  found  and  corrected.  As  a result,  a 
general  review  of  the  linearised  analysis  of  three-dimensional  compressor  flows 
was  presented.  Both  the  thickness  and  loading  contributions  to  the  flow  field 
of  a rotor  were  considered. 

Work  was  begun  on  solving  the  governing  integral  equation  using 
numerical  techniques  analogous  to  those  employed  successfully  for  isolated 
wings  and  cascades.  A.  kernel  function  approach  has  been  chosen  which  in- 
corporates the  expected  behavior  of  the  loading  at  the  blade  leading  and  trailing 
edges  and  at  the  duct  walls.  At  present,  the  singularity  in  the  kernel 
function  is  being  isolated  and,  once  this  is  completed,  programming  and  numerical 
evaluations  of  the  solution  will  begin. 
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APPENDIX  A 


DI-.1-  IN ITION  OF  RADIAf.  EIGENFUNCTIONS 

This  appendix  summarizes  the  definition  and  properties  of  the  ortho- 
normal radial  eigenfun  tions,  . As  used  in  Section  III 
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where  Jr,  and  Y„  are  the  Bessel  functions  of  the  first  and  second  kind,  of 
order  o . The  normalization  factor  is  chosen  such  that 
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where  6,,^  is  the  Kronecker  delta.  This  requires 
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The  set  of  eigenvalues  v^m  corresponding  to  each  function  are  solutions  of 
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which  insures  that  Equation  (8)  will  be  satisfied.  The  zeroth  eigenvalue,  vno 
is  zero  so  that  all  the  Ano  are  trivially  zero,  save  for  the  case  n ..  o ■ Then 
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The  set  of  functions  Rn>>,  is  complete  for  any  a , so  that  a function 
can  be  represented  on  the  interval  A to  1 by  the  series 
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and  a can  be  chosen  For  convenience.  These  functions  are  discussed  at  more 

28  4 

length  in  the  papers  by  McCune  ami  Tyler  and  Sofrin.  (Note  that  in  the 
lifting  surface  theory  of  Section  IV,  a somewhat  different  notation  is  used. 
There,  4 rather  than  in  , is  used  as  the  radial  mode  index,  and  the  eigen- 
values are  referred  to  as  K ^ rather  than  Vnrn.) 
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APPENDIX  B 

EXPRESSIONS  FOR  ROTOR-STATOR  INTERACTION  FORCES 

This  appendix  quotes  the  results  of  Osborne's  analysis  of  the  un- 
steady  compressible  interactions  within  a stator-rotor  pair,"  as  shown  in 
Figure  12a.  Some  typographical  errors  appearing  in  Reference  24,  and  checked 
with  its  author,  are  also  corrected.  The  notation  used  here  is  that,  of  Osborne, 
except  where  otherwise  noted.  Note  in  particular  that  "A  " was  used  for  the 
load  harmonic  index,  in  place  of  his  "rr\'\  to  keep  consistent  with  the  main 
text . 


a]  Rotor  Unsteady  Lift  Resulting  from  Passage  of  Steady 
Stator  Loads  — — — - — • 
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shape,  That  is,  if  ,Y  <xs)  describes  the  stator  blade  shape,  then 
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An  analogous  expression  defines  the  „fl, 
discussed  further  below. 


The  quantities  Au  and  rR  will  be 
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The  function  KJ  represents  the  airfoil  response  to  a gen- 

eralized Kemp-type  upwash  gust.  In  Equation  (6)  of  Reference  24.  it  was  given 
as  the  sum  of  an  infinite  series  of  products  of  Bessel  functions.  Kemp^ 
later  showed  that  the  series  could  be  summed  analytically,  and  gave  the  closed 
form  expression: 
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with  either  r or  s subscripts,  as  appropriate  above.  Here 


J (X)  = J.(X)  - i J,  (X) 

and  C(x)  is  Theodorsen's  function. 
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where  K„  and  K,  are  modified  Bessel  functions  of  the  second  kind  with  imaginary 
arguments.  This  rs  the  form  for  Ku  used  in  the  present  study;  using  it, 
Osborne's  original  numerical  results  were  reproduced  to  within  the  reading 
accuracy  of  his  curves. 
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Furthermore,  it  has  recently  been  noted'  that  Osborne's  original 

asymptotic  aerodynamic  analysis  does  not  appear  to  be  consistent  to  first  order 

in  his  small  parameter  M uo/2 n /3  , as  cl  aimed. ^ 1 After  re-doing  the  analysis 

starting  i'rou  the  integral  equation  for  blade  loading  in  unsteady  flow,  it  was 

found  that  all  ot  Osborne's  expressions  remained  valid  provided  that  X was 
, ~g 

replaced  by  X where 

X ' = [M  2 - f(M)  J a 

f (M)  « JL^.  (i  + /3)  + (i-fl)jUtM  - J!s>\  2 ■ f (o)  = 0 

again  with  either  fl  cr  i subscripts.  Accordingly,  this  revised  definition  of 

Z was  used  in  both  potential  interactions  above  for  the  present  study.  Note 

that  this  is  no  way  alters  the  fact  that  tne  equations  still  reduce  to  those 
2 

ot  Kemp  and  Sears  at  incompressible  speeds. 
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and  2,,  is  some  chosen  value  cf  XR  between  - c,  and  * c*  , usually  the  quarter- 
chord  point.  The  other  quantities  arc  given  in  a)  above.  As  originally 
derived  by  Osborne,  T , which  represents  the  airfoil  response  to  a sinusoidal 
gust,  was  given  as 
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where  J oo  is  defined  above  and  6‘ oO  is  Sears'  original  incompressible 
response  function: 

5 ( X ) = 

l X [Ko(l  X)  + K/iX)] 

However,  Amiet^  has  shown  that  if  this  expression  is  multiplied  by  the  phase 
factor  e ^ the  result  more  closely  approximates  numerical  solutions 

for  the  loading.  Thus,  in  the  present  study 

i 4 o*>Clk 

T(jLxK,XaK)  = J a?„)  5 fi  n„)  <? 

where  f(M)  is  defined  in  fc)  above.  Note,  however,  that  use  of  Osborne's 
original  x , and  not  x‘  , is  retained  in  this  expression;  also,  the  magnitude 
remains  unaffected.  Again,  the  above  results  reduce  to  these  of  Kemp  and 
sears  in  the  limit  M -o  . (The  factor  in  Qi  is  a phase  shift  reflecting 

Osborne's  alternate  choice  of  the  blade  row  positions  at  the  instant  t-o  ). 


Having  quoted  the  aerodynamic  expressions  appropriate  to  a stator- 

rotor  pair  shown  in  Figure  12a,  we  now  consider  how  these  can  be  applied  to 

24 

the  rotor-stator  geometry  of  Figure  I2b.  As  noted  by  Osborne,  the  only  thing 
in  his  two-dimensional  model  which  distinguishes  the  rotor  from  the  stator  is 
that  by  definition  the  rotor  is  the  blade  row  which  moves  with  respect  to  the 
duct.  Since  the  row  stagger  angles,  blade  profiles,  etc.  all  appear  as  input 
parameters  n Osborne's  analysis,  these  are  readily  switched.  But  the  face 
that  it  is  now  the  upstream  blade  row  which  is  moving  is  not  consistent  with 
his  model.  This  is  easily  remedied  by  a Galilean  transformation  to  coordinates 
moving  with  the  rotor,  as  in  Figure  12c,  so  that  the  stator  now  appears  in  mo- 
tion. Noting  the  sense  in  which  a„  and  as  are  defined  positive  in  Figure  12a, 
a comparison  of  Figure  12c  with  Figure  12a  suggests  that  the  interaction  forces 
in  the  reverse  geometry  of  Figure  12b  are  directly  related  to  those  of  Figure  12a 
if  one  simply  substitutes 

aR  ■»  - ce5  cts  « - a*  U - - U 

in  the  abrve  equations  and  reverses  the  subscripts  « and  S wherever  appropriate. 
At  f. rst  glance,  the  reader  might  be  bothered  that  we  appear  to  be  computing 
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the  forces  in  the  wrong  coordinate  frame,  so  to  speak,  i.e.,  with  the  stator 
moving  and  the  rotor  stationary.  But  we  are  only  concerned  with  the  force 
magnitudes  and  their  relative  phasing  between  blades,  which,  of  course,  re- 
main the  same  provided  we  always  measure  in  a Galilean  frame.  In  using  these 
forces  in  the  subsequent  acoustic  calculations,  we  revert  back  to  duct-fixed 
coordinates  £see  Equations  (13a-c)  in  the  main  text]  . 
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APPENDIX  C 


MASS  AND  MOMENTUM  BALANCES  FOR  SINGULARITY  AND 
ROTOR  FLOW  FIELDS 


Tests  performed  or.  the  source,  dipole,  and  rotor  flow  field  solu- 
tions arrived  in  Section  IV  to  determine  if  they  satisfy  global  mass  and  mo- 
mentum conservation  are  described  here.  The  conLrol  volume  employed  is  bounded 
by  the  duct  walls  and  the  annular  areas  perpendicular  to  the  duct  axis  at  up- 
stream and  downstream  infinity  (sec  Figure  19).  Since  this  control  volume 
rotates  with  angular  velocity  ..o  (in  the  negative  Q -direction) , the  conserva- 
tion laws  for  a noninertial  reference  frame  must  be  used.  The  appropriate 
integral  fonts  of  the  equations  for  conservation  of  mass,  momentum,  and 
angular  momentum  for  such  a control  volume  have  been  taken  from  Ref.  61. 

Then,  these  equations  have  been  linearized  so  that  they  are  expressed  in 
terms  of  the  undisturbed  flow  properties  , U , and  ur  , and  the  perturbation 
quantities  -p  t fi  , xrr  , ve  , and  vz  . These  equations  were  specialized 
further  to  account,  for  the  fact  that,  all  the  flow  field  solutions  decay  at 
upstream  infinity  and  satisfy  the  boundary  condition  of  no  flow  through  the 
duct  walls.  Attention  has  been  confined  to  the  expressions  for  mass  conserva- 
tion and  to  the  axial  components  of  the  momentum  and  angular  momentum  balances. 
Under  the  present  set  of  assumptions,  these  expressions  become 
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where  "W denotes  the  rate  of  introduction  of  mass  into  the  control  volume 
while  fi  and  Ji  arc  the  axial  components  of  the  net  force  ana  torque  on  the 
fluid.  In  the  discussion  which  follows,  the  flow  field  solutions  obtained 
for  a mass  source,  pressure  dipole,  rotor  blade  thickness  distribution,  and 
rotor  loading  distribution  are  shown  in  each  case  to  give  the  values  of  , 

7^  , and  J,  required  to  satisfy  these  conservation  laws. 

1 . Mass  Source 

The  flow  field  due  to  a point  source  which  is  located  at  the  point 
( fu  , B0  , ) was  given  in  Section  IV-B.  The  mass  addition  rate  into  the 

control  volume,  afYl  , due  to  this  source  is  just  the  following  volume  integral: 
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This  rate  of  introduction  of  mass  must  be  balanced  by  the  mass  iiow  out  of 
the  control  volume,  as  expressed  by  Eq.  (C-l).  Using  Eq.  (54)  to  eliminate 
p in  favor  of  -p  , that  expression  for  becomes 
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From  Eqs.  (83b)  and  (83c),  the  limiting  values  of  Vt  and  p for  £ — ■ « are 
found  to  be 
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Substituting  these  results  into  Eq.  (C-S),  it  is  easy  to  show  that  the  inte- 
gration yields  ~ p„  Q.  and  so  the  source  solution  satisfies  global  con- 
servation of  mass. 


The  presence  of  a source  in  a mean  flow  results  in  a force  on  the 
fluid.  To  first  order,  the  force  due  to  a point  source  of  strength  Q_  in  an 
undisturbed  stream  with  velocity  is 
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The  axial  component  of  this  force,  , is  pnUQ.  . This  force  on  the  fluid 
must  be  balanced  by  the  momentum  flux  and  the  pressure  forces  at  the  surfaces 
of  the  control  volume.  For  the  control  volume  being  considered  here,  this 
balance  is  expressed  by  Eu.  (C-2).  Again  eliminating  p in  terms  of  jo  , that 
equation  becomes 
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If  the  results  tor  (f*)s  and  -ps  given,  in  F.q.  (C-6)  are  substituted  in  Eq.  (C-8) , 
it  is  straightforward  to  show  that  this  expression  reduces  to  ?e  = p„UQ  , 
and  hence  the  source  solution  also  satisfies  the  axial  component  of  global 
momentum  conservation. 


In  the  present  application,  the  presence  of  a source  also  results  in 
a torque  on  the  fluid  because  the  undisturbed  stream  has  a rotational  component. 
Since  the  0 -component  of  the  undisturbed  velocity  is  ^ the  axial  component 
of  the  torque  on  the  fluid  is 
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Conservation  of  angular  momentum  requires  that  this  also  be  the  result  of 
the  integrations  over  the  source  solution  indicated  in  Eq.  (C-3) . first  we 
write  Eq.  (C-3)  in  the  form 


(C-10) 


where  j , is  the  volume  integral  end  is  the  surface  integral.  From  Eq.  (83a] 
it  can  be  seen  that  (Vr&)s  decays  for  £—•■».  h'hen  Eq.  (C-6)  is  used  for 
the  limiting  values  of  (lrt)s  and  -pf  , the  surface  integral  becomes  similar  to 
those  already  evaluated  in  the  Ojj  and  7^  expressions.  The  result  for  is 
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in  Eq.  (81),  the  following  intermediate  result  needed  in  «$,  can  be  obtained. 
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Also  using  the  definition  of  from  Eq.  (71),  the  expression  forxj,  at  this 
stage  is 
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Next  an  integration  by  parts  is  carried  Out  and  liqs.  (70)  and  (74)  are  used 
in  the  result  so  that  xS,  becomes 
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The  summation  in  (C-14),  but  tor  the  A-0  term,  is  the  Fourier-Bessel  expansion 
for  r0'V.:’.  in  terms  of  the  zeroth  order  eigenfunctions,  (<re~)  . The 
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final  result  forJ,  is 
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Using  Has.  (C-ll)  and  fC- 15)  in  Eq.  (C-10)  yields 

'Ji  * Pc*^  re  & (C-  16) 

which  is  the  result  required  for  the  source  solution  to  satisfy  conservation 
of  angular  momentum  (see  Eq.  (C-9)). 


2 . Pressure  Dipole 

The  pressure  dipole  singularity  should  result  in  a net  force  on  the 
fluid,  but  it  should  not  introduce  any  mass.  Hence,  in  this  case  the  rate 
of  mass  addition  to  the  control  volume,  , must  be  zero.  Before  substituting 
the  results  for  the  dipole  flow  field  into  kq.  (C-i)  to  check  for  conservation 
of  mass,  these  results  are  rewritten  in  a somewhat  different  form.  The  in- 
tegral of  the  dipole  pressure  field,  -pD  , along  the  stream  direction  as  given 
in  Eq.  (100a)  is  labeled  • Then,  from  Eqs.  (99)  and  (100),  the 

noimal  velocity  component  is 
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Now  if  -3—  is  written  in  terms  of  -J-) „ and  J—  ) , as  in  Eq.  (115),  and  it 
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is  recalled  that  t,  is  held  constant  in  the  integration  along  streamlines, 
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Using  Eq.  (S3)  to  relate  (lrf)D  to  -p0  and  then  using  Eq.  (101b)  to  obtain 
(irt)0  from  (Vn  )D  and  (Vf^p  , 
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In  order  to  evaluate  the  integral  in  Eq.  (C-l)  we  need  the  limiting  forms 
3 1 

of  -p0  and  j~~  for  ?.  -»  <*>  . These  quantities  arc  found  to  be 
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Upon  substituting  for(tr2)0  ar.d  pD  in  Eq.  (C-l),  the  integral  expression  for 
is  written 
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Now,  using  Eq.  (C-20a)  in  Eq.  (C-22),  ^ is  found  to  be 
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As  a result  of  the  (9 -integration,  J,  vanishes  because  the  r.  = £>  term  is  absent. 
Also,  since  for  fixed  Z,d8-oL£,  we  have 
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When  the  results  for  J , , and  are  added  together,^  vanishes  as  it 

should. 


As  discussed  in  Section  IV-B,  the  pressure  dipole  exerts  a force 
of  magnitude  D in  th ' negative  n direction  (see  Figure  20].  The  H -component 
of  this  force,  .is  D sin  yra  . The  global  momentum  balance  as  stated  in 
Eq.  (C-2)  provides  a check  on  whether  the  dipole  flow  field  is  consistent 
with  this  force  on  the  fluid.  Eq.  (C-2)  can  be  expressed  in  terms  of  the 
integrals  , and  d?r  defined  above. 

~ ^^3  + z + iV'Js  (C -27) 

Using  the  results  in  Eqs,  (C-25)  and  (C-26) , and  the  fact  that  vanishes, 
the  dipole  flow  field  is  seen  to  satisfy  momentum  conservation. 


The  force  which  the  dipole  exerts  also  results  in  a torque  acting 
on  the  fluid.  The  axial  component  of  this  torque,  J£  , is  -Dra  cos  . In 
order  for  the  dipole  flow  field  to  conserve  angular  momentum,  this  also  must 
be  the  value  of  the  integrals  on  the  right-hand  side  of  Eq.  (C-3).  In  order 
to  simplify  the  description  of  these  integrations,  Eq.  (C-3)  is  written 


whe  re 
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The  radial  velocity  component,  for  the  dipole  field  is  given  by  Eqs.  (98)  and 
(100a).  If  the  9 and  £ integrations  in  Uu  are  carried  out,  the  following 
result  is  obtained. 
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Using  Eq.  (032)  in  Eq.  (029)  and  substituting  for  A**  from  F.q.  (71),  the 
followj.ng  intermediate  expression  for  results. 
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This  same  integral  was  evaluated  in  the  angular  momentum  balance  for  the  mass 
source,  see  Eqs.  (013)  through  (015).  Using  that  result 


c u 


* £)  - 


u 


0 S&ST'  Tyo 


(0  34) 


The  evaluation  of  4-j  requires  the  result  for  ( Vg)D  , which  can  be 
obtained  from  <ys)p  and  (v-„)0by  Eq.  (101a).  Using  Eqs.  (91),  (53),  and 
(0.18)  in  Eq.  (101a), 
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The  limit  of  for  £ — * is  given  by  Eq.  (O20b)  ar.d,  since  the  0-inte- 

gration  of  this  quantity  vanishes,  J7  reduces  to 
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The  integral  d?e  is  similar  to  the  integrals  J3  , Jl4  , and  Jg  done 
in  connection  with  the  mass  balance.  Substituting  into  Eq.  (C-31)  using 
Eqs.  (54),  (C- 19) , ar.l  (C-20) 
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Collecting  the  results  for  ^ , and  Jle  from  Eqs.  (C-34),  (C-36),  and 

CC- 37)  according  to  Eq.  (C-28)  leads  to  the  result  =~Dr,  cos  y„  , and  so 

the  solution  for  the  flow  field  of  a pressure  dipole  also  conserves  angular 
momentum. 


3.  Rotor  Blade  Thickness 


The  flow  field  produced  by  the  thickness  of  a rotor  has  been  ob- 
tained by  the  superposition  of  source  solutions.  For  closed  blade  surfaces, 
the  net  source  strength  must  vanish.  Hence,  conservation  of  mass  requires 
that  there  be  no  net  introduction  of  mass  into  the  control  volume.  It  should 
be  mentioned  that  blade  shapes  having  flat-faced  trailing  edges  are  excluded 
from  consideration.  Blades  having  nonzero  thickness  at  the  trailing  dge  would 
require  infinite  source  strength  because  of  the  discontinuity  in  surface 
s lope . 


In  order  to  separate  the  rotor  flow  field  into  thickness  and  leading 
contributions,  it  has  been  assumed  that  the  rotor  is  unloaded  in  the  thick- 
ness problem.  In  Section  IV-C-3,  it  was  shown  tiiat  the  flow  field  solution 
displays  the  correct  behavior  of  zero  pressure  difference  across  the  blade 
surfaces.  Consequently,  there  should  be  no  net  force  or  torque  exerted  on 
the  fluid  by  the  rotor,  provided  thoie  is  no  blunt  trailing  edge  which  re- 
sults in  a base  drag.  However,  as  discussed  in  the  previous  paragraph  such 
a blade  shape  is  excluded  from  consideration  because  of  the  singularity  in 
the  source  strength  which  results. 

When  the  quantities'^,  7,,  , and  are  computed  from  the  mass  and 
momentum  balances  expressed  in  Eqs.  (C-l)  through  (C-3),  ail  ol  these 


quantities  must  vanish  if  the  flow  field  produced  by  rotor  thickness  is  to 
satisfy  the  global  conservation  laws.  The  evaluation  of  the  expressions  for 
‘M , and  0 i parallels  that  for  the  field  of  a mass  source.  The  flow  fieLd 
results  given  in  Section  IV-C-2  are  used  to  find  the  limiting  values  of  , 
1r£  , and  p , which  are  then  substituted  in  Eqs.  (C-5) . (C-S) , and  (C- 103 . 
The  results  obtained,  in  this  case,  are 
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These  integrals  are  rewritten  by  substituting  the  definition  of  Ura)  ) 
in  terms  of  the  blade  thickness  distribution  from  Eq.  (110). 
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1 4 

But  along  the  blade  surfaces,  UH  ~j-r-  t ( r.  , 5„)  can  be  written  as  U jj- 
and  so  the  z,  integration  which  occurs  in  both  ^ and  ^ becomes 
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Now,  for  closed  blade  surfaces,  t(r0,o)  = t ( r~c  , ) - 0 and  so  ^ - 0 . 

Hence,  from  F:q.  (C-38},1^  , ^ , and  Ji  all  vanish  as  required. 

4 . Rotor  Blade  Loading 

The  results  obtained  for  the  flow  field  of  a lifting  rotor  can  be 
showT.  to  satisfy  the  global  conservation  laws  by  following  much  the  same  steps 
as  for  the  dipc le  field.  Again  there  should  be  no  net  introduction  of  mass 
into  the  control  volume,  i.e.,  'frl-0  . The  force  which  each  rotor  blade  exerts 
on  the  fluid  is  given  by  Eq.  (138).  The  axial  component  of  the  force  of  the 
whole  rotor  on  the  fluid  is 

y*  - J J *?(''•  z>rdrdz  (c-42) 

where  Eq.  (30b)  has  been  used  for  the  unit  vector  a and  Eq.  (108)  has  been 
used  to  project  an  element  of  the  blade  surface  area  on  the  r , z plane. 
Similarly,  if  the  # -component  of  the  blade  force  at  each  radius  is  determined, 
the  axial  component  of  the  torque  exerted  by  the  rotor  is  found  to  be 

rc*  rrr 

-7,  - -8  j I A-pCr,  z ) r'oL  raz  (C-43) 
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In  order  to  evaluate  the  expressions  for  th,  ?t  , and  72  derived 
from  conservation  of  mass  and  momentum,  Eqs.  (C-l)  to  C-3),  the  rotor  flow 
field  variables  appearing  in  these  integrals  are  written  in  the  same  way  as  they 
were  for  the  dipole  field.  First  of  all,  the  axial  component  of  the  mass  flux 
is  written  in  the  form 
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(C-44) 


where  the  expressions  for  p am!  I are  given  by  Eqs.  (134)  and  (144)  respectively . 
The  integrand  in  the  relation  for  “fr? g Lven  in  Eq.  (C-l)  contains  the  limit  of 
this  quantity  for  £ « , and  so  the  limiting  values  of  p and  H / J £ are 

needed . 
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where  &p(ft)  is  defined  by  Eq.  (135). 


By  direct  analogy  with  the  dipole  results,  the  expression  for  is 
rewritten  as  the  sum  of  three  integrals  corresponding  to  integration  over 
each  of  the  three  terms  on  the  right-hand  side  of  Eq.  (C-44).  These  integrals 
are  labeled  , and  Jf  , respectively.  The  first  of  these  is 
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where  Eqs.  (135)  and  (C-45a)  have  been  used.  The  second  integral  vanishes  as 
a consequence  of  the  integration  over  $ , i.e., 


^4 


l 


rdddr 

-+0O 


* 0 


(C-47) 


The  third  integral  is 
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where  again  we  have  used  the  tact  that  at  fixed  Z , dO'd/?-  If  these  three 
integrals  are  added  together  then  the  result  ‘yn-o  is  obtained,  confirming 
that  the  rotor  flow  field  solution  conserves  mass. 

As  with  the  dipole  field,  the  expression  for  the  axial  component 
of  the  blade  force  computed  from  the  momentum  balance,  F.q.  (C-2),  can  be  ex- 
pressed in  terms  of  these  same  three  integrals.  If  Eqs.  (C-46]  to  (C-48}  are 
combined  according  to  Eq.  (C-27),  then  the  expression  obtained  for  the  axial 
component  of  the  rotor  force  agrees  with  the  result  of  direct  integration  of 
the  blade  forces  given  in  Eq.  (C-42) . Hence,  our  results  for  the  flow  field 
produced  by  rotor  loading  also  obey  global  momentum  conservation. 


The  expression  for  the  axial  component  of  the  torque  on  the  fluid 
obtained  from  conservation  of  angular  momentum,  Eq.  (C-3) , can  be  written  as 
in  Eq.  (C-28)  except  that  now  the  integrals  defined  by  Eqs.  (C-29)  to  (C-31) 
contain  the  rotor  rather  than  the  dipole  flow-field  perturbations.  When  the  ex- 
pression for  vr  given  in  Eq.  (146]  is  used  in  F.q.  (C-29),  the  integral 
becomes 
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This  result  is  obtained  by  carrying  out  the  integrations  over  9,2  • ana  r 
in  that  order.  The  first  two  are  straight  forward;  the  integration  over  r uses 
the  same  technique  required  to  evaluate  the  corresponding  integral  in  the 
source  and  dipole  cases. 


The  evaluation  of  J,  in  Eq.  (030)  requires  tne  9 component  of 
the  perturbation  velocity.  This  component  of  the  rotor  velocity  field  can 
be  obtained  from  the  normal  and  streamline  components  using  Eq.  (101a). 
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Using  the  limiting  form  of  31/3%  from  Eq.  (C-4Sb),  and  noting  that  the 
integration  of  this  quantity  over  a vanishes,  wc  obtain 
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Finally,  the  integrand  in  is  just  u>r  times  the  mass  flux  given  in  Eq.  (C-44). 
Tne  integrations  over  r and  6>  are  essentially  the  same  as  those  done  above 
for  *Yn , and  the  result  is 


e <*>* 5 rT*r*  j ^2B  _ 

“ IT  ~~  4 U*  d a CC-53) 

Combining  Jt  , J1  , and  according  to  Eq.  (C-28)  gives  the  same  result  for 
as  Eq.  (C-43) . Therefore,  the  solution  for  the  flow  field  of  a lifting 
rotor  satisfies  the  torque  and  angular  momentum  balance.  This  completes  the 
demonstration  that  the  flow  field  results  of  Section  IV  obey  the  global 
conservation  laws. 


137 


8&K  CONDENSER 
/ MICROPHONE 


ACOUSTIC  PROBE 
CALIBRATOR 


/ 


B&K  1/2  INCH 
MICROPHONE 

(position  m. 

PROBE  RESPONSE! 


TO  AUDIO  OSCILLATOR  OR 
PINK  NOISE  GENERATOR 


NOTE:  WHEN  MICROPHONE  IS  IN  POSITION  #1,  A DUMMY  MICROPHONE  IS  USED  IN 
POSITION  02  AND  VICE  VERSA 
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Figure  4 CALIBRATION  OF  MICROPHONE  IN  ACOUSTIC  PROBE 


ROTOR  SET  NO  1,  ROTOR  STAGGER  ANGLE  ’ 40.0  deg 
STATOR  SET  NO.  1.  STATOR  STAGGER  ANGLE  = 28.2  deg 
ROTOR  RPM  = 1435,  AUDIO  OSCILLATOR  FREQUENCY  = 4380  Hz 
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Figure  b OSC  LLOSCOPE  RECORDS  OF  NOISE  SIGNAL 
FROM  ROTOR  STATOR  INTERACTION 
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Figure  8 EXPERIMENTAL  OUTER  WALL  SOUND  PRESSURE  LEVELS  FROM  ROTOR-STATOR 
INTERACTION,  STATOR  STAGGER  ANGLE  = 28.2  deg 
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Figure  9 BLADE  TIP  PASSAGE  PRESSURE  SIGNALS  ON  ISOLATED  ROTOR 
AT  VARIOUS  ROTOR  SPEEDS,  CHORDWISE  LOCATION  = 13%, 
AXIAL  VELOCITY.  U0  = 60  ft/sec 
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Figuro  9 (Cont.)  BLADE  TIP  PASSAGE  PRESSURE  SIGNALS  ON  ISOLATED  ROTOR  AT 
VARIOUS  ROTOR  SPEEDS,  CHORDWISE  LOCATION  “ 18%,  AXIAL 
VELOCITY,  UQ  = 60  ft/sec 
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Figuis  10  BLADE  TIP  PASSAGE  PRESSURE  SIGNALS  AT  VARIOUS  GHORDWISE 
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jgurfc  IE  RMS  WALL  PRESSURE  AT  FOURTH  HARMONIC  OF  BLADE  PASSAGE  FREOUENCY 
VS  ROTOR  PPM  MEAN  STATOR  STAGGER  ANGLE  = 37.2  DEG 


on 


IKO  1200  noc  1400 

ROTOR  RPM 

Figure  16  RMS  tVAl PRESSURE  AT  FOURTH  HARMONIC  OF  BLADE  PASSAGE  FREQUENCY 
VS  RO'l'CM  P.PM;  MEAf-.'  STATOR  STAGGER  ANGI  E - 28.2  DEG 


Figure  18  THEORETICAL  SPL.  AMD  TOTAL  RADIATED  POWER  AT  RLAD7  PASSAGE 
FREQUENCY  VERSUS  ROTOR  RPM;  MEAN  S’  ATOR  STaGQER 
ANGLE  “ S? .7.  DEG 
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